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RESEARCH  ACCOMPLISHMENTS 


During  the  grant  period  worthwhile  contributions  have  been  made  on  the 
qualitative  and  quantitative  understanding  of  flames.  These  contributions  can 
be  approximately  grouped  into  four  categories:  (1)  Kinetic  structure  of 
laminar  flames;  (2)  Aerodynamic  structure  of  laminar  flames;  (3)  Structure  of 
turbulent  premixed  flames;  (4)  Soot  formation  in  laminar  diffusion  flames. 
The  accomplishments  are  summarized  in  the  following,  while  the  details  can  be 
found  in  the  published  journal  papers. 

1.  Kinetic  Structure  of  Laminar  Flames 
The  emphasis  in  this  area  is  in  understanding  and  quantifying  the 
detailed  and  semi -global  kinetics  of  fuel/oxidizer  systems  undergoing 
combustion. 

1.1  Asymptotic  Simulation  of  Hydrocarbon  Oxidation  Kinetics 

In  this  study  we  have  used  large  activation  energy  matched  asymptotic 
analysis  to  resolve  the  detailed  reaction  process  of  paraffin/air  oxidation  as 
represented  by  the  four-step  semi-global  reaction  mechanism  determined  from 
Princeton's  flow  reactor.  Four  major  regimes  respectively  representing  fuel 
pyrolysis,  hydrogen  and  carbon  monoxide  production,  hydrogen  oxidation,  and 
carbon  monoxide  oxidation  were  identified.  Embedded  within  these  four  major 
regimes  are  fifteen  sub-regimes.  Unlike  previous  asymptotic  studies  of 
reaction  systems  with  complex  chemistry,  transiti  us  between  these  sub-regimes 
are  not  governed  by  the  different  orders  of  activation  energies  but  rather  by 
the  creation/destruction  and  thereby  availability  of  the  various 
intermediates.  Explicit  expressions  were  derived  for  the  characteristic  time 
scales  associated  with  each  sub-regime.  Analytical  results  agree  well  with 
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the  experimental  flow  reactor  data. 

This  work  is  reported  in  Publication  No.  1. 

1.2  Nonadiabatic  Flame  Propagation  in  Dissociation  Equilibrium 

An  intriguing  question  in  the  fundamental  study  of  the  structure  and 
extinction  of  flames  is  the  influence  of  product  dissociation  on  the 
extinguishability  of  flames.  Recently  Williams  and  Peters  have  analyzed  the 
diffusion  flame  structure  and  found  that  product  dissociation  facilitates 
flame  extinction.  In  the  present  study  we  have  analyzed  the  premixed  flame 
structure  allowing  for  product  dissociation.  Analytical  solutions  have  been 
obtained  with  an  explicit  expression  for  the  burning  rate  eigenvalue  as  a 
function  of  the  dissociation  and  heat  loss  parameters.  These  reveal  the 
interesting  results  that  dissociation  alone  cannot  cause  extinction  of  an 
adiabatic  flame.  In  retrospect  this  is  reasonable  because  the  effect  of 
dissociation  is  equivalent  to  reducing  the  heat  of  combustion,  which  lowers 
the  flame  temperature  but  does  not  alter  the  basic  structure  of  the  laminar 
flame  propagation  which  does  not  exhibit  any  extinction  behavior. 

Analysis  has  also  been  performed  allowing  for  volumetric  heat  loss.  In 
this  case  it  was  found  that  extinction  is  facilitated  in  the  simultaneous 
presence  of  dissociation  and  heat  loss. 

This  work  is  reported  in  Publication  No.  2. 

1.3  Kinetic  and  Thermal  Structure  of  Diffusion  Flames 

We  have  analyzed  the  structure  of  a  diffusion  flame  supported  by  a 
representative  branching- termination  reaction  scheme 

F  +  Rx  -+  2R2 
0  +  R2  -  2Rj 
Rt  +  R2  +  M  -  2P  +  M 
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where  F  is  the  fuel,  0  the  oxidizer,  M  a  third  body,  and  Rj  and  R2  are 
regenerated  radicals.  The  analysis  yields  three  types  of  flame  behavior, 
depending  on  the  efficiency  of  the  recombination  reaction.  Of  particular 
interest  is  the  fast  recombination  regime  in  which  branching  and 
recombination  take  place  in  the  same,  thin  reaction  zone.  An  explicit 
extinction  criterion  has  been  derived  which  respectively  specialize  to  a 
kinetic  extinction  limit  and  a  thermal  extinction  limit;  in  the  latter  case  we 
recover  Linan's  extinction  criterion.  This  work  is  reported  in  Publication 
No.  3. 

1.4  Experimental  and  Numerical  Determination  of  Laminar  Flame  Speeds  of 

Methane/Air  Mixtures  under  Reduced  and  Elevated  Pressures 

The  laminar  flame  speed  S°  is  an  important  physico-chemical  parameter  of 

.Lj 

a  combustible  mixture  because  it  contains  the  basic  information  regarding  its 
diffusivity,  exothermicity  and  reactivity.  However,  in  spite  of  the  extensive 
efforts  expended  to  accurately  determine  their  values,  wide  systematic  spreads 
in  the  reported  experimental  data  still  exist  even  though  in  many  instances 
the  experiments  appear  to  have  been  carefully  executed.  Recently,  it  has  been 
suggested  that  these  systematic  discrepancies  are  likely  caused  by  the  coupled 
effects  of  flame  stretch  and  preferential  diffusion.  Specifically,  it  has 
been  demonstrated  both  theoretically  and  experimentally  that  the  flame 
response  can  be  qualitatively  reversed  when  the  flame  stretch  changes  from 
positive  to  negative  (e.g.  expanding  spherical  flame  versus  the  Bunsen  flame), 
and  when  the  mixture's  effective  Lewis  number  crosses  a  critical  value 
typically  around  unity  (e.g.  lean  methane/air  and  rich  propane/air  mixtures 
versus  rich  methane/air  and  lean  propane/air  mixtures). 

A  particularly  serious  implication  of  the  stretch- induced  flame  response 
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is  the  potential  falsification  of  the  kinetic  information  determined  or 
validated  through  comparison  between  the  numerically-calculated  and 
experimentally-determined  results.  In  view  of  this  concern,  we  have  recently 
proposed  a  stagnation- flow  based  methodology  through  which  stretch  effects  can 
be  systematically  subtracted  out  such  the  S°  can  be  unambiguously  determined. 

In  the  present  investigation  we  have  accurately  determined  the  S°  of 

methane/air  mixtures  under  reduced  and  elevated  pressures.  While  these  data 
are  of  practical  interest  in  its  own  right  in  terms  of  high-pressure 
combustion,  we  note  in  addition  that  these  data  also  carry  significant  kinetic 
information  because  pressure  not  only  influences  the  frequency  of  molecular 
collision  but  also  differentiates  the  relative  efficiencies  of  the  two-body 
branching  reaction  versus  the  three-body  termination  reaction.  Thus  they 
provide  additional  constraints  for  the  validation  of  the  kinetic  schemes. 

The  experimental  data  have  been  compared  with  the  numerically-calculated 
values  obtained  by  assuming  different  kinetic  mechanisms.  Consequently,  the 
Cx  mechanism  of  Miller  and  C2  mechanism  of  Warnatz  both  accurately  reproduce 
our  experimental  data  over  the  extensive  concentration  and  pressure  ranges 
tested.  The  relative  superiority  of  these  two  schemes  cannot  be  further 
resolved. 

This  work  is  reported  in  Publication  No.  4. 

1.5  Laminar  Flame  Speeds  of  Methane/" Air"  Mixtures  as  a  Function  of 

Stoichiometry ,  Pressure  and  Flame  Temperature 

The  work  of  1.4  has  been  extended  to  the  determination  of  laminar  flame 
speeds  of  methane/"air"  mixtures  over  the  stoichiometric  range  from  very  lean 
to  very  rich,  the  pressure  range  from  0.25  to  2  atm,  and  flame  temperature 
range  from  1,550  and  2,250  K;  independent  variation  in  the  flame  temperature 
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is  achieved  by  substituting  nitrogen  in  the  air  by  an  equal  amount  of  argon  or 
carbon  dioxide.  These  data  are  expected  to  be  useful  for  the  partial 
validation  of  proposed  kinetic  mechanisms.  In  the  present  study  numerical 
simulation  of  the  experimental  flame  speeds  has  been  conducted  by  using  a  Cj 
mechanism  and  a  full  C2  mechanism.  The  calculated  results  agree  well  with  the 
experimental  data,  for  both  the  and  C2  mechanisms,  except  for  the  very  rich 
mixtures  for  which  there  is  substantial  over-prediction  by  the  C2  mechanism. 
Sensitivity  analyses  have  also  been  performed  where  appropriate  for  enhanced 
insight  into  the  controlling  elementary  reactions. 

This  work  is  reported  in  Publication  No.  5. 

1.6  Experiments  on  Extinction  of  Diffusion  Flames 

In  this  study  we  experimentally  investigate  the  theoretical  results  of 
1.3.  The  theory  predicts  that,  because  of  the  simultaneous  presence  of 
branching  and  termination  reactions,  the  density-weighted  extinction  strain 
rate  of  a  diffusion  flame  increases  only  linearly  with  pressure,  while  the 
conventional  formulation  assuming  only  a  one-step  second-order  overall 
reaction  predicts  a  quadratic  dependence.  The  experiment  was  conducted  by 
using  the  counterflow  diffusion  flame  configuration.  Experimental  results 
indeed  reveal  a  linear  dependence. 

This  work  is  reported  in  Publication  No.  6. 

2.  Aerodynamic  Structure  of  Laminar  Flames 
The  emphases  here  are  on  the  influence  of  stretch  and  preferential 
diffusion  on  flames. 

2.1  Propagation  and  Extinction  of  Stretched  Premixed  Flames 

Practically  all  realistic  flames  are  subjected  to  aerodynamic  stretching, 
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manifest  through  flow  nonuniformity,  f i  amo  fiurv’a<~'lga-i - ar - SXmmm  if  1  <w 

unsteadiness.  These  processes  can  profoundly  affect  the  flame  from  those 
prescribed  by  the  classical  one -dimensional  planar  flame  propagation  theory, 
which  is  frequently  the  one  used  in  modeling  practical  combustion  phenomena. 

The  two  most  important  parameters  characterizing  stretched  flames  are  the 
stretched  flame  speeds  and  the  states  of  extinction  when  certain  critical 
stretch  rates  are  exceeded.  Unfortunately,  these  data  have  not  been  obtained. 
Therefore,  the  first  contribution  of  the  present  study  is  an  accurate 
determination  of  these  values  as  functions  of  stoichiometry  and  stretch  for 
both  methane/air  and  propane/air  mixtures.  These  results  further  show  that 
the  flame  speed  is  finite  at  the  state  of  extinction  and  that  the  extinction 
limits  for  nearly  adiabatic,  stretchless,  and  planar  flames  appear  to 
correspond  to  the  flammability  limits  of  the  respective  mixtures.  It  is 
suggested  that  kinetic  termination  could  be  the  dominant  extinction  mechanism 
at  the  flammability  limits. 

This  work  is  reported  in  Publication  No.  7. 

2.2  Interaction  Between  Laminar  Flames 

Modeling  of  turbulent  combustion  can  be  frequently  facilitated  by 
approximating  the  turbulent  flame  as  an  ensemble  of  laminar  flamelets.  It  is 
clear  that,  contrary  to  previous  modeling  efforts  in  which  such  flamelets  are 
assumed  to  be  isolated  ones,  they  are  actually  in  close  proximity  and 
therefore  should  interact  strongly  with  each  other. 

In  a  previous  study  we  have  experimentally  mapped  the  extinction  limits 
of  two  interacting  laminar  premixed  flames.  Subsequently  a  theory  has  been 
formulated  for  the  extinction  of  this  interacting  system.  The  results  show 
good  agreement  with  the  experimental  data  for  lean  methane/air  mixtures  but 
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In  the  present  study  we  first  argued  that  such  a  discrepancy  is  a 
consequence  of  the  unity  Lewis  number  assumption  in  the  theory.  The  theory 
was  subsequently  extended  to  allow  deviation  of  the  Lewis  numbers  from  unity. 
Satisfactory  agreement  with  all  experimental  data  was  then  obtained. 

This  work  is  reported  in  Publication  No.  8. 

2.3  Dynamics  of  Stretched  Flames 

This  is  an  invited  review  paper  presented  at  the  22nd  Combustion 
Symposium.  Here  recent  advances  in  the  understanding  of  the  structure, 
propagation  and  extinction  of  laminar  flames  under  the  influence  of  stretch, 
as  manifested  by  the  existence  of  flame  curvature,  flow  nonuniformity,  and 
flame  motion,  are  reviewed.  The  emphasis  is  on  premixed  flames  because  of  the 
richness  and  subtlety  of  the  phenomena  involved.  The  review  distinguishes  the 
influences  of  the  tangential  and  normal  velocity  gradients  on  the  flame 
response,  both  at  the  hydrodynamic  scale  and  within  the  flame  structure,  and 
emphasizes  the  importance  of  the  preferential  diffusion  nature  of  heat  and 
mass  transport,  as  well  as  the  extent  to  which  the  flame  can  freely  adjust  its 
location  in  response  to  stretch  in  order  to  achieve  complete  reaction.  It  is 
then  demonstrated  that  stretch  has  only  minimal  effect  on  an  adiabatic, 
unrestrained,  diffusionally-neutral  flame  with  complete  reaction  in  that  the 
temperature,  propagation  rate,  and  thickness  of  the  flame  are  invariant  to 
stretch,  and  that  stretch  alone  cannot  extinguish  such  a  flame.  In  the 
presence  of  preferential  diffusion  and/or  when  the  flame  movement  is 
restrained,  the  response  of  the  flame  to  stretch  becomes  more  sensitive  and 
extinction  is  also  possible.  The  concept  of  flame  stretch  is  applied  to 
interpret  such  practical  flame  phenomena  as  flame  stabilization  and  flame - 
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front  instability,  determination  of  laminar  flame  speeds  and  flammability 
limits,  concentration  and  temperature  modifications  in  flame  chemistry,  and 
modeling  of  turbulent  flames.  The  properties  of  stretched  diffusion  flames 
are  then  briefly  discussed.  The  review  closes  with  suggestions  for  further 
research. 

This  paper  is  designated  as  Publication  No.  9. 

3.  Structure  of  Turbulent  Preaixed  Flames 

3.1  Structure  and  Propagation  of  Turbulent  Premixed  Flames 

Previous  experimental  studies  on  turbulent  premixed  flames  have  used 
either  the  Bunsen  flame  configuration  or  the  V- flame  stabilized  by  a  rod. 
Data  obtained  from  these  two  flame  configurations,  however,  are  complicated  by 
the  uncertainty  in  the  flame  orientation  and  the  eddies  generated  in  the  flame 
stabilization  region.  In  the  present  investigation  we  have  adopted  the 
stagnation  flame  as  our  model  configuration.  The  advantages  here  are  the 
planar  nature  of  the  flame  which  is  also  normal  to  the  bulk  flow,  and  the 
absence  of  fluctuations  generated  by  the  holding  region. 

Using  this  flame  configuration,  experiments  have  been  conducted  covering 
equivalence  ratios  of  0.67  to  1.16,  incident  turbulence  intensities  of  0.1  to 
0.7  m/s,  and  global  stretch  rates  of  100  to  230/sec.  Results  on  turbulence 
statistics  show  that  in  crossing  the  flame  the  increase  in  fluctuation  is  more 
pronounced  in  the  normal  component,  while  the  Reynolds  stress  remains  almost 
zero,  and  that  the  velocity  joint  probability  density  function  is  bimodal  for 
rich  flames  and  mono-modal  for  lean  flames.  Turbulent  flame  speeds  determined 
as  functions  of  the  equivalence  ratio  and  turbulence  intensity  agree  well  with 
those  in  the  literature. 
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This  work  is  reported  in  Publication  No.  10. 

3.2  Velocity  and  Scalar  Fields  of  Turbulent  Premixed  Flames  in  Stagnation  Flow 
Detailed  experimental  measurements  of  the  scalar  and  velocity  statistics 
of  premixed  methane/air  flames  stabilized  by  a  stagnation  plate  are  reported. 
Conditioned  and  unconditioned  velocity  of  two  components  and  the  reaction 
progress  variable  are  measured  by  using  a  two-component  laser  Doppler 
velocimetry  technique  and  Mie  scattering  technique,  respectively. 
Experimental  conditions  cover  equivalence  ratios  of  0.9  and  1.0,  incident 
turbulence  intensities  of  0.3  tc  0.45  m/s,  and  global  stretch  rates  of  100  to 
150  sec'-*-.  The  expei  imental  results  are  analyzed  in  the  context  of  the  Bray- 
Moss-Libby  flamelet  model  of  these  flames.  The  results  indicate  that  there  is 
no  turbulence  production  within  the  turbulent  flame  brush  and  the  second  and 
third  order  turbulent  transport  terms  are  reduced  to  functions  of  the 
difference  between  the  conditioned  mean  velocities  of  the  axial  and  transverse 
components.  The  result  of  normalization  of  these  relative  velocities  by  the 
respective  velocity  increase  across  laminar  flames  suggest  that  the  mean 
unconditioned  velocity  profiles  are  self -similar . 

This  work  is  reported  in  Publication  No.  11. 

4.  Soot  Formation  in  Laminar  Diffusion  Flames 
In  this  series  of  studies  the  effects  of  dilution,  preferential 
diffusion,  inert  addition,  and  aerodynamic  straining  on  soot  formation  in 
laminar  diffusion  flames  were  quantified. 

4.1  Dilution  and  Temperature  Effects  on  Inert  Addition  on  Soot  Formation  in  a 
Counterflow  Diffusion  Flame 

The  isolated  effects  of  fuel  dilution  and  flame  temperature  variation  on 
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soot  formation  have  been  investigated  experimentally  for  counterflow  diffusion 
flames.  The  methodology  of  isolation  through  temperature  adjustment  involves 
changing  the  concentration  of  the  fuel  by  diluting  it  with  nitrogen  and  then 
increasing  the  maximum  temperature  of  the  diluted  flame  back  to  that  of  the 
undiluted  flame  by  replacing  a  portion  of  the  nitrogen  in  the  oxidizer  flow 
with  an  equal  portion  of  argon.  Soot  quantities  are  determined  by  using 
light-scattering  and  extinction  techniques,  flame  temperatures  by  thermocouple 
measurements,  and  velocities  by  LDV  mapping.  Results  on  ethylene  show  that 
both  temperature  and  dilution  contribute  to  a  reduction  in  the  soot  volume 
fraction  and  dilution  accounts  for  a  substantial  amount  of  this  reduction.  It 
is  further  found  that  the  soot  formation  rate  and  specific  surface  area  are 
linearly  dependent  on  che  initial  fuel  concentration,  while  the  specific 
surface  growth  rate  is  insensitive  to  its  variation.  Finally,  the  soot 
formation  process  is  insensitive  to  a  wide  range  of  fuel  ejection  rates  for 
the  present  counterflow  diffusion  flame  established  in  the  forward  stagnation 
region  of  a  porous  cylinder. 

The  work  is  reported  in  Publication  No.  12. 

4.2  Preferential  Diffusion  and  Concentration  Modification  in  Sooting 

Counterflow  Diffusion  Flames 

An  experimental  investigation  has  been  conducted  on  the  influence  of  the 
mobility  of  inert  additives  on  soot  formation  in  propane  and  ethylene 
counterflow  diffusion  flames.  Inerts  used  were  helium,  neon,  argon,  or 
krypton,  and  the  results  show  that  while  the  mobility  of  the  inert  has 
practically  no  effect  when  a  small  amount  of  inert  is  added  to  the  oxidizer 
side,  the  influence  is  significant  when  added  to  the  fuel  side  in  that 
krypton,  being  the  least  mobile  inert,  yields  the  greatest  soot  loading  while 
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helium,  being  the  most  mobile,  yields  the  least.  By  relating  the  spatially- 
resolved  soot  volume  fractions  to  the  corresponding  profiles  of  temperature, 
velocity  and  species  concentrations,  it  is  demonstrated  that  this  influence  on 
soot  loading  is  likely  caused  by  concentration  modifications  of  the  fuel  and 
the  soot  precursors  due  to  the  different  mobilities  of  the  inert  additives. 

This  work  is  reported  in  Publication  No.  13. 

4 .3  Experiments  on  the  Sooting  Limits  of  Aerodynamics lly- Strained  Diffusion 

Flames 

An  experimental  study  has  been  performed  with  axisymmetric  counterflow 
diffusion  flames  to  investigate  the  influence  of  aerodynamic  straining  on  the 
relevant  sooting  limits  of  the  lower  alkanes.  The  limits  are  defined  by  the 
critical  strain  rate  at  which  either  soot  luminosity,  soot  particle 
scattering,  or  fluorescence  is  negligible  compared  to  the  appropriate 
background  signal.  These  critical  strain  rates  are  found  to  be  greatest  for 
the  sooting  limit  based  on  the  fluorescence  signal,  and  least  for  that  based 
on  luminosity.  The  fluorescence  signal,  if  attributed  to  polycyclic  aromatic 
precursors,  yields  a  limit  that  can  be  interpreted  as  the  extinction  of  soot 
precursors  and  is  suggested  to  be  a  possible  limit  for  identifying  a 
completely  nonsooting  flame  condition.  The  separate  effects  of  flame 
temperature  and  fuel  concentration  on  the  critical  strain  rates  for  soot 
extinction  have  also  been  studied.  The  results  are  indicative  of  how 
temperature  and  concentration  influence  the  soot  particle  inception  process 
and  they  show  that  both  are  potentially  important  parameters.  The  critical 
strain  rates  display  an  Arrhenius  temperature  dependence  and  this  dependence 
is  similar  for  all  alkanes  considered. 

This  work  is  reported  in  Publication  No.  14. 
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by  B.H.  Chao  and  C.K.  Law,  Twenty-First  Combustion  Symposium.  The 
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4.  "Laminar  Flame  Speeds  of  Me thane /Air  Mixtures  Under  Reduced  and  Elevated 
Pressures,"  by  F.  Egolfopoulos ,  P.  Cho,  and  C.K.  Law,  Combustion  and 
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PRESENTATIONS  AND  SIGNIFICANT  INTERACTIONS 


During  the  reporting  period  the  Principal  Investigator  gave  the  following 
invited  and  contributed  talks. 


Invited  Talks : 

1.  "Combustion  of  Multicomponent  Droplets,"  Ballistic  Reserach  Laboratory, 
Aberdeen  Proving  Ground,  Aberdeen,  MD,  Mar. 5  ,  1985. 

2.  "Dynamics  of  Stretched  Premixed  Flames,"  Department  of  Mechanical 
Engineering,  McGill  University,  Montreal,  Quebec,  Canada,  Mar.  6,  1985. 

3.  "Behavior  of  Non-Dilute  Sprays,"  AFOSR/ARO  Workshop  on  Non-Dilute  Sprays, 
Sandia  Laboratories,  Livermore,  CA,  Mar.  20-21,  1985.  **Lead  Talk** 

4.  "Recent  Advances  in  Droplet  Combustion,"  Department  of  Mechanical 

Engineering,  Cornell  University,  Ithaca,  NY,  Apr.  30,  1985. 

5.  "Heat  and  Mass  Transfer  in  Combustion,"  Department  of  Applied  Mechanics 
and  Engineering  Sciences,  University  of  California,  San  Diego,  CA,  Oct. 
11,  1985. 

6.  "Dynamics  of  Stretched  Flames,"  Department  of  Mechanical  Engineering, 

University  of  California,  Davis,  CA,  Jan.  30,  1986. 

7.  "Recent  Advances  in  Droplet  Combustion,"  Department  of  Mechanical 

Engineering,  Louisianna  State  University,  Baton  Rogue,  LA,  Feb.  18,  1986. 

8.  "Dynamics  of  Stretched  Flames,"  Department  of  Mechanical  Engineering, 

University  of  Oklahoma,  Norman,  OK,  Feb.  19,  1986. 

9.  "Recent  Advances  in  Droplet  Combustion,"  Department  of  Mechanical 

Engineering,  Texas  A&M  University,  College  Station,  TX,  Feb.  20,  1986. 

10.  "Recent  Advances  in  Droplet  Combustion,"  Department  of  Mechanical 

Engineering,  University  of  Texas,  Arlington,  TX,  Feb.  21,  1986. 

11.  "Recent  Advances  in  Droplet  Combustion,"  Department  of  Mechanical 

Engineering  and  Aerospace  Engineering,  Princeton  University,  Princeton, 
NJ,  April  1,  1986. 

12.  "Dynamics  of  Stretched  Flames,"  Department  of  Mechanical  Engineering, 

University  of  Colorado,  Boulder,  CO,  April  14,  1986. 

13.  "Combustion  of  High  Energy  Propellants,"  Martin  Marietta  Research  Lab., 
Baltimore,  MD,  July  28,  1986. 

14.  "Recent  Advances  in  Droplet  Combustion,"  Department  of  Mechanical 

Engineering,  Washington  State  University,  Pullman,  WA,  Sept.  10,  1986. 
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15.  "Recent  Studies  on  the  Structure  of  Laminar  Premixed  Flames,"  Department 
of  Mechanical  Engineering,  University  of  California,  Berkeley,  CA,  Oct. 
22,  1986. 

16.  Recent  Advances  in  Multicomponent  and  Propellant  Droplet  Vaporization  and 
Combustion,"  ASME  Winter  Annual  Meeting,  Anaheim,  CA,  Dec.  8-12,  1986. 
**Invited  Talk** 

17.  "Recent  Advances  in  Multicomponent  Droplet  Combustion,"  Department  of 
Mechanical  Engineering,  University  of  California,  Davis,  CA,  April  9, 
1987. 

18.  "Incineration  of  Chlorinated  Hydrocarbons,"  Workshop  on  Research  in 
Engineering  and  System  Analysis  for  the  Control  of  Toxic  Substances," 
University  of  California,  Los  Angeles,  CA,  Oct.  21-22,  1987. 

19.  "Structure,  Propagation,  and  Stability  of  Laminar  Premixed  Flames," 
Department  of  Mechanical  Engineering,  Stanford  University,  Palo  Alto,  CA, 
Nov.  4,  1987. 

20.  "Single  Droplet  Considerations  in  Spray  Combustion  Processes- 
Experimental  Overview,"  Workshop  on  Mass,  Momentum  and  Energy  Exchange  in 
Combusting  Sprays:  Droplet  Studies,  Sandia  National  Laboratories, 
Livermore,  CA,  March  28-29,  1988.  **Lead  Talk** 

21.  "Blending  Strategy  for  the  Incineration  of  Chlorinated  Hydrocarbons," 
Workshop  on  Incineration  Technologies  for  Municipal,  Toxic  and  Hospital 
Wastes,  National  Science  Foundation,  Washington,  DC,  April  18-19,  1988. 

22.  "Dynamics  of  Stretched  Flames,"  Twenty-Second  Symposium  (International) 
on  Combustion,  Seattle,  WA,  Aug.  15-19,  1988.  **Invited  Paper** 

23.  "Some  Recent  Advances  in  Droplet  Combustion,"  Third  International 

Colloquium  on  Drops  and  Bubbles,  Monterey,  CA,  Sept.  19-21,  1988. 

**Keynote  Paper** 

24.  "On  Flame  Extinction  and  Flammability  Limits  of  Combustibles,"  Department 
of  Mechanical  and  Aerospace  Engineering,  Rutgers  University,  New 
Brunswick,  NJ,  Nov.  2,  1988. 

25.  "Vaporization,  Combustion,  and  Collisional  Dynamics  of  Fuel  Droplets," 
Department  of  Mechanical  Engineering,  University  of  Pennsylvania, 
Philadelphia,  PA,  Feb.  9,  1989. 

26.  "Extinction  and  Flammability  of  Gaseous  Combustibles,"  Department  of 
Mechanical  Engineering,  Yale  University,  New  Haven,  CT,  Feb.  22,  1989. 
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Contributed  Talks 


1.  "Laminar  Flame  Speeds  of  Hydrocarbon/Air  Mixtures  with  Hydrogen 
Addition,"  Central  and  Western  States  Sections  of  the  Combustion 
Institute  Meeting,  San  Antonio,  TX,  Apr.  22-23,  1985. 

2.  "Droplet  Combustion  of  Organic  Azides,"  AIAA/SAE/ASME/ASEE  21st  Joint 
Propulsion  Conference,  Monterey,  CA,  July  8-10,  1985. 

3.  "Aerodynamic  and  Kinetic  Processes  in  Flames,"  AFOSR  Contractors' 
Meeting,  California  Institute  of  Technology,  Pasadena,  CA,  July  23-25, 

1985. 

4.  "Combustion  of  Liquid  Propellant  Droplets,"  JANNAF  Regenerative  Liquid 
Gun  Propellant  Spray  Combustion  Workshop,  Aberdeen  Proving  Ground,  MD, 
Aug.  26-27,  1985. 

5.  "Stability  and  Stabilization  of  Premixed  Flames,"  DOE  Third  Symposium  on 
Energy  Engineering  Sciences,  Pennsylvania  State  University,  University 
Park,  PA,  Oct.  8-10,  1985. 

6.  "Organic  Azides  as  Jet  Fuel  Additives,"  ONR  Workshop  on  Combustion 
Instability  in  Compact  Ramjets,  Johns  Hopkins  University,  Baltimore,  MD, 
Oct.  16-17,  1985. 

7.  "Propagation  and  Extinction  of  Stretched  Premixed  Flames,"  Fall  Technical 
Meeting  of  the  Western  States  Section  of  the  Combustion  Institute, 
University  of  California,  Davis,  CA,  Oct.  21-22,  1985. 

8.  "Combustion  and  Micro-Explosion  of  Water/Oil  Emulsions  in  High  Pressure 
Environments,"  Sixth  ARO  Combustion/Fuels  Workshop,  University  of 
California,  Davis,  CA,  Mar.  17-18,  1986. 

9.  "Droplet  Combustion  of  HAN-Based  Liquid  Propellants,"  Army  Workshop  on 
HAN-Based  Liquid  Propellants,  Aberdeen  Proving  Ground,  MD,  July  29,  1986. 

10.  "Liquid- Phase  Diffusional  Resistance  in  Multicomponent  Droplet 
Gasification,"  21st  Combustion  Symposium,  Munich,  Germany,  Aug.  3-8, 

1986. 

11.  "Time -Resolved  Gasification  and  Sooting  Characteristics  of  Droplets  of 
Alcohol/Oil  Blends  and  Water/Oil  Emulsions,"  21st  Combustion  Symposium, 
Munich,  Germany,  Aug.  3-8,  1986. 

12.  "Flame  Curvature  and  Preferential  Diffusion  in  the  Burning  Intensity  of 
Bunsen  Flames,"  21st  Combustion  Symposium,  Munich,  Germany,  Aug.  3-8, 
1986. 
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21st 


13.  "Propagation  and  Extinction  of  Stretched  Premixed  Flames," 
Combustion  Symposium,  Munich,  Germany,  Aug.  3-8,  1986. 

14.  "Vaporization,  Combustion,  and  Microexplosion  of  Droplets  of  Organic 
Azides,"  23rd  JANNAF  Combustion  Meeting,  NASA-Langley ,  Oct.  20-24,  1986. 

15.  "Propagation  and  Extinction  of  Stretched  Premixed  Flames,"  Sandia 
Turbulent  Flame  Workshop,  University  of  California,  Davis,  CA,  April  21- 
22,  1987. 

16.  "Aerodynamic  and  Kinetic  Processes  in  Flames,"  AF0SR/0NR  Contractor's 
Workshop,  Pennsylvania  State  University,  College  Station,  PA,  June  22-26, 
1987 . 

17.  "Some  Issues  Related  to  Stretched  Flames,"  US-Japan  Workshop  on 
Mathematical  Combustion,  Juneau,  Alaska,  Aug.  17-21,  1987. 

18.  "Kinetic  Effects  in  Flame  Analysis,"  US-Japan  Workshop  on  Mathematical 
Combustion,  Juneau,  Alaska,  Aug.  17-21,  1987. 

19.  "Droplet  Combustion  of  HAN-Based  Liquid  Propellants,"  Worksnop  on  Liquid 
Propellants,  Aberdeen  Proving  Ground,  Aberdeen,  MD,  Aug.  25-27,  1987. 

20.  "Droplet  Extinction  Studies,"  Workshop  on  Hazardous  Waste  Incineration, 
University  of  California,  Davis,  CA,  Sept.  14-15,  1987. 

21.  "Ignition  and  Combustion  of  liquid  Gun  Propellants,"  Seventh  Army 
Engine/Fuels  Workshop,  Detroit,  MI,  Oct.  13-14,  1987. 

22.  "Combustion  Studies  of  Energetic  Liquid  Materials,"  ONR  Contractor's 
Meeting,  University  fo  California,  Berkeley,  CA,  Nov.  17-19,  1987. 


18 


Appendix  A 


Twenty-Second  Symposium  (International)  on  Combustion/The  Combustion  Institute, 


19S8//5'/ 


EXPERIMENTAL  AND  NUMERICAL  DETERMINATION  OF  LAMINAR 
FLAME  SPEEDS  OF  METHANE/(Ar,  N2,  C02)-AIR  MIXTURES  AS 
FUNCTION  OF  STOICHIOMETRY,  PRESSURE,  AND  FLAME 

TEMPERATURE 


D  L.  ZHU,  F.  N  ECOLFOPOULOS  and  C.  K.  LAW 
Department  of  Mechanical  Enninrerine 
University  oj  California 
Davis.  California  95616 


By  using  the  countertlow  method,  the  laminar  flame  speeds  of  methane/(Ar.  N,.  C'Oj-uir 
mixtures  have  been  accurately  and  extensively  determined  over  the  stoichiometric  range  from 
very  lean  to  very  rich,  the  pressure  range  from  0.25  to  2  atm.  and  flame  temperature  range 
from  1,550  to  2,250  K;  independent  variation  in  the  flame  temperature  is  achieved  In  sub¬ 
stituting  nitrogen  in  the  air  by  an  equal  amount  of  argon  or  carbon  dioxide.  These  data  are 
expected  to  be  useful  for  the  partial  validation  of  proposed  kinetic  mechanisms.  In  the  pres¬ 
ent  study  numerical  simulation  of  the  experimental  flame  speeds  has  been  conducted  In 
using  a  C,  mechanism  and  a  full  C,  mechanism.  The  calculated  results  agree  well  with  the 
experimental  data,  for  both  the  C,  and  C\  mechanisms,  except  for  very  rich  mixtures  for 
which  there  is  substantial  overprediction  by  the  C.t  mechanism.  Sensitivity  analyses  have  also 
been  performed  where  appropriate  for  enhanced  insight  into  the  controlling  elementary  re¬ 
actions. 


Introduction 

Recent  studies  of  the  chemical  kinetics  of  com¬ 
plex  reactants  of  interest  to  combustion  have  been 
greatly  aided  by  the  ability  to  simulate  numerically 
the  progress  of  a  large  number  of  possible  reactions 
involving  many  chemical  species.12  Such  a  simu¬ 
lation  can  serve  two  purposes.  That  is,  if  the  ki¬ 
netics  of  the  elementary  steps  arc  accurately  known, 
then  the  simulation  allows  the  study  of  the  influ¬ 
ence  of  chemistry  on  the  combustion  phenomena  of 
interest  with  confidence.  If,  however,  part  of  the 
kinetic  mechanism  is  either  unknown  or  not  well 
established,  then  the  numerical  simulation  can  be 
used  to  either  identify  and/or  validate  the  needed 
kinetics  information  through  comparisons  with  re¬ 
sults  obtained  through  independent  experimenta¬ 
tion  It  is  clear  that  such  a  comparison  is  useful, 
and  meaningful,  only  if  the  experimental  situation 
closely  simulates  that  being  calculated,  and  if  the 
experimental  data  are  of  the  same  degree  ol  accu¬ 
racy,  and  the  same  level  of  detail,  as  the  calculated 
results. 

The  model  flame  configuration  that  is  perhaps 
most  suited  for  such  a  simulation,  for  the  study  of 
chemical  kinetics,  is  the  steady  propagation  of  the 
adiabatic,  one-dimensional  planar  premixed  flame 
in  the  doubly-infinite  domain.  Since  this  situation 
is  totally  independent  of  external  influences  such  as 


those  of  aerodynamics  and  heat  and  radical  ex¬ 
changes,  the  structure  and  the  propagation  rate, 
St,  of  the  flame  are  therefore  unique  properties  of 
the  combustible  mixture,  representing  its  reactivity- 
in  the  presence  of  diffusive  transport.  Thus  an 
agreement  between  the  calculated  and  measured 
flame  propagation  speed,  as  well  as  the  profiles  of 
temperature  and  species  concentrations,  should  be 
a  strong  indication  that  the  kinetic  mechanism  used 
in  the  calculation  is  a  sufficiently  accurate,  and 
thereby  adequate,  representation  of  the  complete 
mechanism. 

Such  a  requirement,  while  desirable,  is  at  pres¬ 
ent  not  achievable  owing  to  the  lack  of  accurate  and 
sufficient  experimental  information  on  the  flame 
structure.  On  the  other  hand,  extensive  data  on  the 
laminar  flame  propagation  speed  have  been  accu¬ 
mulated  over  the  years.  Therefore  much  of  the  fo¬ 
cus  on  the  (partial)  validation  of  kinetic  mechanisms 
has  been  on  the  comparison  between  calculated  and 
measured  flame  speeds.1  3  Obviously  failure  to 
achieve  agreement  for  even  such  a  global  flame 
property  is  a  clear  indication  that  the  mechanism 
requires  revision.  The  usefulness  of  such  a  com¬ 
parison  has  recently  been  further  enhanced  with  the 
development  of  the  counterflow  method  for  the  de¬ 
termination  of  laminar  flame  speeds.4  5  That  is,  ex¬ 
cept  for  the  flat-flame  burner  method,  all  other 
techniques  adopted  for  measuring  flame  speeds  uti- 
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lize  flames  which  suffer  aerodynamic  straining  man¬ 
ifested  through  either  flow  nonuniformity,  or  flame 
curvature,  or  flame  motion.  When  coupled  to  the 
inherent  preferential  diffusion  effects,  the  flame  be¬ 
havior  such  as  its  temperature  and  propagation  speed 
can  be  qualitatively  and  systematically  affected.6  ‘ 
This  is  probably  one  of  the  major  reasons  that  pre¬ 
viously  reported  flame  speeds,  determined  by  var¬ 
ious  investigators  with  great  care  but  by  using  dif¬ 
ferent  flame  configurations,  can  deviate  from  each 
other  by  substantial  amounts.  For  example,  values 
of  25  cm/s  to  -15  cm/s  have  been  reported  as  the 
laminar  flame  speed  of  stoichiometric  methane/air 
mixtures  at  1  atm.'s  ft  is  therefore  clear  that  such 
strain-induced  effects  will  have  to  be  subtracted  out 
from  the  experimentally-determined  flame  speed  SL 
to  yield  the  proper,  strain-free,  one-dimensional 
flame  speed  S£  before  comparison  with  the  nu¬ 
merical  results  is  valid.  Failure  to  do  so  can  lead 
to  falsification  of  the  kinetic  mechanism  "validated" 
through  comparison.  The  counterflow  method,  to 
be  discussed  in  the  next  section,  offers  such  a  sys¬ 
tematic  deduction.  We  mention  in  passing  that  the 
strainless  flat-flame  burner  method,  with  extrapo¬ 
lation  to  vanishing  heat  loss,9  still  cannot  totally 
simulate  the  model  laminar  flame  because  of  the 
potential  radical  loss  at  the  burner  surface. 

The  general  objective  of  the  present  study  is  to 
provide  experimental  data  of  high  fidelity  and  ac¬ 
curacy  for  the  laminar  flame  speeds  S'i  of  methane/ 
inert/oxvgen  mixtures  by  using  the  counterflow 
method.  These  results  will  then  be  compared  with 
the  numerically  calculated  values  by  assuming  dif¬ 
ferent  kinetic  mechanisms.  Because  of  the  global 
and  somewhat  limited  nature  of  the  flame  speed  as 
a  representation  of  the  detailed  kinetics,  we  shall 
introduce  as  much  kinetics  information  as  possible 
into  the  reported  values  of  the  flame  speed.  Spe¬ 
cifically.  wc  shall  extensively  determine  the  laminar 
flame  speed  not  only  as  a  function  of  the  mixture 
equivalence  ratio  <t>  and  pressure  p,  which  arc  the 
traditional  independent  parameters  used,  but  we 
shall  also  map  it  as  an  explicit  function  of  the  adi¬ 
abatic  flame  temperature  Taj.  which  has  not  been 
attempted  before.  Thus 

sy.  =  sy  (4).  P.  T„i). 

Since  T„|  of  a  methane/air  mixture  is  fixed  for  a 
given  <t>,  in  the  present  investigation  Ta(|  is  to  be 
independently  varied  by  substituting  N2  in  the  air 
by  equal  amounts  of  either  Ar  or  C02  such  that  the 
02  concentration  in  the  (02  +  inert)  mixture  is  fixed 
at  21  volume  percent.  Thus  the  inert  in  our  "air 
can  consist  of  either  N2,  or  (Ar  +  N2),  or  (C02  + 
N2). 

The  fact  that  a  comprehensive  representation  of 
the  kinetics  information  of  a  certain  fuel/air  system 
requires  the  three  independent  global  parameters 
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of  d>.  p  and  T,,|  adopted  herein  can  be  appreciated 
by  recognizing  that  the  rate  of  an  elementary  re¬ 
action  also  requires  at  least  three  fundamental  ther¬ 
modynamic  parameters.  These  parameters  are 
needed  to  represent  the  local  values  ol  the  concen¬ 
tration  of  at  least  one  ol  the  reacting  species,  the 
density  and  thereby  pressure  to  indicate  the  fre¬ 
quency  of  collision,  and  the  temperature  to  indicate 
the  energetics  of  collision  Therefore  m  the  pre¬ 
vious  representation  of 

sy  =  sy.  Kb.  p.  t„i  id)!1. 

the  concentration  and  temperature  sensitivities  ol 
the  elementary  reactions  base  not  been  indepen¬ 
dently  tested. 

The  above  discussion  basically  provides  the  mo¬ 
tivation  and  approach  of  the  present  investigation 
In  the  next  two  sections  the  experimental  meth¬ 
odology  and  results  will  be  presented  These  will 
be  followed  by  a  comparison  with  results  obtained 
from  detailed  numerical  simulation  by  assuming  dif¬ 
ferent  kinetic  schemes. 


Experimental  Methodology 

The  various  flow  variables  are  defined  in  Fig  la 
The  experiment  basically  involved  the  establish¬ 
ment  of  two  symmetrical,  planar,  nearlv-adiabatic 
flames  in  a  nozzle-generated  counterflow,  and  the 
determination  of  the  axial  velocity  profile  along  the 
centerline  of  the  flow  by  using  laser  Doppler  ve- 
locinietrv.  By  identifying  the  minimum  point  of  the 
velocity  profile  as  a  reference  upstream  flame  speed 
corresponding  to  the  imposed  strain  rate  k  which 
is  simply  the  velocity  gradient  ahead  of  it.  and  by 
plotting  SL  versus  k,  the  laminar  flame  speed  with¬ 
out  straining  S'l.  can  be  determined  through  linear 
extrapolation  to  K  =  0  (Fig.  lb)  in  accordance  with 
the  theoretical  prescription. 7  The  linear  extrapo¬ 
lation  is  meaningful  and  accurate  only  for  small  x al¬ 
lies  of  the  nondimensional  strain  rate 

DK 

ka  —  - -  <  <  1 . 

«sy.r 

where  ka  is  the  Karlovitz  number  and  I)  a  rep¬ 
resentative  mass  diffusivity. 

The  diameters  of  the  nozzles  used  m  the  present 
investigation  were  7.  10  and  14  mm.  They  had  high 
contraction  ratios  and  were  water-cooled  and  nitro¬ 
gen-shrouded.5  The  entire  burner  assembly  was 
housed  in  a  large  stainless  steel  chamber  with  ad¬ 
justable  pressure  and  continuous  ventilation  The 
cross-section  of  the  test  portion  of  the  chamber  was 
14.5  cm  by  14.5  cm.  Ignition  was  achieved  by  a 
continuous  high  energy  spark  produced  by  re  tract - 
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Ctg/Oi/Ni/Ar.  Ttd-22J0K,  P-1. Own.  0-0.76,  Ar/(AM-N2H>.8 


Cm/Oj/Nj/Ar.  Tad=2250K.  P=l.0alm.  0-0.76.  Ar/( Ar+N  i)=0.8 


(b) 

Fig.  1.  (a)  Typical  axial  velocity  profile  across  a 
stagnation  flame,  showing  the  definitions  of  K  and 
SL.  (b)  Dependence  of  the  flame  speed  SL  (K)  on 
the  strain  rate  K. 


able  electrodes.  The  LDV  operated  in  the  forward 
scattering  mode  with  0.3  pm  alumina  particle  seed¬ 
ing 

The  flame  temperature  was  independently  varied 
by  substituting  part  of  the  N2  in  air  by  equal  amount 
of  Ar  and  C02  for  higher  and  lower  temperatures 
respectively.  Figure  2  shows  the  different  amounts 
of  substitution  used  to  produce  the  mixtures  for  fixed 
adiabatic  flame  temperatures. 

It  is  important  to  mention  at  this  point  that  this 
inert  substitution  affects  not  only  the  flame  tem¬ 
perature  but  also  the  transport  properties  of  the 
mixture.  Specifically,  for  the  same  flame  tempera¬ 
ture  and  similar  equivalence  ratios,  Ar  and  C02 
substitutions  tend  to  increase  and  decrease  the  flame 
speed  respectively  through  changes  in  the  effective 
(X/cJ,  where  A  and  cp  are  the  mixture  thermal 
conductivity  and  specific  heat  respectively.  Since 
transport  properties  of  mixtures  are  known  with 
much  greater  accuracy  than  the  kinetics,  effects  due 
to  transport  properties  can  be  accurately  described 
such  that  kinetics  is  the  only  component  to  be 
studied. 

The  major  source  of  inaccuracy  in  the  present 
experiment  comes  from  the  difficulty  in  accurately 
controlling  the  individual  concentrations  of  the  var¬ 
ious  gaseous  components  in  the  mixture  prepara¬ 


tion  stage.  When  the  preparation  involves  only  three 
components,  as  in  the  case  of  CH4/Ar/02  in  Fig 
3,  the  control  is  quite  tight  and  the  data  have  higher 
accuracy.  However,  when  four  gaseous  components 
(i. e. ,  CHj,  02,  N2,  Ar  or  C02)  are  involved  as  in 
the  constant  flame  temperature  experiments,  the 
simultaneous  control  is  more  delicate.  Slight  un¬ 
certainties  in  the  mixture  concentration  could  lead 
to  corresponding  uncertainties  in  the  flame  tem¬ 
perature  and  equivalence  ratio.  It  is  estimated  that 
these  uncertainties  could  result  in  maximum  scat¬ 
ters  of  ±1.5  cm/s  in  the  measured  flame  speed  and 
±0.01  in  4>.  These  uncertainties  are  quite  small  and 
represent  an  improvement  over  previously  reported 
data  in  addition  to  the  elimination  of  systematic- 
strain  rate  effects. 


Experimental  Results 

Figure  3  is  the  traditional  plot  of  flame  speed 
S‘l  versus  equivalence  ratio  <f>  for  methane/Ar-air, 
methane/N2-air,  and  methane/C02-air  mixtures  at 
1  atm  pressure.  The  symbol  l-air  here,  as  in  Ar- 
air,  designates  an  air  mixture  in  which  all  the  Nj 
has  been  substituted  by  Ar.  All  of  the  present  re¬ 
sults  are  bounded  by  the  Ar-air  and  C02-air  data, 
from  the  top  and  bottom,  respectively.  We  further 
note  that  for  the  same  <f> ,  flames  above  and  below 
the  N2-air  data  have  higher  and  lower  temperatures 
respectively.  The  data  for  Ar-air  flames  obtained  by 


Flo,  2.  Relative  molar  amounts  of  the  various  in¬ 
erts  required  to  achieve  different  adiabatic  flame 
temperatures,  as  (unctions  of  4> 
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Fic.  3.  Experimental  laminar  flame  speeds  St  (<j>) 
for  methane/Ar-air,  methane/N2-air,  and  methane/ 
CO,-air  mixtures,  at  1  atm. 


0 

Fic.  4.  Experimental  laminar  flame  speeds  S’f  (<b) 
for  constant  flame  temperatures,  T^,  at  I  atm. 


Clingman  et  al.,lu  using  the  area  method,  show 
maximum  flame  speeds  nearly  40%  higher  than  those 
obtained  in  the  present  work. 

Figure  4  shows  the  measured  Sl  as  a  function  of 
<f>,  on  constant  flame  temperatures  Ta,i  and  at  l  atm 
The  experimental  results  of  Fig.  3  are  shown  as 
dashed  curves.  The  data  set  of  Taj  =  2,250  K  cor¬ 
responds  to  the  lowest  Ta,i  for  which  all  the  Sl  are 
higher  than  the  methane/N2-air  curve.  For  lower 
values  of  Ta<j,  C02  substitution  is  needed  to  obtain 
flame  speeds  below  those  of  the  methane/N2-air 
curve.  It  is  significant  to  note  that  Sl  varies  quite 
non-monotonically  with  4);  this  behavior  will  be  ex¬ 
plained  later. 

Figure  5  shows  the  experimental  S'l  (<t>)  for  var¬ 
ious  pressures  p  with  a  constant  Taj  =  2,250  K.  It 
is  seen  that,  for  the  same  <f>,  S£  decreases  with  in¬ 
creasing  pressure.  Since  the  eigenvalue  for  laminar 
flame  propagation  is  the  mass  burning  rate  (per  unit 
area)  p„Sl  instead  of  Sl,  where  pu  is  the  density 
of  the  unburnt  mixture,  the  data  of  Fig.  5  are  re¬ 
plotted  in  Fig.  6  as  puSl-  Now  we  see  that  puSt 
increases  with  increasing  pressure.  Since  for  a  given 
4>  one  can  define  an  overall  reaction  order  n(<J>)  such 
that  p„Sl  ~  pn/z,  the  increasing  trend  of  p„S£  in¬ 
dicates  positive  values  of  n.  Figure  7  shows  n  de¬ 
termined  as  a  function  of  p  for  a  given  <}>.  It  is  seen 
that  while  n  is  positive,  it  decreases  with  increasing 
p  for  a  given  <J>.  This  illustrates  the  progressive  im- 
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Fic.  5.  Experimental  laminar  flame  speeds  S|  (<t>. 
p),  for  T„i  =  2,250  K  and  pressures  between  0.25 
and  2  atm. 
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Fig.  6.  Experimental  mass  burning  rates  P„S'l  (<}>. 
p),  lor  T,,i  =  2,250  K  and  pressures  between  0.25 
and  2  atm 


portance  ol  three-bodv  termination  reactions  rela¬ 
tive  to  two-body  branching  reactions  us  pressure  in¬ 
creases  and  overall  reaction  rate  decreases.  It  is 
especially  ol  interest  to  note  that  the  importance  ol 
the  termination  reactions  is  manifested  even  for  the 
present  situations  involving  relatively  low  pressures 
and  strongly  burning  Haines. 

Numerical  Results 

Numerical  simulation  of  the  one-dimensional  flame 
propagation  was  conducted  by  using  the  flame  code 
developed  by  Kee  and  co-workers,11  14  with  the 
minor  modification  that  diffusion  of  species  k  is  given 
by  the  multicomponent  diffusion  relation1’ 

TO‘  =  i  (») «  -  vtl  a, 

without  the  assumption  of  trace  species  diffusion, 
where  X^.V.-  and  Dj*  are  respectively  the  molar 
fraction  and  diffusion  velocity  of  species  j  and  its 
mass  diffusivity  relative  to  species  k.  Overall  mass 
conservation  is  achieved  by  applying  Eq.  (1)  to  the 
(N  -  1)  species  while  using 


Fic.  7.  Experimentally  determined  overall  reac¬ 
tion  orders  for  lean,  stoichiometric  and  rich  CH ,/ 
Oj/N2/Ar  flames  for  =  2,250  K. 


N 

2  Y>  v,  =  0  (2) 

i 

for  the  Nth  species,  where  Y,  is  the  mass  fraction. 

In  adopting  the  kinetic  mechanisms  of  methane 
oxidation  needed  for  simulation,  we  note  that  in  a 
previous  study  on  methane/air  flames, lh  an  exten¬ 
sive  evaluation  was  conducted  for  the  various  ki¬ 
netic  mechanisms  reported  in  the  literature.  Through 
such  an  evaluation,  two  mechanisms  emerge  as  being 
highly  predictive  of  the  experimentally-determined 
laminar  flame  speeds.  These  are  the  Ci  mechanism 
of  Kee  et  al.13  and  the  full  C2  mechanism  of  War- 
natz,2  which  are  predictive  up  to  4>  =  1.2  to  1.3. 
These  comparisons  are  reproduced  for  the  metli- 
ane/N2-air  data  in  Fig.  8.  The  C|  mechanism  has 
the  advantages  of  being  predictive  up  to  moder¬ 
ately  rich  conditions  and  is  also  considerably  less 
time  consuming  for  simulation.  Furthermore,  it  has 
also  been  suggested  '  that  the  methane/air  laminar 
flame  speeds  may  not  be  sensitive  to  C2-reactions 
up  to  <f>  =  1.3.  On  the  other  hand  the  more  com¬ 
plex  C2  mechanism  has  the  advantage  of  extended 
applicability  over  the  entire  concentration  range  It 
may  also  be  noted  that  since  development  of  these 
two  mechanisms  have  probably  involved  some  im¬ 
plicit  calibration  around  the  stoichiometric  condi¬ 
tion,  the  C|  part  of  the  full  C2  mechanism  ol 
Warnatz2  slightly  underpredicts  the  experimental 
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(lata  and  therefore  does  not  degenerate  to  the  C| 
nit  rhanism  ol  kec  et  al.l! 

In  the  present  simulation  these  two  mechanisms 
are  subjected  to  further  validation  by  comparing  the 
calculated  results  with  the  present  data  which  in¬ 
corporate  variations  in  concentration,  pressure  and 
tlame  temperature.  The  constraints  imposed  on  (In¬ 
validation  criterion  are  therefore  more  stringent  than 
those  ol  the  previous  validation  studies. 

Figure  >S  shows  the  calculated  tlame  speeds  of 
methane/Ar-air  mixtures  by  using  both  the  Ct  and 
Cs  mechanisms.  The  flame  speeds  ol  methane /(X)2- 
air  have  not  been  calculated  because  they  have  very 
low  values  and  the  calculations  are  very  difficult  to 
converge.  These  calculated  results  are  compared  with 
the  experimental  data,  shown  in  Fig.  8  as  dashed 
curves  which  are  reproduced  from  the  solid  curves 
of  Fig.  3.  The  comparison  shows  that  the  (fo  mech¬ 
anism  predicts  well  the  flame  speeds  of  methane/ 
\ 2-air  and  methane/Ar-air  mixtures  even  up  to 
about  (Jr  =  1  2  to  1.3.  The  C2  mechanism  slightly 
underpredict::  the  flame  speed  lor  Argon  mixtures 
up  to  d>  ~  1.0.  although  the  diflerenee  is  too  small 
to  be  significant.  For  richer  flames,  however,  it 
overpredicts  the  flame  speed  by  somewhat  larger 
(mounts. 

Figure  9  compares  the  constant  T,„i  data  of  Fig. 
-1,  shown  as  the  dashed  curves  here.  Most  ol  the 
calculations  were  conducted  with  the  C  i ,  mecha¬ 
nism  because  numerical  solution  with  the  C>  mccli- 
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Fic.  H.  Comparison  between  experimentally  and 
numerically  determined  S’,’  (<t>)  lor  methane/Ar-air 
mixtures  at  1  atm.  using  C,  and  C,  kinetic  schemes. 
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I'll.  9  Comparison  between  cxperimcntalh  and 
numerically  determined  Si'  ld>)  lor  various  1 ,,  at  I 
atm.  using  C,  and  C\  kinetic  schemes 

anism  is  very  time  consuming.  It  is  seen  that  the 
calculated  results  satisfactorily  reproduce  the  highly 
non-monotonic  behavior  of  the  experimental  data 
The  C|  mechanism  again  shows  good  agreement 
while  the  limited  amount  ol  solutions  with  (.'< 
mechanism  yield  higher  flame  speeds  lor  the  richer 
mixtures.  It  is  perhaps  somewhat  surprising  that  the 
agreement  with  the  C|  mechanism  is  better  lor  the 
methane/Ns-air  and  methane/Ar-air  flames  ol  Fig. 
8  than  those  ol  the  present  constant  T„|  flames  which 
are  Ixnmded  by  them.  This  is  believed  to  be  caused 
by  the  increased  scatter  in  the  data  (or  the  con- 
stant-Tal|  experiments,  as  mentioned  earlier. 

Finally,  Fig.  10  shows  the  comparison  for  the 
constant  pressure  data,  at  T„|  =  2.250  k  The 
agreement  can  again  be  considered  to  be  mostly 
satisfactory. 


Discussions 

In  order  to  explain  the  present  experimental  and 
numerical  results,  which  largely  agree  with  each 
other  and  are  therefore  consistent,  it  is  necessary 
to  consider  both  transport  (-fleets  due  to  inert  sub¬ 
stitution  and  kinetic  (-fleets  due  to  changes  in  stoi¬ 
chiometry .  Let  us  lirst  discuss  the  non-monotonic 
behavior  ol,  say  ,  the  I  atm  and  T,„i  =  2.25(1  k  data 
ol  Fig  (i.  p„S'|  is  adopted  lor  discussion  because  it 
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Fic.  10.  Comparison  between  experimentally  and 
numerically  determined  S£  (<}>,  p)  for  T„  =  2,250 
K,  using  the  C,  kinetic  scheme. 


is  the  actual  eigenvalue  of  flame  propagation.  The 
results  show  that  p„S'l  increases  in  both  the  lean 
and  rich  directions  from  the  stoichiometric  regime. 
Since  flame  temperature  is  held  constant  on  a  given 
curve,  and  since  p„Sl  should  decrease  based  on 
concentration  effect  alone,  the  observed  increase  is 
therefore  due  to  the  higher  values  of  (\/c,J  of  the 
mixture  with  increasing  concentration  of  Ar.  The 
same  explanation  can  be  extended  to  C02-diluted 
flames. 

The  behavior  of  p„Sl  with  further  departures  from 
stoichiometry  becomes  quite  asymmetrical  in  that  it 
continues  to  increase  for  leaner  mixtures  but  de¬ 
creases  for  richer  mixtures.  The  increase  on  the  lean 
side  is  of  course  partly  still  due  to  transport  effects. 
We  have,  however,  found  that  kinetic  effects  could 
also  be  partly  responsible  in  that  among  the  three 
leading-order  reactions 

H  +  02  =  OH  +  0  (A) 

H  +  02  +  M  =  HOz  +  M  (B) 

CH,  +  H  +  M=  CH4  +  M,  (Cl 

the  reaction  rates  of  (A)  and  (B)  arc  minimally  af¬ 
fected  with  decreasing  <J>.  while  the  recombination 
rate  of  reaction  (C)  is  substantially  reduced,  by  about 


40%,  from  <f>  =  1.0  to  0.7.  Thus  the  overall  reac¬ 
tion  rate  is  enhanced 

The  decreasing  trend  for  rich  flames  is  kinetic  in 
nature.  For  moderately  rich  flames  (1  <  d>  S  1.3) 
the  forward  reaction  rates  of  (A)  and  (B)  respec¬ 
tively  decrease  and  increase  with  increasing  <f>.  This 
eventually  causes  a  significant  decrease  in  puS'f  in 
spite  of  the  positive  influence  of  (\/cp).  For  very 
rich  flames  (<t>  >  1.3)  the  forward  reaction  rate  of 
(A)  increases  again,  although  the  forward  reaction 
of 

C2H5  +  H  =  CH3  +  CH3  (D) 

also  becomes  important  and  tends  to  slow  down  (A) 
by  competing  for  the  H  radicals. 

Figure  6  also  shows  that  the  curves  “flatten  out’’ 
with  decreasing  pressure.  This  can  be  simply  ex¬ 
plained  by  considering  the  relations 

,,  dlpuSV ) 

p„s?.  ~  p"'2/«i>).  -  ,r-/W.  1.31 

<7<J) 

Thus  for  a  given  4>  and  for  the  range  of  n  encoun¬ 
tered  in  the  present  study  (Fig.  7),  the  slopes  of 
the  curves  in  Fig.  6  decrease  with  decreasing  p. 

It  is  also  of  interest  to  assess  whether  C02  can 
have  any  significant  direct  kinetic  effect  resulting 
from  its  substitution.  A  comparison  between  the 
CHj/NVair  flame  and  a  CH4/02/N’2/C02  flame 
with  Taj  -  2,150  K,  while  keeping  p  =  1  atm  and 
4>  =  l  for  both  flames,  shows  that  the  sensitivity 
on  S'l  of  all  the  crucial  reactions  is  minimally  af¬ 
fected  and  the  maximum  CO  concentrations  also 
remain  almost  the  same.  Thus  C02  substitution  does 
not  directly  contribute  to  kinetic  effects. 

In  view  of  the  higher  flame  speeds  predicted  by 
the  C2  mechanism  for  rich  flames,  a  sensitivity- 
analysis  was  conducted  for  the  T,l(j  =  2,250  K  and 
p  =  1  atm  flames  at  tf)  =  1.15  and  1.435.  The  re¬ 
sults  show  that  reactions  of  leading  order  impor¬ 
tance  are  (A)  and  (C)  as  would  be  expected  For 
the  reactions  at  the  next  order  of  importance,  it  is 
found  that  the  sensitivity  of  reaction  (D)  increases 
by  a  factor  of  3.2  from  <t>  =  1.15  to  d>  =  1.435. 
while  the  sensitivities  of  all  the  other  reactions  arc 
not  affected  much.  By  slightly  modifying  the  kinetic 
constants  of  this  reaction,  the  calculated  flame  speeds 
can  be  substantially  lowered  and  improved  agree¬ 
ment  with  the  experimental  data  is  achieved.  How¬ 
ever  we  have  not  attempted  to  calibrate  these  ki¬ 
netic  constants  against  our  data.  Rather  we  feel  it 
is  more  appropriate  for  a  study  of  the  present  na¬ 
ture  to  recommend  that  the  rate  constants  of  this 
elementary  reaction  may  need  to  be  either  re-ex¬ 
amined  and/or  re-determined  through  basic  kinet¬ 
ics  studies. 


LAMINAR  FLAMES:  PROPAGATION  AND  EXTINCTION 


Concluding  Remarks 

There  are  two  aspects  of  contributions  from  the 
present  study.  First,  we  have  performed  the  tedi¬ 
ous  task  of  accurately  determining  the  stretch-free, 
laminar  flame  speeds  of  methane/(Ar  N2,  C02)-air 
mixtures  over  extensive  parametric  ranges  of  stoi¬ 
chiometry,  pressure,  and  flame  temperature.  These 
data  are  expected  to  be  useful  for  the  partial  vali¬ 
dation  of  proposed  kinetic  mechanisms.  Second,  we 
have  also  performed  detailed  numerical  calculations 
which  partially  validated  a  Ci  mechanism  and  a  full 
C2  mechanism  for  lean  to  moderately  rich  mixtures. 
The  need  for  possible  modification  of  the  C2  mech¬ 
anism  has  also  been  identified.  It  is  clear  that  such 
concrete  conclusions  cannot  be  made  without  an 
accurate  experimental  data  set  as  a  reference 

The  present  study  also  emphasizes  that  agree¬ 
ment  in  the  flame  speed  is  only  a  necessary  but  not 
sufficient  criterion  for  the  validation  of  a  kinetic 
mechanism  For  example,  we  have  shown  that  cal¬ 
culated  flame  speeds  using  the  Cj  and  C2  mecha¬ 
nisms  both  agree  well  with  the  experimental  data 
even  for  moderately  rich  mixtures,  while  there  ex¬ 
ist  ample  kineticallv-related  experimental  phenom¬ 
ena  which  can  only  be  explained  by  including  the 
C2  reactions.  Thus  there  exists  an  urgent  need  for 
accurate  experimentally-determined  temperature  and 
concentration  profiles  of  laminar  flames  which  can 
be  used  for  full  validation.  In  conducting  such  ex¬ 
periments  or  validations,  we  again  caution  the  need 
to  either  remove  or  account  for  the  effects  due  to 
flame  straining.  These  effects  can  be  especially  sig¬ 
nificant  in  the  potential  modification  of  the  local 
concentrations  of  such  highly  mobile  species  as  the 
important  hydrogen  radicals. 
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COMMENTS 


J.  Warnatz,  Physikalisch-Chemisches  Institut,  Fed. 
Rep.  of  Germany.  I  do  not  agree  with  your  con¬ 
clusion  that  the  high  flame  velocities  calculated  on 
the  rich  side  with  the  C2-mechanism  are  due  to  de¬ 
ficiencies  in  the  C2  chemistry.  To  our  present 
knowledge,  the  weakest  point  in  these  mechanisms 
is  CHj  oxidation  though  we  do  not  know  the  re- 
action(s)  responsible.  Removal  of  this  deficiency 
(hopefully  in  the  near  future)  should  lower  both  the 


results  with  C,  and  C2  mechanism  and  low  results 
for  the  C,  mechanism. 

Authors  Reply.  We  appreciate  very  much  your 
comment.  We  have  adopted  a  neutral  tone  in  re¬ 
plying  because  there  just  do  not  exist  enough  in¬ 
formation  for  a  definitive  statement.  We  hope  this 
reply  is  satisfactory  to  you.  Let  us  know  when  you 
get  additional  insight/results  regarding  this  issue. 
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Recent  advances  in  the  understanding  of  the  structure,  propagation,  and  extinction  of  lam¬ 
inar  (lames  under  the  influence  of  stretch,  as  manifested  by  the  existence  of  (lame  curvature, 
(low  nonuniformity,  and  (lame  motion,  are  reviewed.  The  emphasis  is  on  premixed  flames 
because  of  the  richness  and  subtlety  of  the  phenomena  involved  The  review  distinguishes 
the  influences  of  the  tangential  and  normal  velocity  gradients  on  the  flame  response,  both 
at  the  hydrodynamic  scale  and  within  the  flame  structure,  and  emphasizes  the  importance 
of  the  preferential  diffusion  nature  of  heat  and  mass  transport,  as  well  as  the  extent  to  which 
the  flame  can  freely  adjust  its  location  in  response  to  stretch  in  order  to  achieve  complete 
reaction.  It  is  then  demonstrated  that  stretch  has  only  minimal  effect  on  an  adiabatic,  unre¬ 
strained,  diffusionally-neutral  flame  with  complete  reaction  in  that  the  temperature,  propa¬ 
gation  rate,  and  thickness  of  the  flame  are  invariant  to  stretch,  and  that  stretch  alone  cannot 
extinguish  such  a  flame.  In  the  presence  of  preferential  diffusion  and/or  when  the  flame 
movement  is  restrained,  the  response  of  the  flame  to  stretch  becomes  more  sensitive  and 
extinction  is  also  possible.  The  concept  of  flame  stretch  is  applied  to  interpret  such  practical 
flame  phenomena  as  flame  stabilization  and  flame-front  instability,  determination  of  laminar 
flame  speeds  and  flammability  limits,  concentration  and  temperature  modifications  in  flame 
chemistry,  and  modeling  of  turbulent  flames.  The  properties  of  stretched  diffusion  flames  a  re 
then  briefly  discussed.  The  review  closes  with  suggestions  for  further  research 


Introduction 

Historically,  the  concept  of  flame  stretch  was  first 
introduced  by  Karlovitz1  to  describe  flame  extinc¬ 
tion  in  velocity  gradients.  It  was  subsequently 
adopted  by  Lewis  and  von  Elbe*  to  explain  and 
quantify  the  various  phenomena  associated  with 
flame  stabilization.  The  flame  curvature  aspects  of 
stretch  also  formed  the  basis  for  the  study  of  flame- 
front  instability  by  Markstein.3  In  the  past  fifteen 
years  or  so  significant  advances  have  been  made  in 
our  understanding  of  flame  stretch  and  its  influence 
on  the  flame  structure.  These  include  a  fundamen¬ 
tal  generalization  of  the  definition  of  stretch,  math¬ 
ematically  rigorous  analyses  of  the  structure  and 
propagation  of  stretched  flames,  experimental  and 
numerical  verification  and  quantification  of  pre¬ 
dicted  properties  via  model  flames,  and  re-inter- 
pretation  of  certain  flame  phenomena.  Some  recent 
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reviews4  1 1  have  summarized  the  various  aspects  of 
these  advances. 

The  objectives  of  the  present  review  are  to  syn¬ 
thesize  the  current  understanding  of  stretched  flames 
from  a  unified  viewpoint,  and  to  identify  areas  for 
further  research.  In  the  next  section  the  mathe¬ 
matical  definition  of  stretch  and  its  physical  inter¬ 
pretation  will  be  presented.  We  then  study  the  in¬ 
fluence  of  stretch  on  the  flame  structure  through 
its  coupled  effect  with  heat  and  mass  diffusion.  Since 
the  flame  response  turns  out  to  depend  critically  on 
the  unequal  rates  of  heat  and  species  diffusion,  as 
signified  by  Le,j  t4  1,  the  flame  behavior  without, 
and  with,  preferential  diffusion  arc  separately  dis¬ 
cussed.  We  then  review  the  role  of  stretch  in  such 
practical  flame  problems  as  flame  stabilization,  flame 
speed  determination,  flammability  limits,  flame-front 
instability,  concentration  and  temperature  modifi¬ 
cations  as  related  to  studies  of  chemical  kinetics, 
and  modeling  of  turbulent  flames. 

While  much  of  the  present  review  is  on  pre¬ 
mixed  flames  because  of  the  richness  and  subtlety 
of  the  response  of  a  premixed  flame  to  stretch,  the 
properties  of  stretched  diflusion  flames  are  also 
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briefly  discussed  and  contrasted  with  those  of  pre¬ 
mised  flames.  The  paper  closes  with  suggestions  for 
further  research. 


Definition  of  Stretch  and  Preliminary 
Considerations 

Understanding  of  the  structure  and  propagation 
of  premised  flames,  both  at  the  fundamental  level 
as  well  as  for  the  analysis  of  practical  combustion 
phenomena,  is  frequently  based  on  the  classical 
model  involving  the  steady  propagation  of  an  adi¬ 
abatic,  one-dimensional,  planar  flame  into  a  com¬ 
bustible  medium  of  temperature  T„  and  reactant 
concentrations  Y,  as  shown  in  Fig.  la  in  the  flame¬ 
stationary  frame.  This  flame  is  characterized  by  its 
final,  equilibrium  temperature  T£,  which  is  simply 
the  adiabatic  flame  temperature  T„j  of  the  mixture, 
and  by  its  mass  buining  rate  (per  unit  area)  in", 
which  is  a  function  of  the  exothermicity,  diffusivitv, 
and  reactivity  of  the  mixture;  the  superscript  “o’ 
designates  properties  pertaining  to  this  specific  flame 
propagation  mode.  Thus  if  we  symbolize  the  de¬ 
pendence  on  exothermicity  by  the  specific  heat  re¬ 
lease  q,  diffusivitv  by  the  Lewis  number  Ley  =  X/ 
(c^pD.j),  and  reactivity  by  the  reaction  rate  tc*  of' 
reaction  k,  then 


-3 

'g 

II 

(Lit 

tj:  =  n«i). 

(1.2) 

Furthermore,  if  the  characteristic  chemical  times 
ol  the  major  heat  release  reactions  are  much  shorter 
than  those  of  heat  and  mass  diffusion,  then  the  flame 


structure  can  be  considered  to  consist  of  a  thin  re¬ 
action  zone  preceded  by  a  broader  transport-dom¬ 
inated,  diffusive-convective  zone,  which  are  respec¬ 
tively  characterized  by  thicknesses  5|<  and  8,, 
satisfying  8,,  <<  8|>.  From  in  =  pu  we  can  then 
define  an  unburnt,  upstream  flame  speed  S"  -  in"/ 
p„  at  the  upstream  boundary  of  the  transport  zone 
and  a  burnt,  downstream  flame  speed  SJJ 
in"/ p)’,  for  the  entire  reaction  zone  because  velocits 
does  not  change  much  owing  to  the  thinness  of  tins 
zone 

The  flame  configuration  of  Fig.  la  is  only  an 
idealization,  in  realistic  situations  the  upstream  floss 
can  be  nonuniform  while  the  flame  can  also  be 
curved  and  nonstationary  (Fig.  lb),  producing  what 
are  collectively  known  as  stretch  effects  Further¬ 
more.  the  flame  can  also  be  globally  nonadiabatic 
through  exchanges  in  heat  and  radicals  svith  the 
surrounding.  Thus  if  sve  symbolize  influences  ol 
stretch  by  a  Karlovitz  number  Ka,  to  be  defined 
later,  and  system  nonadiahaticity  by  a  "loss  param¬ 
eter  L,  then  the  flame  behavior  is  expected  to  be 
described  by 

m  =  m(q.  Ley,  u*;  Ka,  L)  i  1 . 3 ‘ 

Tj.  =  T,,(</.  fa-,,.  Ka.  L1  I  f 

which  should  deviate  from  in"  and  17,  of  tin'  ideal¬ 
ized  situation.  It  may  be  noted  that  although  m  <1.  3 
and  ( 1  .-If  sve  have  purposely  separated  out  the 
functional  dependence  into  the  two  groups  of  pa¬ 
rameters  (q, Lev,tc\)  and  (Ka.L),  strong  coupling  be¬ 
tween  the  individual  effects  is  expected.  For  ex¬ 
ample.  while  T/’,  depends  only  on  </.  the  presence 
of  stretch  and/or  loss  mechanism  can  cause  T/,  to 


Fit.  1  Schematics  of  the  structure  of  (a)  a  one -duucnsioiiul  planar  flame, 
a  nonunilorm  flow. 
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V 


Fig.  2.  Schematic  of  surface  element  s  with  ve¬ 
locity  V  and  unit  normal  vector  n  in  a  flow  field  v. 


also  depend  on  the  diffusive  and  reactive  aspects  of 
the  problem.  We  also  note  that  the  presence  of  flow 
nonuniformity  and  flame  curvature  introduces  an 
additional,  hydrodynamic,  length  scale,  8H,  into  the 
problem. 

To  define  stretch,  we  first  identify  an  arbitrary  sur¬ 
face  s,  say  an  isotherm,  within  the  flame  structure 
(Fig.  2).  This  surface  has  a  velocity  V  while  the 
fluid  has  a  velocity  v.  A  general  definition8  10  12-14 
of  stretch  at  any  point  on  this  surface  is  the  time 
derivative  of  the  logarithm  of  the  area  of  an  infin¬ 
itesimal  element  of  the  surface, 


with  the  boundary  of  this  surface  element  moving 
tangentially  to  the  surface  at  the  local  tangential 
component  of  the  fluid  velocity.  In  terms  of  the  flow 


variables,  it  can  be  shown1313,16  that 

K  =  |V,-v,  +  (V-nXV-n)},.  (2.2) 

where  V,  and  v,  are  the  tangential  components  of 

V  and  v  evaluated  at  the  surface,  and  n  is  the  unit 
normal  vector  of  the  surface,  pointed  in  the  direc¬ 
tion  of  the  unburnt  gas.  The  stretch  rate  K  has  the 
unit  of  sec’ '. 

Equation  (2.2)  shows  the  two  sources  of  stretch 
a  flame  can  be  subjected  to.  The  first  term  repre¬ 
sents  the  influence  of  flow  nonuniformitv  along  the 
flame  surface  Since16 

{V,-v,},  =  {-n*V  X  (v  x  n)}„  (2.3) 

this  term  embodies  the  effects  due  to  flow  non¬ 
uniformity  through  v  and  flame  curvature  through 

V  ■  n.  Furthermore,  this  term  exists  only  if  the  flow 
is  oblique  to  the  surface  (v  x  n  ¥>  0).  The  second- 
term  in  Eq.  (2.2)  represents  stretch  experienced  by 
a  nonstationary  flame  through  V,  although  the  flame 
also  has  to  be  curved  because  (V  •  n)  vanishes  oth¬ 
erwise.  These  three  stretch-induced  effects  can  be 
separately  referred  to  as  those  due  to  aerodynamic 
straining,  flame  curvature,  and  flame  motion.  We 
further  note  that  since  (heat  and  mass)  diffusion  is 
parallel  to  n,  the  non-orthogonality  requirement  of 
(v  x  n)  leads  us  to  anticipate  the  importance  of  dif¬ 
fusive  transport  in  the  dynamics  of  stretched  flames, 
even  though  the  discussion  so  far  has  been  kine¬ 
matic  in  nature 

As  examples,  let  us  compute  the  stretch  rate  K 
for  the  three  common  flame  configurations15  shown 
in  Fig.  3.  The  flames  are  assumed  to  be  infinitely 
thin  for  simplicity,  thus  the  stretched  surface  is  the 
flame.  Figure  3a  shows  a  planar  flame  situated  in 
a  divergent  stagnation  flow.  Assuming  potential  flow, 
the  velocity  vector  is  v  =  {[a/(a  +  l)]x,  -tfi/,  0}, 


u„ 


(a)  Stagnation  Flame  (b)  Propagating  Spherical  (c)  Axisymmetric  Bunsen 

Flame  Flame 

Fig.  3.  Examples  of  stretched  (lames:  (a)  Stagnation  flame,  (!>)  outwardly  propagating  spherical  flame, 
and  (c)  axisymmetric  Bunsen  flame. 
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where  a  is  the  strain  rate  of  the  flow,  a  =  0, 1  for 
cartesian  and  cylindrical  coordinates  respectively, 
and  the  x-  and  {/-velocities  in  the  cylindrical  coor¬ 
dinate  are  those  in  the  radial  and  axial  directions 
respectively.  Since  V  =  0,  Eq.  (2.2)  readily  yields 

K  =  u.  (2.4) 


Similarly  it  can  be  shown  that  for  an  outwardly 
propagating  spherical  flame  (Fig.  3b), 


2  r% 

%  dt  ' 


(2.5) 


where  Ry  is  the  instantaneous  flame  radius,  while 
for  the  conical  surface  of  an  axisymmetric  Bunsen 
flame  (Fig.  3c), 


K  =  - 


sin  9 
2R, 


(2.6) 


where  u*  is  the  uniform  upstream  velocity,  0  the 
apex  angle,  and  Ry  the  flame  radius  at  a  given  axial 
cross  section.  Note  that  K  >  0  for  the  positively- 
stretched  forward  stagnation  flame  and  outwardly- 
propagating  spherical  flame,  while  K  <  0  for  the 
negatively-stretched,  or  compressed,  Bunsen  flame. 
The  signs  of  Figs.  3a,  b  and  c  are  also  reversed  for 
the  corresponding  flame  configurations  of  the  rear¬ 
ward  stagnation  flame,  the  inwardly-propagating 
flame,  and  the  inverted  flame  if  only  the  effect  of 
flame  structure  is  considered.  We  shall  demon¬ 
strate  subsequently  that  the  flame  response  can  be 
qualitatively  reversed  for  positive  and  negative  val¬ 
ues  of  K 

Recently  the  flame  behavior  in  rotating  flow  fields 
has  also  been  studied  by  using  a  rotating  stagnation 
flow  system17'20  and  via  tubular  flames.21  25  Stretch 
rates  experienced  by  these  flames  can  be  individ¬ 
ually  calculated. 

There  are  also  stretchless  flames.  Examples  are 
the  steady  and  unsteady  one-dimensional  planar 
flame,  and  the  stationary  cylindrical  and  spherical 
flames  sustained  by  line  and  point  sources  respec¬ 
tively. 

The  stretch  rate  K  can  be  nondimensionalized  by 
the  characteristic  flow  time  within  the  flame. 


Ka  =  (—|K  = 
\S" 


D„ 

,(S")2J 


K, 


(2.7) 


which  is  the  Karlovitz  number  mentioned  earlier. 
The  significance  of  the  functional  form  of  Ka  will 
be  discussed  in  relation  to  Eq.  (4.4).  We  further 
note  that  since  stretch  is  associated  with  a  given 
surface,  for  a  flame  of  finite  thickness  we  can  define 
an  upstream  and  a  downstream  stretch  rate,  K„  and 
K i,.  for  the  unburnt  and  burnt  states  of  the  mixture 


respectively.  From  the  viewpoint  of  asymptotic 
analysis  with  one-step  reactions,  K;,  is  a  more  nat¬ 
ural  representation  of  the  stretch  rate  because  the 
reaction  zone  is  thin  and  is  therefore  indicative  of 
a  well-defined  'surface.'  In  realistic  situations, 
however,  the  multi-step  nature  of  the  chemical  re¬ 
action  schemes,  and  the  possibility  that  complete 
reaction  may  not  be  achievable,  render  this  reac¬ 
tion  surface  less  sharply  defined.  On  the  other  hand, 
both  numerical  and  experimental  results  have  shown 
that  K„,  which  is  not  directly  affected  by  diffusion 
and  reaction,  can  be  easily  identified.  Furthermore. 
K„  can  also  be  considered  as  an  imposed  parameter 
while  K/,  is  strongly  influenced  by  the  flame  re¬ 
sponse.  Thus  we  shall  adopt  Ku  in  subsequent  dis¬ 
cussions. 

The  influence  of  stretch  on  the  flame  response 
can  be  discussed  based  on  the  scales  as  w’ell  as  the 
tangential  and  normal  components  of  the  flow  field 
at  the  flame.  Let  us  first  consider  stretch  at  the 
hydrodynamic  scale  (Fig.  lb).  Here  diffusion  is  not 
important  and  the  entire  flame  can  be  considered 
to  collapse  into  a  flame  sheet.  Thus  the  burning 
intensity  of  the  flame,  such  as  its  temperature  T;, 
and  burning  rate  m„,  are  not  affected  by  stretch. 
The  presence  of  stretch,  through  the  tangential  ve¬ 
locity  gradient  at  the  flame  and  quantified  by  the 
Karlovitz  number  Ka,  changes  the  flame  surface  area 
A  and  consequently  the  volumetric  burning  rate 
m„A.  Thus  positive  stretch  increases  the  volumetric 
burning  rate  while  negative  stretch  decreases  it  We 
shall  refer  to  this  stretch  as  hydrodynamic  stretch 
The  role  of  the  normal  velocity  gradient  is  to  allow 
adjustment  of  the  flame  location  in  the  normal  di¬ 
rection  so  that  the  flame  situates  where  the  local 
flame  speed  balances  the  local  normal  velocity,  that 
is  S„  =  —  v„.  Thus  the  combined  effects  of  the  tan¬ 
gential  and  normal  velocity  gradients  are  the  dis¬ 
placement  of  the  flame  surface,  distortion  of  its  to¬ 
pography,  and  modification  of  the  volumetric  burning 
rate. 

In  the  transport  zone  the  tangential  velocity  gra¬ 
dient  directly  affect  the  normal  mass  flux  nq,  en¬ 
tering  the  reaction  zone.  Furthermore,  through  in¬ 
teraction  with  heat  and  mass  diffusion,  it  can  also 
modify  the  temperature  and  concentration  profiles 
in  this  zone  and  consequently  the  burning  inten¬ 
sity.  T/,  and  m/„  in  the  reaction  zone.  We  shall  re¬ 
fer  to  this  stretch  as  flame  stretch  The  normal  ve¬ 
locity  gradient  in  this  zone  affects  the  residence  time 
within  the  reaction  zone  and  consequently  T/,  and 
the  completeness  of  reaction.  A  measure  of  this  in¬ 
fluence  is  the  Damkohler  number  Da,  defined  as 
the  ratio  of  the  characteristic  flow  time  to  the  char¬ 
acteristic  reaction  time. 

Finally,  because  of  the  secondary  importance  of 
convective  transport  in  the  reaction  zone,  stretch 
has  little  influence  here  except  tor  highly  stretched 
situations 
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It  is  important  to  emphasize  the  different  roles 
played  by  the  tangential  and  normal  velocity  com¬ 
ponents  in  affecting  the  flame  response.  Since  these 
two  components  are  coupled  through  continuity, 
ditiusiunal  interaction,  and  thermal  expansion,  there 
has  been  some  contusion  in  the  literature  in  distin¬ 
guishing  effects  caused  by  stretch  Irom  those  caused 
by  residence  time,  or  Damkohler  number,  consid¬ 
erations. 

We  also  note  that  hydrodynamic  stretch  and  (lame 
stretch  can  he  strongly  coupled  in  that  the  hydro- 
dynamic  stretch  imposes  the  stretch  intensity  at  the 
(lame  needed  for  flame  stretch  calculations,  while 
the  flame  stretch  not  only  yields  the  burning  rate 
at  the  flame  surface  needed  for  the  calculation  of 
the  hydrodynamic  flow  field,  but  it  also  predicts  such 
critical  phenomena  as  flame  extinction. 

The  final  point  to  emphasize  is  the  flexibility  with 
which  a  premixed  flame  can  adjust  its  location  to 
accommodate  changes  in  the  normal  velocity  gra¬ 
dient  Thus  a  change  in  the  stretch  rate  does  not 
necessarily  lead  to  a  change  of  corresponding  extent 
in  the  residence  time.  We  shall  call  a  flame  with 
total  freedom  of  adjustment  an  unrestrained  flame. 

The  influence  of  hydrodynamic  stretch  is  quite 
straightforward  and  will  not  be  further  discussed. 
The  influence  of  flame  stretch  will  be  studied  in 
the  next  two  sections. 


Flame  Stretch  in  Diffusionally-. Neutral  Mixtures 

Adiabatic  l  unstrained  Flames: 

This  is  perhaps  the  simplest  situation  involving 
stretched  flames,  with  L  =  0,  Lc^  =  1.  The  flame 
is  also  assumed  to  be  able  to  freely  adjust  its  po¬ 


sition  in  response  to  changes  in  stretch  such  that 
complete  reaction  can  be  achieved  Simple  as  it  mas 
he  considerable  disagreement  exists  ill  the  uitei- 
prctation  of  its  basic  phenomena  For  example 
Karlov  it'/1  and  Lewis  and  son  Kibe"  argued  that 
with  increasing  flow  divergence  the  gradients  nl 
temperature  and  concentration  become  Mattel  Midi 
that  extinction  cxcntualh  occurs  as  a  result  of  the 
continuous  reduction  in  the  temperature  at  the  re¬ 
action  zone.  On  the  other  hand.  Klimov-1’  and  then 
Williams*  suggested  that  (positive)  stretch  tends  to 
thin  the  flame  sheet  convectixeb  and  to  steepen 
the  temperature  and  concentration  gradients  within 
it.  Conflicting  statements  abound  in  the  literature 
To  demonstrate  the  influence  of  stretch,  let  us 
first  consider  the  response  of  a  planar  adiabatic  flame 
to  aerodynamic  stretching  (Fig  far  A  simple  phe¬ 
nomenological  analvsis-'  -S  for  the  control  volume 
bounded  In  the  streamlines  and  the  boundaries  of 
the  fl  ame  vields  il  lations  tor  the  overall  mass  eon- 
solvation,  overall  energv  conservation.  mass  hal 
ance  at  the  reaction  zone  and  energv  balance  at  tin 
reaction  zone,  which  are  rexpcctivciv  given  In 

m„A„  -  mi.  A/,  i  l  l 

(T,.  -T„t 

I ii  —  I  „)  *  A -  tA i,  -  A,,’ 
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-  ^  nk/,ftu</  A/.  3  2 
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\  -  -  \nki,hnq  i-h  L 

bn 
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Km;  4  Phenotncmology  of  an  aerodynamic-ally -stretched  planar  flame  iai  (Correct  representation  allow  mu 
for  the  orthogonality  of  the  diffusive  flux,  (hi  meorreet  representation  suppressing  the  effeet  of  stn-teh 
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where  e,,,  \  and  (pD)  are  assumed  to  be  constants. 
ki,  =  B  exp(  —  T„/T/J  is  the  reaction  rate  constant 
evaluated  at  the  reaction  /.one,  and  Y|t  is  a  char¬ 
acteristic  reactant  concentration  within  the  reaction 
/one. 

Approximating  YR  =  eY„  and  8R  =  e8D,  where 
e  =  T^/T,,,  Eqs.  (3.1)  to  (3.4)  constitute  four  al¬ 
gebraic  relations  from  which  we  can  solve  for  the 
four  tlame  properties,  namely  the  flame  tempera¬ 
ture  Tj„  the  upstream  burning  rate  mu,  the  down¬ 
stream  bunting  rate  m/„  and  the  flame  thickness  8D: 


Tfc  =  T„  +  (Wc„)  =  Tfc  =  T 

(3.5) 

li 

Ol 

o? 

'"a. 

-< 

If 

(3.6) 

mi,  =  nt«(A„/ At)  =  mub{  Au/Afc) 

(3.7) 

So  =  (h/c,,)/m'u  =  8|j. 

(3.8) 

Equations  (3.5)  to  (3.8)  show  the  interesting  re¬ 
sults  that  the  important  parameters  T/„  m„  and  8|3 
of  this  adiabatic,  diffusionally-neutral,  unrestrained 
flame  arc  not  affected  by  stretch  and  therefore  as¬ 
sume  the  corresponding  values  of  the  one-dimen¬ 
sional  unstretched  flame.  The  downstream  burning 
rate  m i,  is  modified  front  the  unstretched  value 
m)',  by  the  area  ratio  (A„/A /,),  which  can  be  consid¬ 
ered  as  a  measure  of  the  extent  of  stretch  due  to 
effects  of  both  flow  divergence  and  thermal  expan¬ 
sion  ft  is  also  clear  that,  because  of  the  invariance 
of  the  flame  temperature,  extinction  cannot  be 
achieved  for  this  flame.  For  a  curved  flame  the  above 


results  remain  the  same  except  m„  is  modified  ac¬ 
cording  to  Eq.  (4.2) 

There  appear  to  lie  different  causes  leading  to 
the  various  earlier  interpretations  of  the  effect  ol 
stretch.  For  example,  the  argument  in  favor  ot  re¬ 
duced  flame  temperature  discusses  only  the  nega¬ 
tive  effect  of  stretch  on  flame  temperature  due  to 
the  diffusive  transport  of  heat  and  radicals  away  from 
the  reaction  zone,  but  it  fails  to  consider  the  pos¬ 
itive  effect  of  reactant  diffusion  into  the  reaction 
zone.  The  argument  based  on  reduced  flame  thick¬ 
ness  does  not  account  for  the  fact  that  the  flame, 
being  unrestrained,  can  move  freely  in  response  to 
changes  in  stretch.  There  are  also  interpretations 
based  on  some  nondinicnsional  and/or  transformed 
coordinates  which  however  depend  on  the  stretch 
rate. 

Asymptotic  analyses’1'1*'”  have  confirmed  that  T /, 
is  invariant  to  stretch  while  in/,  is  reduced  from  the 
unstrctched  value  in)',  by  a  fractional  amount  Ka. 
However,  the  discussions  of  Kefs.  29  and  30  on  the 
reduction  of  the  total  extent  ol  reaction  across  the 
tlame  per  unit  surface  area,  and  on  the  extinction 
behavior  of  such  a  flame,  have  overlooked  the  lacts 
that  flow  divergence  is  actually  the  major  cause  ol 
this  reduction,  and  that  extinction  is  brought  about 
with  only  a  small  amount  ol  incomplete  reaction  11 

To  further  substantiate  the  insensitivity  of  flame 
thickness  to  stretch,  the  temperature  profile  across 
an  almost  adiabatic,  diffusionally-neutral.  ethane/air 
flame  in  a  symmetrical  counterflow  has  been  ex¬ 
perimentally  measured.  ’1  In  Fig.  5  these  temper¬ 
ature  profiles  are  superimposed  by  shifting  their 
spatial  locations  such  that  the  upstream  boundaries 


/Air,  0  =  0.65 


Shifted  Spatial  Coordinate,  cm 

Fl<;  5  [Experimental  temperature  profiles  ol  ethauo/air  flames,  with  coordinate  shifted  to  demonstrate 
the  minimal  influence  of  stretch  on  the  flame  properties,  especially  the  flame  thickness,  lor  adiabatic, 
unrestrained,  diffusionally-neutral  flames.*1 
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ol  the  luminous  zones  coincide.  It  is  seen  that  these 
temperature  profiles,  and  thereby  the  tlaine  tem¬ 
peratures  and  thicknesses,  are  basically  unaffected 
by  stretch  The  thicknesses  of  the  luminous  zones 
are  also  not  affected.  It  is  significant  to  note  that 
at  the  state  of  extinction,  caused  by  incomplete  re¬ 
action  to  be  discussed  next,  the  stretch  rate  for  this 
mixture  has  been  determined  to  be  about  320  s'1 
Thus  the  flame  temperature  profile  remains  insen¬ 
sitive  to  stretch  up  to  stretch  rates  sufficiently  close 
to  that  of  extinction.  This  experimental  result  has 
also  been  corroborated  by  numerical  solutions. 11 

Using  D„  =  0.2  cnr/s  as  given  in  Ref.  32,  and 
S"  =  17.5  cni/s  determined  in  the  manner  of  Fig. 
16  which  will  be  discussed  later,  the  karlovitz 
number  fva  for  the  situations  of  Fig.  5  can  also  be 
estimated,  it  is  seen  that  these  values  are  all  suf¬ 
ficiently  smaller  than  unity,  specifically,  the  extinc¬ 
tion  value  is  only  lxarv  =  0.2. 

We  mention  in  passing  that  in  writing  the  overall 
energy  balance  of  Eq.  (3.2).  it  is  essential  that  the 
diffusive  flux  is  normal  to  the  reaction  zone  as  shown 
m  Fig  -fa.  It  the  diffusive  flux  is  forced  to  lie  along 
the  streamlines  iFig.  *11)'.  then  the  diffusion  term 
in  Ei).  (3.2'  is  suppressed,  resulting  in 

m„  =  in" tA(,/A„),  m/,  =  m)'„  (3.9) 

which  are  completely  contrary  to  Eqs.  (3.6)  and  (3.7). 
and  are  incorrect.  Thus  it  is  necessary  to  preserve 
the  orthogonality  relation  between  the  flame  sur- 
lace  and  the  diffusive  fluxes  when  performing  "quasi- 
one-dimcnsional"  analysis.  In  fact,  the  planar  flames 
m  the  divergent  flow  fields  of  Refs.  33  and  34  are 
stretch!  ess,  according  to  the  definition  of  K  in  Eq 
(2.2),  because  diffusive  transport  has  been  assumed 
to  occur  along  the  streamlines  (v  x  n  =  O'. 

Adiabatic  Restrained  Flames 

We  next  consider  situations  under  which  the  free 
movement  of  a  flame  in  response  to  stretch  is  re¬ 
strained,  as  typified  by  the  strongly-stretched  stag¬ 
nation  flame. 's  50  The  flame  is  now  situated  close 
to  the  stagnation  surface  and  thereby  loses  its  flex¬ 
ibility  to  move  freely.  If  the  stagnation  surface  is 
adiabatic  and  effectively  impenetrable  to  net  mass 
transfer,  then  with  further  increase  in  stretch  the 
residence  time  decreases,  the  flame  becomes  thin¬ 
ner,  complete  reaction  cannot  be  achieved,  and  ex¬ 
tinction  eventually  occurs.  Extinction  in  this  case  is 
caused  by  stretch-induced  incomplete  reaction,  not 
by  stretch  alone.  The  Damkohler  number  now  ex¬ 
erts  a  much  stronger  role  in  the  flame  behavior, 
especially  the  state  of  extinction. 

If  the  stagnation  surface  is  penetrable,  as  in  the 
case  of  impinging  the  combustible  mixture  against 
a  hot  product  gas  stream,  then  the  reaction  zone 
can  actually  migrate  across  the  stagnation  surface 


and  thereby  assume  a  negative  flame  speed.'1’ The 
fuel  consumption  and  heat  release  rates  are  ol  course 
still  positive.  Combustion  m  this  situation  is  sup¬ 
ported  by  diffusion  of  the  reactants  against  the  un¬ 
favorable  convection  The  existence  of  negative  flame 
speeds  has  been  experimentally  verified’’  for  su- 
peradiabatie  flames  whose  downstream  temperature 
is  higher  than  the  adiabatic  flame  temperature. 
Adiabatic  negative  flame  speeds  have  not  been  ob¬ 
served. 

Another  example  ol  a  restrained  flame  is  the  tip 
segment  of  a  Bunsen  flame  because  the  orientation 
of  the  shoulder  segment  is  somewhat  fixed  for  a 
given  mixture  flow,  and  because  the  upstream  flame 
speed  at  the  tip  is  also  fixed  at  the  value  of  the 
local  (low  velocity  Thus  the  flame  speed  here  is  a 
given  parameter  instead  of  a  flame  response  to  hi’ 
solved.  Indeed,  the  exceptionally  large  value  ol  the 
"flame  speed"  at  the  flame  tip2  has  long  been  an 
anomaly  in  combustion  science  Some  preliminarv 
theoretical  and  experimental  investigations’1'  "  have 
been  conducted,  although  the  detailed  structure  ol 
this  restrained,  highly -curved  flame  segment  in¬ 
quires  further  study 

Sonadiabatic  Flames: 

The  effects  of  heat  loss  are  quite  obvious,  it  slows 
down  the  burning  rate  and  can  eventually  cause  ex¬ 
tinction.  Aerodynamic  stretching  changes  the  prox¬ 
imity  of  a  flame  to  the  heat  sink  while  tlamc  cur¬ 
vature  also  changes  the  intensity  of  the  heat  flux 
over  the  llamc  surface.  Thus  the  influence  ol  stretch 
and  heat  loss  are  intimately  coupled.  Libby  and 
Williams'*1'40  have  performed  detailed  analysis  of 
uonadiabatic  counterflow  flames  Experimentally. 
I -aw  et  al."  have  found  that  downstream  heat  loss 
does  not  strongly  affect  the  flame  behavior  as  long 
as  the  (lame  is  not  too  close  to  the  cold  stagnation 
surface. 


Flame  Stretch  in  Diflusionally-lmbalanced 
Mixtures 

/’/leimmewi/ogi/ 

While  pure  stretch  basically  has  minimal  effect 
on  the  propagation  of  an  adiabatic,  dillusiouallv  - 
neutral,  unrestrained  flame  the  flame  response  ran 
be  significantly  modified  in  the  presence  ol  picl- 
erential  diffusion  To  demonstrate  this  influence  we 
first  note  dial  there  are  at  least  three  diffusiv  ities 
of  interest  lor  an  inert-abundant  mixture,  nainclv 
those  associated  with  heat  (1)T).  the  deficient  reac¬ 
tant  S 0,1.  and  the  excess  reactant  (D().  From  these 
three  diffusivities  two  interpretations  for  the  effects 
of  preferential  diffusion  have  been  developed,  com¬ 
paring  I),  with  1>|.  and  D,  with  D,.  lor  sufficiently 
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ott-  and  near-stoichiometric  situations  respectively. 
These  two  interpretations  can  be  respectively  termed 
non-unity  Lewis  number  (Le  =  DT/D,)  effect  and 
differential  diffusion  (D,/D;)  effect. 

We  again  take  the  stagnation  flow  as  an  example 
(Fig  6a),  and  draw  a  control  volume  enclosing  the 
transport  zone  and  the  divergent  streamline  as 
shown  Then  with  the  Lewis  number  interpreta¬ 
tion,12'13  it  is  clear  that  the  control  volume  loses 
thermal  energy  to  the  external  streamlines  while  it 
gains  chemical  energy  from  them  due  to  an  in¬ 
crease  of  the  deficient  reactant  concentration.  Thus 
the  flame  behavior,  especially  the  flame  tempera¬ 
ture,  depends  on  the  relative  rates  of  heat  and  mass 
diffusion.  If  the  diflusivities  are  equal  such  that  Lc, 
=  1,  then  total  energy  conservation  is  maintained 
and  the  flame  temperature  is  the  adiabatic  flame 
temperature.  This  is  the  case  studied  in  the  pre¬ 
vious  section.  However,  if  Le  >  1,  heat  loss  ex¬ 
ceeds  mass  gain  and  we  expect  T/,  <  T„,/.  Con¬ 
versely  T/,  >  T if  Le  <  1. 

We  next  take  the  differentia)  diffusion  interpre¬ 
tation.  Then  if  the  leaner  reactant  is  also  the  more 
diffusive  one,  the  reactant  concentration  in  the  re¬ 
action  zone  will  become  more  stoichiometric  such 
that  the  flame  temperature  is  higher  and  the  burn¬ 
ing  more  intense.  The  converse  holds  if  the  leaner 
reactant  is  less  diffusive.  In  the  following  we  shall 
adopt  the  Lewis  number  interpretation  because  off- 
stoichiometric  burning  is  more  prevalent  and  be- 
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Ft<;  6  Schematics  of  (a)  the  stagnation  flow,  and 
ibl  the  Bunsen  flame,  illustrating  the  coupled  influ¬ 
ence  of  flame  stretch  and  preferential  diffusion  on 
the  flame  response. 


cause  heat  diffusion  is  a  crucial  mechanism  in  flame 
propagation.  Experimental  results44  also  indicate  that 
the  Lewis  number  interpretation  seems  to  have  a 
wider  range  of  applicability. 

If  we  now  attempt  to  extinguish  this  stagnation 
flame  by  increasing  the  stretch  rate,  the  flame  will 
be  pushed  closer  to  the  stagnation  surface  At  the 
same  time  it  will  suffer  stronger  flame  stretch  and 
thereby  preferential  diffusion  effects.  Thus  it  is  clear 
that  for  Le  >  1.  there  exists  a  critical  stretch  rate 
at  which  I/,  will  be  reduced  to  such  an  extent  that 
steady  burning  is  not  possible.  Extinction  occurs 
when  the  flame  is  at  a  finite  distance  awax  from  the 
surface,  with  the  lean  reactant  almost  completed) 
depleted  in  crossing  the  flame.  On  tile  other  hand, 
for  a  Le  <  1  flame,  increasing  stretch  elevates  T(, 
and  therefore  extinction  cannot  occur  until  the 
downstream  boundary  ol  the  reaction  zone  is  pushed 
onto  the  stagnation  surface.  With  further  stretch¬ 
ing.  reaction  cannot  be  completed  because  ol  the 
reduced  residence  time.  Onlx  then  will  the  flame 
temperature  start  to  decrease,  leading  cxcntualK  to 
extinction.  Thus  the  extinction  mechanism  is  simi¬ 
lar  to  that  ol  the  Le  =  1  case  studied  prcxiouslx. 
except  the  extinction  stretch  rate  with  Lc  <  1  is 
higher  in  order  to  compensate  for  the  higher  flame 
temperature  initially  attained. 

The  second  situation  illustrating  flame  stretch  ef¬ 
fects  involves  the  burning  intensity  over  tile  curved 
surface  ol  an  axisvmmetric  Bunsen  cone  (Fig.  fihi. 
Here  if  we  assume  for  simplicity  that  the  lloxx  is 
uniform,  then  flame  stretch  is  mainly  manifested 
through  curvature  effects,  especially  in  the  tip  re¬ 
gion  which  has  the  strongest  curvature.  Thus  lor  a 
closed  tip,  its  concave  nature  towards  the  bosh 
mixture  focuses  the  heat  ahead  of  the  flame  and 
therefore  tends  to  raise  the  flame  temperature  to  a 
value  above  T On  the  other  hand,  this  curvature 
has  a  dcfocusiiig  effect  on  the  reactants  approaching 
the  flame,  and  therefore  tends  to  reduce  the  flame 
temperature.  Thus  the  temperature  of  the  flame 
again  depends  on  the  relative  rates  of  heat  and  mass 
diffusion.  The  deviation  increases  progressively  from 
the  flame  base  towards  the  flame  tip  because  ol  the 
corresponding  increase  in  the  stretch  intensitx  Thus 
for  Le  >  1.  the  tip  will  burn  more  mteusclx  rel¬ 
ative  to  the  shoulder  region  of  the  flame,  xxhile  loi 
Le  <  I  the  horning  is  less  intense  and  can  lead  to 
local  extinction,  eommonlv  known  as  opening,  ol  the 
tip.  Note  that  because  ol  the  compressive  nature  ol 
the  Bunsen  flame,  its  behavior  in  response  to  Lewis 
number  variations  is  completely  opposite  to  that  of 
the  positively-stretched  stagnation  flame. 

Energy  Conservation ; 

The  above  result  that  the  temperature  of  the  burnt 
gas  ol  a  diflusionallx -imbalanced  stretched  flame  can 
deviate  from  the  adiabatic  flame  temperature  of  its 
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unburnt  mixture  appears  to  violate  the  principle  of 
energy  conservation.  To  explore  this  point,  Figs.  7a 
and  7b  show  the  profiles  of  the  stoichiomctrieally- 
weighted  temperature  T.  mass  fraction  Y,  and  the 
total  energy  as  represented  by  the  coupling  func¬ 
tion  P  =  T  +  Y,  for  the  ease  of  Le  <  1  for  which 
i'/»  >  T «,,/  where  T„,y  =  T . ,  4-  Y  . ...  It  is  clear  that 
in  a  steady-state  situation,  as  for  the  stagnation  and 
Bunsen  Haines,  energy  is  not  conserved  as  the  fresh 
mixture  traverses  the  flame.  This  is  a  local  phe¬ 
nomenon  because  only  diffusion  normal  to  the  flame 
has  been  considered.  It  is  reasonable  to  expect  that 
global  energy  conservation  can  again  be  established 
by  considering  diffusion  of  the  excess/ deficit  heat 
or  mass  in  the  tangential  direction  along  the  flame 
within  the  transport  zone  of  the  entire  flame  sur- 
face. 

The  consideration  is  somewhat  different  for  the 
expanding  spherical  flame  because  the  entire  flame 
surface  is  already  accounted  for.  The  profiles  of  Fig. 
7  are  now  instantaneous  values.  Energy  in  this  case 
is  conserved  in  that  the  increase  in  energy  of  the 
burnt  gas  over  T,,,;  is  achieved  by  drawing  excess 
energy  in  terms  of  the  reactant  concentration  from 
the  upstream  mixture.  As  the  flame  expands  the 
amount  of  excess  energy  associated  with  the  burnt 
gas  increases  because  of  the  increase  in  volume. 
However,  at  the  same  time  the  amount  of  energy 


Fit;.  7  .Schematics  of  (a)  the  temperature  and 
concentration  profiles,  and  (b)  total  energy  profile 
lor  a  Le  <  1  positively  stretched  flame,  illustrating 
the  discussions  on  energy  conservation. 


deficit  in  the  transport  zone  also  increases  because 
of  the  corresponding  increase  in  the  volume  of  the 
transport  zone  Thus  at  any  instant  the  increase  in 
excess  energy  in  the  burnt  zone  is  achieved  at  the 
expense  of  the  corresponding  decrease  in  energy  in 
the  transport  zones.  Consequently  the  amount  ol 
excess  and  deficit  energies  exactly  balance. 


Theoretical  Results: 

Most  analytical  studies  on  stretched  flames  arc 
based  on  large  activation  energy  asymptotics  These 
sin'  s  can  be  grouped  into  three  categories  namely 
'.  ulti-scale  analy  ses  of  a  general  wrinkled  flame 
in  a  nomuutorm  flow,’1'1  ‘t‘  (l>)  detailed  analyses  ol 
a  specific-  flame  configuration  such  as  the  stagnation 
flame. ,IM- 4ri  ,H  ,s  the  Bunsen  flame.’1’  "  h~  the 
spherical  flame. h’*  ,1,s  or  the  tubular  flame.'"  "’  and 
ye)  analyses  involving  ttame/ilow  situations  charac¬ 
terized  by  some  prescribed  mathematical  forms  ol 
stretch  "'1  "The  first  category,  while  general  in  tin- 
physical  situations,  yields  only  linearized  results  for 
weak  stretch.  The  second  category,  while  being 
specific  in  the  physical  configuration,  does  predict 
nonlinear  phenomena  such  as  extinction.  The  third 
category  can  also  he  very  uselul  in  identify  ing  some 
specific-  features  of  stretched  flames.  Recently  the 
method  of  integral  analysis  has  also  been  ap¬ 
plied,2,  28  yielding  characteristics  of  the  bulk  flame 
parameters  which  appear  to  be  both  general  and 
capable  of  exhibiting  the  nonlinear  effects  of  strong 
stretch  and  flame  extinction . 

The  linearized  solution  ol  Ref.  27  yields  the  lol- 
lowing  explicit  expressions  of  the  various  flame  re¬ 
sponses  for  a  wrinkled  flame  in  a  nomiiiiform  flow- 
field: 


Ti,  T,, 

—  —  - - 1  +  (  —  —  1  1  Ka 

•n:  t 


Le 


,4  L 


\m„/  S„ 


V  •  n 


Le 


Ixa 


m/,  S,  ka 

— - - ;  =  1 - +  I  —  -  I 

in/,  Si,  Le  \Lc 


<2Tu,//T„ 

ka 

C2T„,,/Tu> 


,4  2: 


,4.31 


where 


Ka  =  —  K  = 
VS” 


-  JPlsL 

(S”)2 


5  K 


IX 


[S”(Le  =  1)1- 


5  K  ,44' 


B-9 


INVITED  LECTURE 


Fic.  8.  Schematic  showing  the  response  of  the  reaction  rate  to  increasing  stretch  in  a  stagnation  flow, 
for  various  Lewis  numbers.71 


\  S"  /  l  S"(Le  =  1)  J 


and  we  have  used  (he  relation  S"(Le)  = 
VLeS"(l). 

The  above  results  show  that  the  (lame  temper¬ 
ature  can  deviate  from  the  adiabatic  flame  temper¬ 
ature  only  iu  the  simultaneous  presence  of  stretch 
(Ka  /  0)  and  preferential  diffusion  (Le  #  1).  The 
modification  in  the  downstream  velocity  results  from 
flow  divergence  and  changes  in  the  flame  temper¬ 
ature  and  thereby  burning  intensity  in  the  presence 
of  preferential  diffusion,  as  given  by  the  second  and 
third  terms  of  Eq.  (4.3).  The  modification  iu  the 
upstream  velocity  results  from  the  flame  curvature 
and  the  preferential  diffusion  effects.  These  results 
essentially  degenerate  to  Eqs.  (3.5)  to  (3.7)  for  a 
planar,  Le  =  1  flame. 

Equation  (4.4)  shows  that  the  effect  of  stretch  as 
represented  by  Ka  can  be  increased  by  increasing 
the  stretch  rate  or  mixture  diffusivity,  or  by  de¬ 
creasing  the  flame  speed.  We  also  note  that  be¬ 
cause  the  laminar  flame  speed  S l  increases  with  in¬ 
creasing  Lewis  number,  the  dependence  of  S/,  on 
the  mixture  diffusivity  can  be  actually  non-mono¬ 
tonic.  It  further  implies  that  the  preheat  zone 
thickness  5:J  depends  on  the  geometric  average  of 
the  thermal  and  mass  diffusivities,  from  which  a 


weighted  dilfusivitv  L>"  can  be  defined,  as  in  Eq. 
(4.5). 

A  generalized  nonlinear  theory  of  wrinkled  flames 
based  on  differential  analysis  has  not  been  formu¬ 
lated,  although  an  integral  theory  has  been  ad¬ 
vanced.2*  The  salient  features  of  the  flame  re¬ 
sponse,  however,  can  best  be  illustrated  b\  the 
numerical  solution  of  the  simple  adiabatic  stagna¬ 
tion  flame.'3  '1  Figure  8  shows  the  qualitative  be¬ 
havior  of  the  reaction  rate  profiles  with  increasing 
stretch  for  various  Lewis  numbers.  1  For  Le  >> 
1,  increasing  stretch  decreases  the  maximum  reac¬ 
tion  rate  and  extinction  occurs  when  the  reaction 
zone  is  away  from  the  stagnation  surface  and  re¬ 
action  is  complete.  For  Le  >  1,  reaction  rate  still 
decreases  with  stretch  except  now  extinction  occurs 
with  the  reaction  zone  pressed  against  the  surface 
and  reaction  is  not  complete  For  Le  =  1.  the  re¬ 
action  rate  remains  unchanged  until  the  reaction 
zone  reaches  the  surface.  Finally,  for  Le  <  1,  the 
reaction  rate  continuously  increases  with  stretch  until 
the  reaction  zone  reaches  the  surface,  then  it  de¬ 
creases.  Thus  except  for  the  case  of  Le  >>  1.  in¬ 
complete  reaction  is  essential  in  causing  flame  ex¬ 
tinction. 

Figure  9  shows  variation  of  the  maximum  flame 
temperature  with  stretch  for  different  Le.  The  be¬ 
havior  basically  corroborates  that  of  the  reaction  rate 
profiles. 

The  influence  of  reaction  intermediates  on  the 
response  ol  stretched  flames  has  not  been  ade¬ 
quately  studied.1  ’  ,'1  Numerical  solutions  of  the 
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Fic.  9.  Schematic  of  the  variation  of  the  maxi¬ 
mum  flame  temperature  in  response  to  increasing 
stretch,  for  various  Lewis  numbers. 


stagnation  flame  with  complex  chemistry'5  ' 7  have 
shown  that  the  behavior  of  the  bulk  flame  param¬ 
eters  are  in  qualitative  agreement  with  interpreta¬ 
tions  based  on  a  single  reactant,  without  consid¬ 
ering  the  possible  preferential  diffusion  nature  of 
the  combustion  intermediates.  Reference  76,  how¬ 
ever,  has  noted  the  role  of  incomplete  reaction  of 
the  II  radical  in  flame  extinction,  while  the  asymp¬ 
totic  results'8  obtained  with  a  two-step  branching- 
termination  chain  reaction  model  have  shown  that 
depending  on  the  parameter  range  the  mass  burn¬ 
ing  rate  can  he  either  increased  or  decreased. 

F.xpcrimcnUil  Results: 

The  above  discussions  show  that  the  response  of 
a  stretched  flame  exhibits  opposite  behavior  when 
the  stretch  changes  from  positive  to  negative,  and 
when  the  mixture’s  effective  Lewis  number  is  greater 
or  less  than  a  critical  value,  which  is  unity  for  the 
flame  temperature.  These  completely  opposite  trends 
should  provide  definitive  verification  of  the  concept 
of  flame  stretch  with  preferential  diffusion. 

Two  groups  of  mixtures  are  especially  suitable  for 
the  study  of  preferential  diffusion  effects.  The  first 
group  consists  of  lean  inethanc/air  and  rich  pro¬ 
pane/air  mixtures.  Rough  estimates32  show  that  their 
effective  Lewis  numbers,  based  on  the  deficient 
species,  are  less  than  unity.  Thus  positive  (nega¬ 
tive)  stretch  is  expected  to  increase  (decrease)  the 
flame  intensity  of  such  mixtures.  This  argument  still 
holds  even  if  we  just  consider  the  relative  mass  dif- 
fusivities  of  the  fuel  and  oxidizer  species.  That  is, 
based  on  molecular  weight  considerations,  we  ex¬ 
pect  the  diffusivities  of  the  various  reactants  rela¬ 
tive  to  nitrogen  should  increase  in  the  order  of  pro¬ 
pane,  oxygen,  and  methane.  Thus  positive  (negative) 
stretch  will  increase  (decrease)  the.  methane  con¬ 
centration  of  a  lean  methane/air  mixture  but  de¬ 
crease  (increase)  the  propane  concentration  of  a  rich 
propane/air  mixture  at  the  flame.  Both  mixtures  are 


consequently  rendered  more  (less)  stoichiometric  at 
the  flame,  leading  to  enhanced  (reduced I  binning 
intensity 

The  second  group  consists  of  rich  methane/air 
and  lean  propane/air  mixtures,  whose  effective  Ix-wis 
numbers  are  greater  than  unity  while  positive  (neg¬ 
ative)  stretch  also  renders  the  mixture  less  (morel 
stoichiometric.  Thus  the  responses  of  these  two 
groups  of  mixtures  to  stretch  are  expeetei.  to  he 
qualitatively  opposite. 

Three  series  of  experiments  have  been  con¬ 
ducted  using  these  mixtures  in  flame  configurations 
exhibiting  positive  and  negative  stretches  For  pos¬ 
itive  stretch,  extensive  experiments  have  been  per¬ 
formed  by  using  the  symmetrical  counterflow  s\s- 
tem. 32  80  81  This  system  has  the  advantage  over  the 
one  created  by  impinging  a  single  stream  onto  a 
solid  surface  in  that,  because  of  symmetry,  the 
stagnation  surface  in  a  counterflow  is  close  to  being 
adiabatic. 

Although  theoretical  results  show  that  both  the 
flame  temperature  and  flame  speed  exhibit  flame 
stretch  and  non-unity  Lewis  number  effects,  ex¬ 
perimental  determination  of  the  flame  temperature 
is  more  definitive  as  compared  with  that  for  the 
flame  speed,  whose  value  depends  on  where  it  is 
evaluated.  Thus  in  Fig.  10  we  show  the  measured 
flame  temperature  T|  and  the  distance  of  the  center 
of  the  luminous  zone  yj  from  the  stagnation  surface 
with  varying  stretch  rate.  If  we  exclude  the  small 
regimes  near  extinction,  the  results  agree  with  pre¬ 
vious  discussions  in  that  with  increasing  stretch  the 
llame  temperature  increases  for  a  lean  methane/air 
flame  and  a  rich  propane/air  flame,  decreases  for 
a  lean  propane/air  flame,  and  slightly  decreses  for 
a  rich  methane/air  flame.  Furthermore,  the  flame 
extinguishes  when  it  is  away  from  the  stagnation 
surface  for  rich  methane/air  and  lean  propane/air 
flames,  and  at  the  stagnation  surface  tor  lean  n icth- 
ane/air  and  rich  propane/air  flames. 

Next  we  examine  the  flame  temperature  re¬ 
sponse  to  negative  stretch,  provided  b>  the  increas¬ 
ing  curvature  along  the  surface  of  an  axisymmetrie 
Bunsen  cone.44  84  85  The  behavior  should  he  com¬ 
pletely  opposite  to  that  of  the  counterflow  flame. 
Figure  11  shows  the  measured  maximum  flame 
temperature  Ty-  along  the  flame  surface.  Excluding 
the  segment  near  the  flame  base  where  burning  is 
weak  due  to  heat  loss  to  the  burner  rim,  it  is  seen 
that  as  the  intensity  of  negative  stretch  increases  by 
moving  along  the  flame  towards  its  tip,  the  flame 
temperature  increases  for  the  rich  methane/air  and 
lean  propane/air  flames,  but  decreases  for  the  lean 
methane/air  and  rich  propane/air  flames.  The  neu¬ 
tral  compositions  are  approximately  <f>  =  1.00  and 
0.94  for  methane/air  and  propane/air  Haines  re¬ 
spectively.  Figure  11  further  shows  that  ethylene, 
whose  molecular  weight  is  intermediate  between 
those  of  methane  and  propane,  shows  almost  neu- 
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Fl<;.  11.  Experimental  results  showing  the  quali¬ 
tatively  different  behavior  for  the  flame  tempera¬ 
ture  of  negatively-stretched  axisymmetric  Bunsen 
flame  for  mixtures  exhibiting  (CHj/Air,  C.,H„/Air) 
and  not  exhibiting  (C2H«/Air)  preferential  diffusion 
effects.*4 


FlC.  12.  Experimental  results  on  the  ii..;nr  »•  m- 
perature  of  two-dimensional  Bunsen  flames  .  1 
methane/air,  and  (b)  propane/air  mixtures,  show¬ 
ing  the  lack  of  influence  of  preferential  diffusion  on 
stretchless  Haines.44 


Eq.  (5.2)  is  essentially  the  criterion  (or  blowoff  (of 
inverted  Haines)  as  originally  advanced  by  Karlovitz1 
and  Lewis  and  von  Elbe.2  This  stabilization  mech¬ 
anism  is  dynamic  in  nature,  and  strictly  speaking  is 
also  not  a  stretch  phenomenon  because  only  the 
normal  velocity  gradient  is  considered.  Indirectly, 
stretch  may  appear  to  influence  flame  stabilization 
because  of  the  relationship  between  the  normal  and 
transverse  velocity  gradients,  and  because  S„  can 
be  affected  by  flame  stretch  through  its  coupled  ef¬ 
fects  with  heat  loss,  flow  nonuniformity,  flame  cur¬ 
vature,  and  preferential  diffusion. 

The  above  stabilization  mechanism  is  physically 
sound.  However,  in  an  effort  to  be  quantitatively 
predictive,  Lewis  and  von  Elbe2  further  proposed 
that  blowoff  should  occur  when  the  velocity  varies 
significantly  within  the  flame,  or 


The  validity  of  Eq.  (5.3),  or  similar  expressions 
with  different  forms  of  velocity  gradients,  has  been 
extensively  explored.2'92'98  However,  the  meaning¬ 
fulness  of  such  a  criterion  and  the  fidelity  of  the 
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Rf ,  mm 


FlC.  13.  Experimental  results  of  the  flame  speeds 
of  the  positively-stretched  outwardly-propagating 
flame  for  propane/air  mixtures  exhibiting  prefer¬ 
ential  diffusion  effects.86  87 
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(b) 

FlC.  14.  Extinction  stretch  rates  Kr,  for  the  sym¬ 
metrical  counterflow  methane/air  and  propane/air 
flames.88 


comparisons  have  been  seriously  questioned.81*1 
Ka warn  ura  and  co-workers91100  have  also  shown  that 
flame  curvature,  instead  of  flow  gradient  as  given 
by  Eq.  (5.3),  yields  better  correlation  with  the 
blowoff  data  of  inverted  flames. 


In  light  of  our  earlier  discussion,  the  possible 
failure  of  Eq.  (5.3)  can  be  appreciated  for  several 
reasons.  First,  the  original  postulate2  of  Eq.  (5.3) 
is  based  on  the  concept  that  the  flame  temperature 
is  reduced  with  increasing  stretch,  which  eventu- 
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.ilh  causes  extinction  We  have,  however,  shown 
that  the  flame  temperature  can  either  increase  or 
decrease  with  stretch  depending  on  the  nature  of 
preferential  diffusion  Second,  Eq.  (5.3)  also  does 
not  account  for  the  influence  of  heat  loss  and  flame 
curvature  on  the  flame  speed.  Third,  the  inverted 
flame  does  not  appear  to  be  sufficiently  restrained 
to  suffer  extensive  incomplete  reaction.  It  is  there¬ 
fore  quite  likely  that  Eq,  (5.3)  is  not  applicable  to 
flame  blowoff. 

Determination  of  Laminar  Flame  Speeds. 

In  spite  of  the  technological  and  fundamental  im¬ 
portance  of  the  laminar  flame  speeds  S"  of  com¬ 
bustible  mixtures,  and  the  extensive  efforts  ex¬ 
pended  to  accurately  determine  their  values, 
substantial  systematic  spreads  still  exist  in  the  re¬ 
ported  data.  While  it  is  not  always  easy  to  individ¬ 
ually  identify  the  causes  for  such  discrepancies,  it 
is  recognized  that  since  most  of  these  determina¬ 
tions  involve  stretched  flames  such  as  the  Bunsen 
and  outwardly-propagating  spherical  flames,  the  re¬ 
ported  data  are  actually  S„  instead  of  S".  A  partic¬ 
ularly  serious  implication  of  the  influence  of  stretch 
is  the  potential  inaccuracy  in  the  kinetic  informa¬ 
tion  determined  or  validated  through  comparison 


Ka 


Kit;  lb  Principle  of  counterflow  mcthodolog)  m 
determining  S”:  (at  Velocity  profile  and  definition  of 
reference  S„  lb)  extrapolation  of  reference  S„  to  yield 

s;: 


between  the  numerically-calculated  S"  and  experi¬ 
mentally-determined  S,,.1111  111,1 

Recently  the  symmetrical  cotmterflow  flame  con¬ 
figuration  has  been  used  through  which  stretch  cl- 
fects  can  be  systematically  subtracted  out.'’9  11,1  ll,s 
The  experiment  basically  involves  determining  the 
axial  velocity  profile  along  the  centerline  of  the  flow 
by  using  LDV  (Fig.  16a),  and  identifying  the  min¬ 
imum  point  of  the  velocity  profile  as  a  reference 
upstream  flame  speed  S„  corresponding  to  the  im¬ 
posed  stretch  rate  1C.  Thus  by  plotting  S„  versus 
K,  S','  can  be  unambiguously  determined  through 
linear  extrapolation  to  K  =  0  (Fig  16b). 

Figure  17  shows  close  agreement  between  the 
experimentally-determined  S"  and  the  numerically- 
calculated  values10,  'IHS  by  using  the  flame  code  ol 
Kec  et  al.109  Such  close  agreement  oxer  extended 
ranges  of  equivalence  ratio  and  pressure  provides 
mutual  support  for  the  validity  of  the  experimental 
methodology  and  the  kinetic  scheme,  as  lar  as  flame 
speed  determination  is  concerned. 

An  accurate  determination  of  S"  also  allows  an 
assessment  ol  the  overall  reaction  order  n  because 
S"  ~  p'"/2)~  .  Figure  IS  shows  experimentally-  and 
numerically-determined  values  of  n  as  a  function  ol 
pressure  lor  stoichiometric  methane/air  flames  with 
various  amount  ol  nitrogen  dilution  1111  It  is  seen 
that  n  decreases  with  increasing  pressure  and  ni¬ 
trogen  dilution  such  that  it  can  actually  assume 
negative  values.  This  behavior  is  caused  by  the  in¬ 
creased  relative  importance  ol  the  three-body  ter¬ 
mination  reactions  over  the  two-body  branching  re¬ 
actions  with  increasing  pressure  and  decreasing  flame 
temperature. 

Concept  of  Flammability  Limits: 

Although  the  term  "flammability  limit ”  has  been 
widely  and  loosely  applied  to  diverse  situations  of 
unsustainable  combustion.9,  !!fs  111  121  the  general 


Flo  17  Comparison  between  experiment. >llv -de¬ 
termined  S‘t)  bv  using  the  coimtcrflovv  technique  and 
numerically-calculated  values  bv  using  the  lull  C- 
kinetics  ol  VVarnatz  lor  methane/air  flames 
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Fic.  18.  Experimentally-  and  numerically-deter¬ 
mined  values  of  the  overall  reaction  order  n  for 
stoichiometric  methane/air  flames  with  various 
amount  of  nitrogen  dilution,  demonstrating  the  de¬ 
creasing  trend  of  n  with  increasing  pressure  and  the 
existence  of  negative  values  of  n. 

consensus28,  is  that  the  flammability  limit  of  a 
combustible  mixture  should  be  a  unique  physico¬ 
chemical  property  of  that  mixture,  being  indepen¬ 
dent  of  external  system  parameters.  As  such,  at  the 
flammability  limit  all  external  extinction  mecha¬ 
nisms  such  as  conductive  heat  loss  and  radical  loss 
should  have  been  eliminated  such  that  it  is  truly 
the  ultimate  concentration  limit  beyond  which 
combustion  is  absolutely  not  possible. 

As  of  now  it  is  not  clear  what  is  the  role  of  stretch 
in  flammability  limits.  Certainly  if  the  earth's  grav¬ 
ity  is  counted  as  an  adjunct  property  of  the  mix¬ 
ture,  then  buoyancy  induced  flame  curvature  and 
thereby  stretch  can  extinguish  a  flame  even  in  the 
absence  of  heat  loss.97-98  119  Furthermore  it  is  con¬ 
ceivable  that  curvature  can  develop  over  the  sur¬ 
face  of  a  flame  in  the  course  of  its  resistance  to  ex¬ 
tinguishment.  Extinction  may  eventually  occur  when 
the  curvature  attains  a  certain  critical  value. 

One  conceptual  difficulty  with  including  stretch 
as  a  viable  flammability  mechanism  is  the  influence 
of  preferential  diffusion.  That  is,  depending  on  the 
nature  ol  the  stretch  and  the  mixture's  Lewis  num¬ 
ber,  flame  stretch  can  either  extinguish  a  flame  or 
support  the  propagation  of  a  sub-limit  mixture.  An 
interesting  example  is  the  phenomenon12'*  of  "self¬ 
extinguishing  flame,  in  which  an  outwardly-prop¬ 
agating  flame  can  be  initially  established  in  a  sub¬ 
limit  Le  <  1  mixture  because  of  the  enhanced 
burning  due  to  preferential  diffusion.  Extinction 
eventually  occurs  as  the  flame  size  increases  and 
the  stretch  effect  decreases. 

Recently  the  influence  of  stretch  on  the  extinc¬ 
tion  of  an  adiabatic  flame  has  been  experimentally 
determined88  by  first  measuring  the  extinction 
stretch  rate  of  the  symmetrical  counterflow  flame  as 
a  function  of  concentration  (Fig.  14),  and  extrapo¬ 
lating  the  data  to  zero  stretch  to  yield  the  corre¬ 


sponding  concentration.  The  values  so  determined 
for  lean  and  rich  methane/air  and  propane/air  mix¬ 
tures  lie  within  the  previously  reported  range  of 
flammability  limits  containing  the  effect  of 
stretch. 112  125  Since  this  range  is  quite  narrow,  these 
results  suggest  that  the  influence  of  stretch,  even 
if  it  does  exist,  is  probably  very  weak  as  compared 
to  other  internal  extinction  mechanisms. 

There  are  two  possible  internal  extinction  mech¬ 
anisms  which  can  cause  the  failure  of  the  stretch¬ 
less  planar  flame  to  propagate,  namely  radiative  heat 
loss12'*  and  the  emergent  importance  of  chain  ter¬ 
mination  over  chain  branching  with  the  gradual  re¬ 
duction  of  the  flame  temperature  as  the  flamma¬ 
bility  limit  is  approached.  Theories  of  flammability 
have  been  formulated  based  on  heat  loss121’  12 '  and 
chain  mechanism,120  128  with  the  theory  of  Ref.  128 
appearing  to  be  quite  quantitatively  predictive 


Flame-Front  I  instability: 

When  a  smooth  flamefront  is  subjected  to  non- 
planar  disturbances  the  resulting  wrinkles  can  either 
grow  or  decay.  This  is  a  stretched  flame  phenom¬ 
enon  because  the  susceptibility  to  instability  de¬ 
pends  on  the  flame  response  to  the  development  of 
the  flame  curvature  over  its  surface. 

There  are  two  intrinsic  modes  of  flamefront  in¬ 
stability.  The  first  is  the  Landau,  or  hydrodynamic 
instability,129  which  has  its  origin  from  thermal  ex¬ 
pansion  of  the  gas  upon  crossing  the  flame.  The 
flame  is  assumed  to  be  structureless,  locally  prop¬ 
agating  with  the  constant  laminar  flame  speed 
(Fig.  19a).  The  gas  densities  assume  constant  but 
different  values  upstream  and  downstream  of  the 
flame  sheet.  It  can  then  be  easily  argued8  that  upon 
being  wrinkled  the  flow  ahead  of  the  convex  seg¬ 
ment  slows  down  because  of  widening  of  the 
streamtube.  However,  since  the  flame  speed  re¬ 
mains  unchanged,  a  dynamic  imbalance  results  which 
leads  to  a  further  protrusion  of  this  flame  segment 
Similar  argument  for  the  concave  segment  shows 
that  it  will  further  recede  into  the  burned  mixture. 
Thus  the  flame  is  absolutely  unstable  according  to 
this  mechanism,  which  involves  only  hydrodynamic 
stretch  without  considering  the  effect  of  flume  stretch 
within  the  flame  structure. 

The  importance  of  ditfusional  transport  on  flame 
stability  can  be  readily  appreciated  by  considering 
the  effect  of  heat  and  mass  diffusion  on  the  ability 
of  the  flame  to  restore  stability  (Fig.  19ln.  Follow¬ 
ing  the  same  discussion  as  that  for  the  intensifica¬ 
tion  and  weakening  of  the  Bunsen  flame  tips,  it  is 
clear  that  a  Le  >  1  flame  is  stable  in  that  the  burn¬ 
ing  is  more  (less)  intense  in  the  receding  (protrud¬ 
ing)  segment,  which  tends  to  advance  (recede)  to 
resume  its  undisturbed  position.  Similar  argument 
shows  that  a  Le  <  1  flame  is  unstable.  Experi- 
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(Le  <  1) 


Fic.  19.  Schematics  demonstrating  mechanisms  of 
(a)  hydrodynamic  instability  and  (b)  dilfusional-ther- 
mal  instability. 


mental  evidence  of  the  above  phenomena  can  be 
found  in  Ref.  3. 

Sivashinskv5  130  has  analyzed  the  diffusional-ther- 
inal  instability  by  neglecting  thermal  expansion. 
Analyses  allowing  for  both  hydrodynamic  and  dif- 
fusional-thermal  instabilities  have  also  been  ad¬ 
vanced  by  Sivashinskv131  assuming  small  heat  re¬ 
lease.  and  by  Pelce  and  Clavin132  and  Matalon  and 
Matkowskv133  through  multi-scale  analyses.  These 
analyses  show  that  flames  are  hydrodynamically  un¬ 
stable  when  subjected  to  long  wavelength  distur¬ 
bances.  At  short  wavelengths  the  diflusional-ther- 
mal  mechanism  can  either  stabilize  or  further 
destabilize  the  flame,  depending  on  the  Lewis 
number. 

Sivashinsky  et  a). 134  have  also  shown  that  posi¬ 
tive  aerodynamic  stretch  has  a  stabilizing  influence 
on  the  flame.  This  effect  is  further  demonstrated  in 
Ref.  66  for  expanding  spherical  flames. 

Concentration  and  Temperature  Modifications  in 
Flame  Chemistry: 

The  study  of  fundamental  kinetics  is  sometimes 
conducted  by  employing  a  particular  flame  as  a 
"chemical  reactor"  fed  by  the  freestream  mixture, 


which  is  assumed  to  react  either  at  the  adiabatic- 
flame  temperature  or  in  a  temperature  environ¬ 
ment  prescribed  by,  say,  the  one-dimensional  lam¬ 
inar  flame.  Recognizing  that  some  of  these  flames 
suffer  stretch,  which  can  cause  modifications  in  both 
the  mixture  composition  as  well  as  the  flame  tem¬ 
perature,  failure  to  account  for  these  modifications 
can  lead  to  inaccuracies  in  the  deduced  chemical 
information.  In  fact,  the  need  for  an  accurate  de¬ 
termination  of  S[J,  as  discussed  earlier,  steins  from 
such  a  concern.  Another  example  of  the  concentra¬ 
tion  modification  is  soot  formation  in  flames, 1 ’’ 1 
as  evidenced  by  the  formation  of  soot  streak  from 
the  tip  of  the  Bunsen  flame.13,  Studies  of  polyhed¬ 
ral  flames  have  also  shown  substantial  composition 
variations  from  the  ridges  to  the  valleys. 11,1  Because 
of  the  nonlinear  nature  of  chemical  kinetics,  these 
effects  obviously  will  not  average  out  in  the  bulk 

Modelling  of  Turbulent  Flames: 

In  a  turbulent  flow  the  fluctuating  uonuniform 
local  velocities  impose  aerodynamic  stretching  on 
the  highlv-convoluted  flame  surfaces  Thus  the  re¬ 
sults  of  stretched  laminar  flames  can  be  applied  to 
certain  situations  involving  turbulent  premised 
flames.  The  minimum  requirement  for  such  an  ap¬ 
plication  is  that  the  flamelet  satisfies  the  Klimov  - 
Williams  criterion,8  139  which  basically  states  that 
the  laminar  flame  thickness  must  be  smaller  than 
the  Kolmogorov  microscale.  Peters11  has  further 
demonstrated  that  this  criterion  is  equivalent  to  Ka 
<  1.  Since  we  have  shown  that  the  extinction  Kar- 
lovitz  number,  Ka„,  is  sufficiently  smaller  than  one. 
it  seems  the  criterion  for  the  existence  ol  laniinui 
flameiets  can  be  further  tightened  to  situations  sat¬ 
isfying  Ka  <  Ka,,,. 

The  properties  of  wrinkled  laminar  flames  can  be 
related  to  turbulent  flame  propagation  In  treating 
the  wrinkled  laminar  flame  simply  as  a  surface  which 
locally  propagates  with  the  upstream  flame  velocity 
Su.  The  rates  of  flame  propagation  and  fuel  con¬ 
sumption  are  then  related  to  the  evolution  of  this 
surface  due  to  the  hydrodynamic  stretch  of  the  tur¬ 
bulent  flow  field,  with  proper  allowance  for  tin- 
variation  of  S„  with  the  local  stretch  rate  due  to 
flame  stretch. 

The  propagation  of  a  weaklv-stramed,  weakly  - 
wrinkled  flame  has  already  been  discussed  Flame 
propagation  in  medium-  and  high-intensity  turbu¬ 
lence  has  recently  been  analyzed140  141  by  using  the 
renormalization  group  method  Predicted  results 
agree  well  with  experimental  data  obtained  in  high- 
Reynolds-number  turbulent  flows,  especially  in  re¬ 
producing  the  "bending"  effect  observed  experi¬ 
mentally.11 

Another  application  of  stretch  in  turbulent  flames 
is  the  flame  dynamics  in  coherent  structures  142  The 
rolling  up  of  the  vortex  increases  the  total  surface 
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area  for  burning,  although  excessive  stretching  can 
lead  to  extinction  in  the  vortex  core. 


Stretched  Diffusion  Flames 

Compared  to  premixed  flames,  the  response  of 
diffusion  flames  to  stretch  is  less  subtle.  First  we 
note  that  even  in  a  neutrally-diffusive,  stretchless 
situation, 143  incomplete  reaction,  reactant  leakage, 
and  consequently  flame  extinction  are  inherent  fea¬ 
tures  of  a  diffusion  flame  with  finite  rate  kinetics. 
Furthermore,  since  a  diffusion  flame  does  not  have 
a  flame  speed  and  thereby  the  extent  of  flexibility 
of  a  premixed  flame  to  adjust  its  location  in  the 
presence  ol  convection,  convection  across  the  flame 
always  reduces  the  residence  time  for  reaction.  This 
leads  to  reduced  flame  thickness,  increased  reactant 
leakage,  and  consequently  increased  vulnerability 
to  extinction. 

A  stretch  factor  K  can  still  be  defined  for  a  given 
non-uniform  flow  field  and  a  wrinkled  flame  sheet, 
with  n  pointed  opposite  to  the  general  direction  of 
the  flow  velocity  v.  Stretch-induced  concentration 
and  flame  temperature  modifications  also  exist  for 
diffusionally-imbalanced  mixtures.  Preferential  dif¬ 
fusion,  while  still  quantitatively  important,  is  only 
an  additional  feature  of  the  flame  and  does  not  have 
the  profound  effect  on  the  flame  response  as  in  the 
case  of  a  premixed  flame 

A  detailed  theory  of  one-dimensional  laminar  dif¬ 
fusion  flames  was  formulated  by  Linan  for  the 
counterflow  configuration,144  with  the  assumption 
of  unity  Lewis  number.  This  theory  has  since  been 
extended  to  non-unity  Lewis  number  mixtures. 145 
Explicit  ignition/extinction  criteria  have  been  de¬ 
rived  as  function  of  the  system  Damkohler  number 
which  depends  on  the  velocity  gradient  of  the  Mow, 
allowing  a  priori  assessment  of  the  possible  exis¬ 
tence  ol  these  critical  states.  A  general  theory  of 
stretched,  wrinkled  diffusion  flames  has  not  been 
formulated.  Modeling  of  turbulent  diffusion  flames 
on  the  basis  of  laminar  diffusion  flamelets  has  been 
reviewed  by  Peters.11 


Summary  and  Suggestions 

Summarizing,  the  basic  phenomena  of  stretched 
flames  are  mostly  understsood,  at  least  qualita¬ 
tively.  Effects  of  stretch  can  be  distinguished  as  those 
caused  by  hydrodynamic  stretch,  which  does  not 
affect  the  flame  intensity  but  changes  the  flame  sur¬ 
face  area  and  thereby  the  volumetric  burning  rate, 
and  by  flame  stretch,  which  modifies  the  flame  in¬ 
tensity  and  can  also  cause  extinction.  Preferential 
diffusion  interacts  strongly  with  stretch  and  can 
qualitatively  affect  the  flame  response  depending  on 


whether  the  effective  mixture  Lewis  number  is 
greater  or  smaller  than  unity,  and  whether  the 
stretch  is  positive  or  negative.  The  extent  with  which 
the  flame  can  freely  adjust  its  location  in  response 
to  stretch,  in  order  to  achieve  complete  reaction 
is  also  an  important  factor  especially  for  extinction 
considerations,  ft  is  further  clarified  that  llame  sta¬ 
bilization  is  an  issue  of  dynamic  stability,  it  is  nei¬ 
ther  a  stretehed-flame  phenomenon  nor  an  extinc¬ 
tion  phenomenon. 

More  research,  however,  is  still  needed  for  com¬ 
plete  understanding  and  further  application  of  the 
stretched  flame  phenomena  Some  of  the  various 
issues  ol  current  interest  and  importance  are 

1)  A  nonlinear  analysis  ol  the  gcncial  stretched 
premixed  flame,  allowing  lor  flame  extinction 
is  needed.  The  formulation  can  lie  facilitated  In 
the  observation  that  Kar,  is  sufficiently  smallei 
than  one.  The  analysis  should  provide  turthei 
insight  into  the  insensitivity  ol  the  properties 
of  the  adiabatic,  unrestrained,  diftusionallx -neu¬ 
tral  flames  to  stretch,  especially  foi  large  stretch 
rates. 

2)  The  anticipation  that  energy  conservation  can 
be  satisfied  for  stagnation  and  Bunsen  Haines 
by  considering  tangential  transport  needs  to  be 
demonstrated. 

3)  The  role  played  by  the  preferential  diffusion  ol 
reaction  intermediates,  especially  such  highly  - 
mobile  and  chemicallv-crucial  species  is  the 
hydrogen  radical,  on  the  dynamics  ol  stretched 
flames  requires  more  detailed  studs  Asymp¬ 
totic  analysis  using  reduced  kinetic  schemes" 
appears  to  be  especially  suitable  for  such  an  in¬ 
vestigation. 

4)  Numerical  modelings  of  stretched  flames,  es¬ 
pecially  those  in  the  counterflow  configuration 
have  employed  different  flow  fields  While  tins 
is  acceptable  for  qualitative  studies  a  tunic  n 
ulistic  flow  field  is  needed  when  pertoimmg 
quantitative  comparison  with  expeiiiuent.il 
data  "  The  concern  is  especially  serious  i!  tin 
chemical  kinetic  mechanism  used  is  also  not 
certain  and  possibly  requires  "adjustment 

5)  Ex|)crimental  stretched  flame  results  are  fre¬ 
quently  reported  in  terms  of  the  global  stretch 
rate.  This  can  be  quite  different  from  the  local 
stretch  rate  determined  by  using,  say.  LDV 
Caution  is  needed  when  comparing  these  data 
with  the  theoretical  results  obtained  by  using 
local  stretch  rates 

6)  Explicit  theoretical  formulas  have  already  been 
derived  for  yveaklv-stretched  flames,  at  least  foi 
one-step  reactions.  Experimental  results,  how¬ 
ever.  have  been  scarce  A  systematic  expen- 
mental  effort  should  be  mounted  to  obtain  data 
of  high  lidelitv  lor  stretched  flames  with  posi¬ 
tive/negative  stretches  and  different  preferen- 
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tial  diffusion  effects,  and  to  correlate  them  ac¬ 
cording  to  the  theoretical  formulas. 

71  More  effort  is  needed  to  integrate  the  theoret¬ 
ical  results  on  stretched  flames  in  the  modeling 
of  turbulent  flames.  Caution  is  needed  to  en¬ 
sure  overall  energy  conservation  lor  diffusion- 
alb -imbalanced  mixtures.  Flame  interaction, 
discussed  next,  also  needs  to  be  accounted  lor 

S)  Most  ol  the  studies  on  stretched  flames  have 
employed  a  solitary  flame,  which  is  either  pre¬ 
mixed  or  nonpremixed.  However,  in  the  inte¬ 
rior  of  a  turbulent  flame,  the  individual  fla 
melets  interact  with  each  other.  Partial 
premixing  also  occurs,  resulting  in  partially 
premised  flames.  Some  research  has  been  con¬ 
ducted  on  these  issues, "!i  but  more  is 
needed 

9i  The  existence  of  flammability  limits  needs  to  be 
settled.  It  they  indeed  exist,  then  a  general 
theory  ol  flammability  limits  is  needed  Fur¬ 
thermore,  since  a  sub-limit  mixture  can  be  ren¬ 
dered  flammable,  through  either  flame  stretch 
and/or  flame  interaction,  the  absolute  signifi¬ 
cance  of  flammability  limits  also  lias  to  be  as¬ 
sessed. 

ltb  A  general  theory  of  stretched  diffusion  flames 
is  needed 
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Nomenclature 

a  strain  rate  ol  stagnation  flow 

A  area  ol  flame  surface 

II  (requeues  factor 

i',,  specific  heat 

Da  Dainkohler  iminbcr 

1)  dittusivitv 

k  reaction  rate  constant 


K  stretch  rate  isec  1 1 

Ka  karlos  itg  number 

E  loss  parameter 

lx  Lewis  number 

in  mass  burning  rate  per  unit  area 

11  overall  reaction  order 

n  unit  normal  vector 

p  pressure 

1/  beat  ol  combustion 

li/  radial  dimension  ol  flame 

t  time 

T  T  temperature  and  stoicliioinetricalls -weighted 

T 

T„  activation  temperature 

T „,/  adiabatic  flame  temperature 

u.u  flow  velocity 

V  flame  velocity 

ice  reaction  rate  ol  kth  reaction 

x.y  spatial  coordinate 

Y.Y  mass  fraction  and  stoichiometi  icalls - 
ssciglitcd  Y 

Greeks 

0  coupling  Itinction 

5  thickness  ol  various  flame  /ones 

\  thermal  conductivity 

<t>  equivalence  ratio 

p  density 

<r  gcometrs  factor  lor  stagnation  floss 

H  angle  ol  Kunseii  cone 

.Nubvrripf.s 
If  blossoll 

D.ll.K  Transport  hydrodynamic.  rcactmn  /ones 

M  l  species  and  thermal  diffusion  /ones 

c.v  extinction 

/  (lame 

i.j  index  for  species 

11. o  imhurut  and  hiirnt  states  ol  mixture 

s  surface,  especially  flames 

/-  ambience 

Superscripts 

o  the  adiabatic,  unstretehod  planar  flame 
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FLAMES  IN  STAGNATION  FLOW 
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Detailed  experimental  measurements  of  the  scalar  and  velocity  statistics  of  premixed  meth¬ 
ane/air  flames  stabilized  by  a  stagnation  plate  are  reported.  Conditioned  and  unconditioned 
velocity  of  two  components  and  the  reaction  progress  variable  are  measured  by  using  a  two- 
component  laser  Doppler  velocimetry  technique  and  Mie  scattering  technique,  respectively. 
Experimental  conditions  cover  equivalence  ratios  of  0.9  and  1.0.  incidence  turbulence  in¬ 
tensities  of  0.3  to  0.45  m/s,  and  global  stretch  rates  of  100  to  150  sec"'.  The  experimental 
results  are  analyzed  in  the  context  of  the  Bray-Moss-Libbv  flamelet  model  of  these  flames. 
The  results  indicate  that  there  is  no  turbulence  production  within  the  turbulent  flame  brush 
and  the  second  and  third  order  turbulent  transport  terms  are  reduced  to  functions  of  the 
difference  between  the  conditioned  mean  velocity.  The  result  of  normalization  of  these  rel¬ 
ative  veleocities  by  the  respective  velocity  increase  across  laminar  flame  suggests  that  the 
mean  unconditioned  velocity  profiles  are  self-similar. 


Introduction 

In  a  recent  paper1  we  introduced  the  use  of  the 
turbulent  stagnation  flow  configuration  for  funda¬ 
mental  studies  of  premixed  turbulent  flames.  The 
turbulent  flame  propagates  in  the  divergent  stag¬ 
nation  flow  region  of  a  turbulent  fuel/air  jet  im¬ 
pinging  on  a  plate.  Depending  on  flow  and  mixture 
conditions  the  flame  brush  is  stabilized  at  various 
distances  away  from  the  plate  surface  where  the 
mean  velocity  balances  the  turbulent  flame  speed. 
Consequently,  flame  stabilization  is  achieved  with¬ 
out  involving  flow  recirculation  as  in  the  conven¬ 
tional  V-flame  and  conical  flame  configurations. 
Moreover,  at  the  centerline,  the  turbulent  flame 
brush  is  planar  and  locally  normal  to  the  approach 
flow  so  that  the  turbulent  flame  speed  can  be  de¬ 
duced  unambiguously  from  the  flow  velocity  enter¬ 
ing  the  reaction  zone.  Therefore,  this  configuration 
provides  a  convenient  and  more  accurate  means  to 
deduce  the  turbulent  flame  speed. 1 

The  most  significant  theoretical  implication  of  this 
configuration,  however,  is  that  it  is  the  closest  ex¬ 
perimental  approximation  to  the  idealized  one-di- 
mcnsional  infinite  planar  flame  assumed  in  many 
theoretical  turbulent  combustion  models.23  Bray, 
Champion  and  Libby1  (BCL)  recognize  that  a  the¬ 


oretical  analysis  for  the  infinite  planar  flame  is  lim¬ 
ited  because  the  idealized  situation  cannot  be  pro¬ 
duced  in  the  laboratory.  As  a  result,  they  have  relied 
on  a  comparison  of  their  theoretical  predictions  with 
experimental  measurements  made  in  2D  oblique 
flames.  Although  these  comparisons  have  been  suc¬ 
cessful  in  validating  many  features  of  the  model,  it 
is  however  difficult  to  ascertain  whether  the  dis¬ 
crepancies  between  predictions  and  measurements 
are  caused  by  variations  of  the  flame  properties  not 
accounted  for  in  the  theory,  or  by  the  actual  de¬ 
ficiencies  in  the  theory.  Hence,  encouraged  in  part 
hv  our  experimental  work,  they  have  recently  ini¬ 
tiated  the  modeling  of  the  stagnation  flow  stabilized 
flame  so  that  a  more  direct  and  compatible  com¬ 
parison  can  be  made  between  theoretical  and  ex¬ 
perimental  studies.  Their  first  analysis'1  is  for  the 
case  of  flames  stabilized  in  two  opposed  jets  and 
closure  is  achieved  by  the  use  of  conventional  gra¬ 
dient  transport.  Subsequent  studies^  will  model  the 
turbulent  transport  using  the  Bray-Moss-Libby 
(BML)  flamelet  concept  which  assumes  infinite  rate 
chemistry  and  zero  flamelet  thickness." 

The  stagnation  flame  configuration  is  also  signif¬ 
icant  for  other  numerical  models  of  turbulent  com¬ 
bustion.  such  as  the  joint  probability  density  func¬ 
tion  (jpdf)  method  developed  by  Pope’'’  and  a 
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deterministic  model  lor  simulating  the  interaction 
of  the  flamelet  with  incident  turbulence  using  the 
vortex  dynamics  technique.7  At  present,  the  jpdf 
model  is  capable  of  predicting  the  turbulent  flame 
speed  and  other  properties  similar  to  those  calcu¬ 
lated  by  the  BML  model.  Vortex  dynamic  simula¬ 
tions  of  the  stagnation  flames  compute  the  dynamic 
movement  of  the  wrinkled  flamelet  and  also  predict 
the  conditioned  velocity  statistics.  This  configura¬ 
tion  is,  therefore,  important  for  the  further  devel¬ 
opment  of  fundamental  experimental  and  theoreti¬ 
cal  turbulent  combustion  research. 

Since  very  little  is  known  about  the  behavior  of 
the  stagnation  flow  stabilized  flames,  it  is  necessary 
to  improve  our  understanding  of  their  fundamental 
properties  through  detailed  experimental  measure¬ 
ment  of  the  scalar  and  velocity  statistics.  These  ex¬ 
perimental  data  will  be  useful  in  the  development 
of  theoretical  models  for  this  configuration.  In  re¬ 
sponse  to  such  a  need,  we  report  herein  results  from 
an  investigation  of  the  scalar  and  velocity  statistics 
which  provide  further  characterization  of  the  flame 
properties.  Measurements  were  made  of  the  con¬ 
ditioned  and  unconditioned  velocity  of  two  com¬ 
ponents  and  the  reaction  progress  variable  c\  a 
nondimensionalized  scalar.  Features  of  the  mean  and 
rms  velocity  profiles  significant  to  modeling  arc  dis¬ 
cussed  In  addition,  implications  tor  the  closure  re¬ 
quirements  of  the  third  order  turbulence  transport 
are  also  considered. 


Experimental  Methodology 

A  uniform  axisymmetrie  flow  of  premixed  meth¬ 
ane/air  is  provided  by  a  50  mm  diameter  nozzle. 
This  jet  of  reactants  is  shielded  by  a  co-flow  of  air 
with  outer  diameter  of  100  mm.  For  this  study,  the 
stagnation  plate  is  placed  at  75  mm  above  the  exit, 
and  the  exit  flow  velocity  is  held  constant  at  5.0 
m/s.  The  incident  turbulence  is  generated  by  plac¬ 
ing  at  50  mm  upstream  of  the  exit  a  square  grid  of 
1.0  mm  elements  with  5.0  nun  openings  or  a  per¬ 
forated  plate  with  3.2  mm  diameter  circular  holes 
•  hich  gives  a  40%  blockage  ratio.  Mixture  and  flow 


conditions  of  the  four  flames  are  listed  m  Table  1 
The  turbulence  Reynolds  numbers  and  Damkohlcr 
numbers  determined  at  the  nozzle  exit  indicate  that 
these  flame  fall  within  the  wrinkled  laminar  flame¬ 
let  regime  (Table  1). 

Flow  velocities  are  measured  by  a  four-beam  two- 
color  laser  Doppler  velocimetry  (LDV)  system  which 
measures  the  velocity  of  the  txvo  components  using 
a  10  (isec  co-validation  criterion  Details  of  the  LDY 
system  and  the  computer  controlled  data  acquisi¬ 
tion  system  can  be  found  in  Ref.  1.  Seeding  lor 
unconditioned  velocity  measurement  is  pros  ided  In 
refractory  aluminum  oxide  partieles  ol  0.3  pm  in¬ 
troduced  into  the  jet  by  a  cyclone  seed  generator 
Conditioned  velocity  measurement  in  the  reactant 
side  is  made  possible  In  using  a  silicon  oil  aerosol 
as  the  LDV  seed.  The  technique  is  based  on  the 
principle  that  the  silicon  oil  droplets  evaporate  and 
burn  through  the  flainelets.  A  blast  atomizer  is  used 
to  generate  the  aerosol.  At  each  measurement  po¬ 
sition,  8192  pairs  of  validated  velocity  data  are  n 
corded.  To  deduce  the  conditioned  veloeitv  statis¬ 
tics  in  the  products,  a  procedure  which  involves 
deeonvoluting  the  velocity  pdf’s  as  described  in 
Cheng  and  Shepherd’’  is  used. 

The  silicone  oil  seeding  technique  is  also  used 
for  measuring  the  conserved  scalar,  c.  by  monitor¬ 
ing  the  Mie  scattering  intensity  using  a  photomul¬ 
tiplier  and  lens  assembly.  Within  the  flame  brush, 
the  Mie  scattering  signal  resembles  a  random  tel¬ 
egraph  signal  with  sharp  transition  between  the  sig¬ 
nal  levels  representing  the  reactant  and  product 
states.  By  specifying  a  threshold  criterion,  the  local 
value  of  the  conserved  scalar  is  determined  by  ev  al¬ 
uating  the  percentage  of  the  time  spent  in  the 
products.  Using  the  zero  thickness  flamelet  approx¬ 
imation.  the  conserved  scalar  so  determined  is 
identical  to  the  normalized  density  ip  -  p,,i/tp,  - 
p,,),  where  p  is  the  mean  density  and  subscripts  r 
and  p  denote  reactants  and  products.  The  Mie  scat¬ 
tering  signal  can  also  be  analyzed  to  obtain  other 
time  dependent  scalar  parameters  such  as  the  scalar 
spectra,  flame  crossing  frequency'1  and  pdf  s  of  the 
flame  passage  times.10  The  turbulent  flame  brush 
thickness  is  obtained  from  the  distance  between  the 
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two  points  on  the  axial  scale  intersected  by  the  tan¬ 
gent  to  the  c  profile  at  c  =  0.5. 

Results  and  Discussions 

As  reported  in  our  previous  study,  typical  flow 
velocities  in  the  llamt  region  are  of  the  same  order 
of  magnitude  as  the  turbulent  flame  speed.  These 
velocities  are  lower  than  those  of  the  other  axisym- 
metric  or  plane-symmetric  flame  configurations  in 
which  the  velocity  of  the  flow  entering  the  flame 
region  is  about  the  same  as  that  of  the  freestream. 
Therefore,  it  becomes  necessary  to  determine 
whether  or  not  the  models  developed  for  high  Rey¬ 
nolds  numbers  are  still  valid  for  the  stagnation 
flames.  The  two  criteria  for  the  models  are  that  Re/, 
the  turbulent  Reynolds  number  based  on  the  in¬ 
tegral  length  scale  and  velocity  fluctuations,  and  Da, 
the  large  eddy  Damkohler  number,  are  both  much 
large  than  unity.  Although  these  conditions  are  sat¬ 
isfied  at  the  exit,  it  is  important  to  determine  if  the 
turbulence  scales  change  close  to  the  plate. 

Shown  in  Fig.  1  are  the  profiles  of  the  time  scales 
ol  the  non-reacting  stagnation  point  flow,  ob¬ 
tained  from  auto-correlations  and  the  associated 
length  scale,  lx ,  deduced  by  Taylor’s  hypothesis  for 
ease  #1.  The  length  scale  is  relatively  constant  ex¬ 
cept  at  x  <  10  mm  where  a  decrease  is  found.  This 
decrease  seems  to  be  an  indication  of  the  break¬ 
down  of  the  Taylor's  hypothesis  in  the  low  velocity 
region.  Similar  features  are  also  found  in  the  tur¬ 
bulence  generated  by  the  perforated  plate.  Al¬ 
though  it  is  not  possible  to  deduce  the  integral 
length  scale  for  the  transverse  component,  the  time 
scales  show  that  the  grid  generated  turbulence  is 
isotropic  whereas  the  turbulence  generated  bv  the 
perforated  plate  is  slightly  anisotropic  with  smaller 
transverse  integral  time  scales.  The  result  that  both 


the  length  scale  and  the  velocity  fluctuations  arc 
constant1  demonstrates  that  the  flow  criteria  tor 
current  turbulent  combustion  models  are  satisfied 
even  close  to  the  stagnation  flame  region. 

The  mean  centerline  axial  profiles  of  the  condi¬ 
tioned  and  unconditioned  velocities  for  ease  #  1  are 
shown  in  Fig.  2.  Also  shown  is  the  profile  of  the 
reaction  progress  variable.  It  shows  that  the  reac¬ 
tant  boundary  of  the  flame  zone,  as  indicated  bv 
the  position  where  c  deviates  from  zero,  occurs  at 
the  minimum  of  the  unconditioned  velocity  profile 
This  provides  an  experimental  v  alidation  of  our  pre¬ 
vious  method  of  determining  S-;-  where  the  mini¬ 
mum  velocity  was  equated  to  the  turbulent  flame 
speed. 1 

The  turbulent  flame  brush  thickness.  5/,  de¬ 
duced  from  the  c  profiles  is  listed  in  Table  1  Also 
listed  is  the  corresponding  laminar  flame  thickness, 
bL,  based  on  the  definition  of  Zeldovich. 11  The  ra¬ 
tio  between  the  turbulent  and  laminar  llame  thick¬ 
nesses  is  a  significant  parameter  in  the  BML,  model 
which  assumes  infinite  rate  chemistry  .  Studies  in 
turbulent  Y-flames  and  Bunsen  flames  have  shown 
that  at  positions  close  to  the  stabilization  region,  the 
zero  thickness  llamelet  assumption  is  not  valid  lie- 
cause  the  density  pdf  s  show  appreciable  contribu¬ 
tions  from  the  intermediate  burning  states.  These 
conditions  are  satisfied  only  at  positions  farther  away 
because  the  flame  brush  thickness  increases  almost 
linearly  w'ith  downstream  distance.  In  contrast,  since 
the  stagnation  flames  do  not  have  recirculating  sta¬ 
bilization  regions,  the  flame  brush  thickness  re¬ 
mains  relatively  constant  along  the  flame  front. 
Therefore,  the  process  of  choosing  a  location  suit¬ 
able  for  comparison  with  theory  is  not  encountered. 
Moreover,  the  smallest  turhulcnt/laminar  flame 
thickness  ratio  of  30  (case  #2)  indicates  that  the  zero 


Distance  trom  plate  (mm) 


Distance  from  plate  (mm) 


FlC.  1  The  time  scales  of  the  non-reacting  stag¬ 
nation  flow,  obtained  from  auto-correlation  and 
the  associated  length  scale,  lx,  obtained  by  Taylor's 
hypothesis. 


Fie.  2.  Mean  centerline  axial  profile  of  the  con¬ 
ditioned  and  unconditioned  velocities  for  case  #  1 
(<t>  =  1  and  u'  =  40  cm/s).  Inset  profile  of  the 
reaction  progress  variable. 
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thickness  flamelet  approximation  is  valid  for  tur¬ 
bulent  stagnation  flames  presently  considered. 

Another  important  aspect  of  the  scalar  field  is  the 
heat  transfer  to  the  stagnation  plate  Since  the  plate 
is  made  of  brass,  the  turbulent  flame  may  be  influ¬ 
enced  by  heat  loss  from  the  llame  /one  To  inves¬ 
tigate  the  possible  effect  of  heat  loss,  thermocouple 
measurements  were  made  of  the  plate  temperature 
and  the  thermal  boundary  layer  over  the  plate.  The 
results  for  case  #4  show  that,  although  the  plate 
temperature  is  only  1500  K,  the  thermal  layer 
thickness  is  about  3.0  mm,  which  is  smaller  than 
the  standoff  distance  of  the  flame  brush.  This  in¬ 
dicates  that  the  wall  effect  is  confined  within  the 
product  region  and  is  not  expected  to  have  any  di¬ 
rect  influence  on  flame  propagation 

The  features  of  the  unconditioned  axial  velocity 
profile  (Fig.  2)  are  typical  of  those  measured  in  the 
stagnation  flames  as  reported  earlier  1  The  uncon¬ 
ditioned  velocity  profile  shows  three  distinct  re¬ 
gions.  The  first  region  extends  from  the  exit  of  the 
nozzle  to  the  upstream  boundary  of  the  turbulent 
flame  brush  It  is  characterized  by  a  steady  de¬ 
crease  in  the  axial  velocity  due  to  the  divergence 
of  the  cold,  unreacted  flow.  The  global  stretch  rate 
imposed  on  the  flow  is  given  by  the  mean  velocity 
gradient  in  this  region.  The  second  region  is  char¬ 
acterized  by  a  rapid  increase  in  velocity,  indicating 
flow  acceleration  due  to  combustion  heat  release. 
The  third  region  is  where  the  product  stream  is  again 
decelerated  by  the  presence  of  negative  pressure 
gradient  and  the  flow  is  eventually  brought  to  rest 
at  the  stagnation  plate. 

The  upstream  stretch  rate,  a.  determined  from 
these  profiles  (Table  1)  shows  that  the  use  of  the 
grid  tends  to  produce  a  slightly  higher  upstream 
stretch  rate.  This  global  stretch  rate  is  of  the  same 
magnitude  as  the  local  stretch  rate  on  the  flamelctx, 
which  can  be  estimated  by  u'/l ,  However,  studies 
on  laminar  flames  “  have  shown  that  such  small 
stretch  rates  have  minimal  effect  on  the  local  burn¬ 
ing  rate. 

Within  the  flame  zone,  the  unconditioned  mean 
velocity  increases  to  form  a  "hump  at  the  product 
boundary  as  shown  in  Fig.  2.  As  discussed  previ¬ 
ously.1  tins  increase  is  due  to  intermittent  mea¬ 
surement  of  the  conditioned  velocities,  U,  and  l’,,. 
’file  relative  velocity  between  the  conditioned  mean 
velocities.  All  —  U(,  -  l',.  which  is  referred  to  as 
the  slip  velocity  in  the  BML  model,  is  a  major  con¬ 
tributor  to  the  unconditioned  fluctuations  and  the 
second,  and  third  order  turbulent  transport  terms. 
The  cause  of  the  so-called  counter  gradient  scalar 
transport  is  a  direct  result  of  the  fact  that  A  U  is 
positive  in  most  flames.  Shown  in  Fig.  3  are  the 
four  A l'  profiles  plotted  against  c  obtained  in  this 
study  The  data  are  normalized  by  the  increase  in 
velocity  An  -  S(  (t  -  1)  across  a  planar  laminar 
flame  where  t  is  the  reactant /product  density  ratio 
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it  can  be  seen  the  Al  s  are  almost  constant  and 
their  values  are  about  half  ol  those  of  An  Previous 
measurements  of  this  ratio  in  V -flames  and  conical 
Bunsen  flames  have  shown  that  they  are  indepen¬ 
dent  of  r  though  they  are  generally  lower  than  unity  . 
As  discussed  by  Cheng  et.  al.,u  flow  divergence  in 
the  products  and  random  flow  deflections  due  to 
the  wrinkled  flame  could  lower  this  ratio.  How¬ 
ever,  the  observation  of  constant  Al'/An  ratio  for 
different  conditions  in  the  stagnation  flames  is  quite 
encouraging  because  it  implies  that  the  conditioned 
and  unconditioned  velocity  profiles  may  be  self- 
similar.  This  feature  ol  the  stagnation  flame  profiles 
is  anticipated  to  facilitate  the  modeling  of  these 
flames.  At  present,  similarity  variables  for  the  dis¬ 
tance  from  the  plate  surface  is  proposed  in  IK. "I,  1 
The  proper  means  to  nondimonsionalize  the  veloc¬ 
ity,  however,  requires  further  studies. 

The  conditioned  and  unconditioned  mis  fluctua¬ 
tion  profiles  for  the  axial  and  transverse  compo¬ 
nents  are  shown  in  Fig.  4.  The  peak  shown  in  the 
unconditioned  u"  profiles  (Fig  4ia'i  represents  con¬ 
tributions  from  AU.  The  more  significant  feature, 
however,  is  that  the  two  axial  conditioned  nils  fluc¬ 
tuations  are  almost  equal.  This  suggests  that  then- 
is  no  flame  generated  turbulence  in  the  axial  di¬ 
rection.  The  lack  of  a  peak  in  the  r’  profiles  is  due 
to  the  fact  that  the  transverse  component  ol  tin- 
relative  velocity  is  zero.  The  ditleroncc  in  the  con¬ 
ditioner.  fluctuations,  r-J,  —  r,\  is  slightly  higher  at 
positions  close  to  the  product  boundary  than  those 
shown  for  the  axial  component.  Recent  analysis  ot 
turbulence  production  across  turbulent  flamelets  have 
shown  that  the  increase  in  turbulence  is  associated 
with  vorticity  production.'  Vorticitv  production  is 
expressed  bv  the  baroelinic  term  which  is  the  cross 
product  between  the  pressure  gradient  normal  to 
the  flamelet  and  the  local  flame  curvature  due  to 
stretch  In  the  stagnation  flame  the  adverse  mean 
axial  pressure  gradient  max.  In  reducing  the  mag- 
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Fit;.  4.  Conditioned  and  unconditioned  rins  ve¬ 
locity  fluctuations  for  (a)  axial,  (b)  transverse  com¬ 
ponents. 

nitude  of  the  baroclinic  tern".,  suppress  production 
of  vorticity.  It  is,  however,  important  to  emphasize 
the  tentative  nature  of  this  interpretation,  pending 
further  experimental  evidence  and  theoretical  anal¬ 
ysis. 

The  lack  of  significant  flame  generated  turbu¬ 
lence,  in  particular  in  the  axial  fluctuations,  also  has 
a  significant  impact  on  the  modeling  of  the  third 
order  turbulent  transport  terms.  For  example,  the 
turbulent  diffusion  of  scalar  is  modeled*  as 

p u"2c"  .  _  „  , 

- -  =  c(l  -  c){(Up  -  Ur)V 

P 

+  (1  -  c)2]  -  ur'2  +  u/,2} 

the  turbulence  flux  is 

pu”3  .  .  —  .  .  — 

-  =  (1  -  c){ 3(Ur  -  I’lu,:2  +  (Ur  -  U)!  +  ur'3} 

p 

+  e{3(C„  -  U)^T2  +  (C'„  -  U)3  +  V3} 

where 

0  =  cUp  +  (1  -  c)Ur 

These  sub-models  consist  of  contributions  from 
AC.  (he  difference  in  the  conditioned  turbulence 


intensities,  and  the  third  order  terms.  However,  as 
the  conditioned  velocity  pdf s  are  essentially  Gaus¬ 
sian,  these  third  order  terms  may  be  neglected  s 
In  previous  theoretical  studies,  Bray  et  al  1  used 
an  empirical  constant  to  relate  the  increase  in  the 
conditioned  velocity  fluctuations.  To  study  the  rel¬ 
ative  contributions  from  the  two  terms,  the  ratios 
(Up2  —  tir2)/(AU)2  and  (Up2  —  t>'2)/(AU  )2  are  plotted 
in  Fig.  5.  As  can  be  seen  the  magnitudes  of  these 
ratios  are  very  small,  although  (v'/  -  c'2)/(AU  '2 
tend  to  be  slightly  higher  than  (« 'p~  -  u,'2)/(AL-r 
for  cases  #3  and  #4.  The  lack  of  significant  tur¬ 
bulence  production  by  the  flame  suggests  that  all 
the  third  order  transport  can  be  modeled  simply  by 
considering  the  contributions  from  the  relative  ve¬ 
locity  AU.  This  again  simplifies  the  modeling 
scheme. 


Concluding  Remarks 

Reported  in  this  paper  is  a  detailed  experimental 
investigation  of  the  characteristics  of  the  velocity  and 
scalar  fields  in  stagnation  flow  stabilized  premixed 
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turbulent  llames.  This  relatively  new  turbulent  flame 
configuration  is  considered  to  be  well  suited  for  di¬ 
rect  comparisons  between  experimental  and  theo¬ 
retical  predictions.  The  measurements  include  con¬ 
ditioned  and  unconditioned  velocity  statistics  of  two 
components,  and  the  profiles  of  the  mean  reaction 
progress  variable  which  is  identical  to  the  normal¬ 
ized  density  in  the  limit  of  zero  thickness  flainelet 
assumption.  The  turbulence  Reynolds  numbers  and 
Damkohler  numbers  for  the  four  methane/air  flames 
indicate  that  these  flames  are  within  the  wrinkled 
flamelet  regime.  The  experimental  results  arc  ana¬ 
lyzed  to  provide  insight  which  will  assist  in  the  de¬ 
velopment  of  the  Bray-Moss-Libbv  flamelet  model 
for  these  flames. 

For  all  the  flow  and  mixture  conditions  consid¬ 
ered  here,  the  relative  velocity  AU  between  the 
reactants  and  products,  i.e.  the  difference  in  the 
conditioned  velocities  is  constant  across  the  flame 
brush.  Moreover,  the  relative  velocities,  when  nor¬ 
malized  by  the  velocity  increase,  A u  —  SL  (t  -  1), 
under  laminar  condition,  have  the  same  value.  This 
suggests  that  the  mean  unconditioned  velocity  pro¬ 
files  across  the  flame  region  are  self-similar.  These 
characteristics,  il  confirmed  by  more  extensive 
studies,  would  greatly  simplify  the  modeling  of  the 
complex  turbulent  flame  flow  field 

Within  the  flame  brush,  the  conditioned  turbu¬ 
lence  intensities  for  the  axial  and  transverse  veloc¬ 
ity  components  are  about  equal  and  constant.  That 
is,  the  stagnation  flames  do  not  create  additional 
turbulence.  This  is  quite  different  form  the  results 
obtained  in  oblique  flames  where  flame  generated 
turbulence  occurs.  The  most  important  implication 
to  modeling  is  that  all  the  second  and  third  order 
turbulence  transport  terms  are  reduced  to  functions 
of  the  relative  velocity  AU.  This  again  implies  self- 
similarity  of  the  turbulence  and  transport  proper¬ 
ties 
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COMMENTS 


A  (ilmniciu.  Mil'.  (’SA.  Wli\  dill  von  expect 
turlmlciK  e  generation  in  tins  stagnation-point  flame:' 
Are  not  tin*  density  and  pressure  gradients  almost 
aligned  in  this  case  especially  when  you  consider 
that  the  thermal  expansion  in  the  radial  direction 
is  almost  enough  to  balance  the  hdrodynainic  pres¬ 


jXulhor  s  liffih/  For  othei  premised  turbulent  « 
flames  ic  g  \ -flames  oi  Bunsen  ll.uues  .  tin  bulein  i 
production  is  shown  in  increases  in  the  eonditionei 
rms  fliictiialions  in  the  product  /ones  The  lack  n 
increase  observed  here  indicates  tii.it  the  centerlm 
ol  the  stagnation  flow  stabilized  flames  is  an  excep 
lion  Whether  or  not  this  is  due  to  the  signiluan 
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streamtube  divergence  as  suggested  would  require 
the  support  of  data  obtained  away  from  the  cen¬ 
terline.  In  particular,  measurements  made  along 
Ilowlines  (see  Cheng,  et  aL,  22nd  Combustion 
Symposium)  would  be  most  useful. 


S  .V  B  Murphy,  Purdue  Vniv..  USA.  A  turbu¬ 
lent  flow  is  governed  by  rather  restricted  similarity 
considerations  at  a  stagnation  point.  What  is  the  basis 
for  the  choice  of  curvature  in  the  experiments?  Are 
there  any  data  to  indicate  an  acceptable  curvature 
for  burning  rate  determination? 

Author’s  Reply.  As  mentioned  in  our  previous  pa¬ 
per  (Cho  et  al.  21st  Symposium,  p.  1493)  the  cri¬ 
terion  for  the  plate  curvature  was  to  enable  the  LDA 
probe  to  measure  flow  velocities  within  3.0  mm  from 
the  plate  surface. 

At  the  centerline  of  the  stagnation  flow,  the  only 
effect  of  the  plate  curvature  may  be  the  overall 
stretch  rate  This  overall  stretch  rate  can  be  de¬ 
duced  from  the  velocity  profiles.  Therefore,  the 
question  regarding  the  “acceptable  range  for  the 
curvature  is  related  to  the  issues  concerning  the  ef¬ 
fects  of  overall  stretch  on  flame  propagation.  Per¬ 
haps  the  effects  of  stretch  can  be  studied  by  vary¬ 


ing  the  plate  curvature  in  contrast  to  varying  the 
separation  distance  between  the  plate  and  the  noz¬ 
zle  exit 

• 

T  A.  Baritaud.  Institut  Fraticais  du  Pet  role. 
France.  You  see  that  u'  remains  unchanged  when 
crossing  the  flame  front,  when  U  is  increased  due 
to  the  expansion  of  the  gas.  Hence,  p„u„J  is  almost 
instantaneously  transformed  in  a  quantity  much 

smaller.  Where  is  this  turbulent  kinetic  energ. 
going:  merest*  of  the  dissipation,  effect  of  the  pres¬ 
sure  gradient,  transfer  to  the  other  velocity  com¬ 
ponents  .  .  .? 

Author’s  Reply.  A<  shown  in  Fig.  4,  some  of  the 
turbulence  kinetic  energy  is  being  transferred  to  the 
v  component.  Nevertheless,  the  total  kinetic  en¬ 
ergy  after  weighing  by  the  local  density  is  indeed 
smaller.  This  is  true  for  most  premixed  turbulent 
flames.  The  present  situation  may  be  an  extreme 
case  where  the  effects  of  streamtube  divergence  are 
also  significant.  As  mentioned  in  the  answer  to  Prof. 
Ghoniem.  conclusions  regarding  the  turbulence 
production  in  these  flame  can  onh  be  drawn  In 
comparing  measurements  made  along  oil -center 
flowlines. 
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An  experimental  investigation  has  been  conducted  on  the  influence  of  the  mobility  of  inert 
additives  on  soot  formation  in  propane  and  ethylene  counterflow  diffusion  flames.  Inerts  used 
were  helium,  neon,  argon,  or  krypton,  and  the  results  show  that  while  the  mobility  of  the 
inert  has  practically  no  effect  when  a  small  amount  of  inert  is  added  to  the  oxidizer  side, 
the  influence  is  significant  when  added  to  the  fuel  side  in  that  krypton,  being  the  least 
mobile  inert,  yields  the  greatest  soot  loading  while  helium,  being  the  most  mobile,  yields 
the  least.  By  relating  the  spatially-resolved  soot  volume  fractions  to  the  corresponding  profiles 
of  temperature,  velocity  and  species  concentrations,  it  is  demonstrated  that  this  influence  on 
soot  loading  is  likely  caused  by  concentration  modifications  of  the  fuel  and  the  soot  precursors 
due  to  the  different  mobilities  of  the  inert  additives. 


Introduction 

Flame  phenomena  are  strongly  influenced  by  the 
diffusive  processes  of  the  reactants.  Since  the  dif- 
fusivities  of  the  various  species  in  a  gaseous  mixture 
are  in  general  not  equal,  it  is  reasonable  to  expect 
that  the  local  reactant  concentration  will  be  depen¬ 
dent  on  the  extent  with  which  the  diffusivities  dif¬ 
fer.  This  concentration  modification  through  pref¬ 
erential  diffusion  can  in  turn  affect  the  flame  behavior 
from  that  anticipated  based  on  the  freestream  con¬ 
centrations.  Indeed,  there  exist  many  flame  phe¬ 
nomena  which  can  only  be  satisfactorily  explained 
on  the  basis  of  preferential  diffusion.  For  premixed 
(lames,  we  cite  the  earlier  works  of  Refs.  1  and  2 
and  the  recent  advances  as  summarized  in  Ref.  3. 
The  amount  of  research  on  nonpremixed  flames, 
however,  has  been  less;  some  recent  contributions 
include  analyses4'5  and  experiments6  on  flame  ex¬ 
tinction. 

The  objective  of  the  present  work  is  to  study  the 
effect  of  preferential  diffusion  and  thereby  concen¬ 
tration  modification  on  an  important  nonpremixed 
flame  phenomenon  in  which  diffusive  transport  and 
chemical  kinetic  processes  are  inherently  coupled. 


*  Present  address:  Department  of  Mechanical  and 
Aerospace  Engineering,  Princeton  University, 
Princeton,  NJ  08544 


namely  soot  formation  when  inerts  of  different  dif- 
fusivities  are  added  to  either  the  fuel  or  the  oxi¬ 
dizer  stream. 

That  diffusional  stratification  may  influence  the 
soot  formation  process  can  be  demonstrated  by  the 
following  experimental  anomaly:  In  his  study  of  the 
influence  of  the  various  diluents  on  the  soot  for¬ 
mation  process,  McLintock'  varied  the  inert  dilu¬ 
ent  on  the  oxidizer  side  of  an  ethylene  coflow  flame 
and  found  that  the  smoke  points  decreased  in  the 
order  of  helium,  nitrogen  and  argon.  That  is,  the 
propensity  to  soot  was  least  when  the  oxidizer  dilu¬ 
ent  was  helium  and  greatest  when  the  diluent  was 
argon.  With  the  heat  capacities  of  argon  and  he¬ 
lium  being  the  same,  it  was  reasoned'  that  the  dif¬ 
ference  was  due  to  the  high  diffusivitv  of  helium, 
causing  the  concentration  of  helium  to  be  greater 
within  the  flame. 

An  O'  uosite  result,  however,  was  reported  by 
Schug  e»  al.,N  who  found  that  when  helium  was 
added  to  the  fuel  side  of  an  ethylene  coflow  dif¬ 
fusion  flame,  it  was  not  as  effective  as  argon  in  sup¬ 
pressing  soot.  In  fact,  helium  was  the  least  effective 
additive  for  soot  suppression  This  behavior  was  at¬ 
tributed  to  helium’s  high  thermal  diffusivitv  and  its 
subsequent  effect  on  the  temperature  distribution 
and  thereby  the  fuel  pyrolysis  rate. 

While  Refs.  7  and  8  have  respectively  used  con¬ 
centration  modification  and  temperature  modifica- 
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tion  to  interpret  the  different  trends  of  their  re¬ 
sults,  it  is  of  interest  to  note  that  these  observations 
can  also  be  explained,  in  a  consistent  manner,  as  a 
consequence  of  either  concentration  modification  or 
temperature  modification  alone.  For  example,  by 
using  the  temperature  modification  interpretation, 
the  difference  in  the  results  for  helium  could  be 
solely  a  consequence  ol  the  broadening  of  the  tem¬ 
perature  distribution  due  to  helium. y  That  is,  when 
helium  is  added  to  the  oxidizer  side,  the  broadened 
temperature  distribution  has  a  greater  influence  on 
the  rate  of  particle  burnup  and  a  lesser  influence 
on  the  rate  of  soot  formation,  implying  that  helium 
would  be  an  effective  soot  suppressor  relative  to  ar¬ 
gon.  On  the  other  hand  when  helium  is  added  to 
the  fuel  side,  the  temperature  broadening  would 
promote  the  soot  formation  rate  but  not  the  burn¬ 
up  rate,  implying  that  helium  would  be  less  effec¬ 
tive  in  suppressing  soot. 

We  next  adopt  the  alternate  viewpoint  that  con¬ 
centration  modification  could  be  the  physical  mech¬ 
anism  responsible  for  the  anomalous  results  for  he¬ 
lium.  Recognizing  that  helium  is  more  mobile  than 
argon,  it  is  clear  that  when  the  inerts  are  added  to 
the  fuel,  helium  will  he  defocused  to  a  greater  ex¬ 
tent  than  argon  due  to  the  compressive,  concave 
curvature  of  the  co-flow  flame.3  Consequently,  in 
all  hut  the  lowest  regions  of  the  flame,  the  fuel 
concentration  of  the  helium-diluted  mixture  will  be 
higher  than  that  of  the  argon-diluted  mixture,  im¬ 
plying  that  the  former  has  a  greater  propensity  to 
soot  By  applying  the  same  reasoning  to  the  oxi¬ 
dizer  side,  and  recognizing  that  helium  is  now  fo¬ 
cused  to  a  greater  extent  due  to  the  convex  flame 
curvature,  the  fuel  concentration  with  the  helium- 
diluted  oxidizer  will  be  relatively  lower,  implying 
that  the  propensity  to  soot  is  reduced. 

In  order  to  further  understand  the  relative  influ¬ 
ences  of  temperature  modification  versus  concen¬ 
tration  modification  on  the  soot  formation  process, 
we  have  adopted  an  entirely  different  flame  config¬ 
uration  to  study  the  soot  response  to  inert  addition 
The  flame  configuration  is  that  ol  the  counterflow 
diffusion  flame.  The  rationale  for  such  a  choice  is 
that  in  this  flame/floxv  configuration  the  effect  of 
concentration  modification  is  opposite  to  that  of 
temperature  broadening,  thus  allowing  unambigu¬ 
ous  conclusions  to  he  drawn  regarding  the  effect  ol 
preferential  diffusion.  We  shall  amplify'  on  this  (xiint 
later  in  the  discussion  of  the  experimental  results 

In  the  next  section  the  experimental  methodol¬ 
ogy  and  arrangement  will  be  specified  Ibis  is  fol¬ 
lowed  by  presentation  and  discussion  of  the  ex¬ 
perimental  results. 

Experimental  Specifications 

In  this  study  the  coimtcrflow  diffusion  flame  es¬ 
tablished  in  the  forward  stagnation  region  ol  ,1  po¬ 


rous  cylinder  has  been  used  The  flame  is  amenable 
to  light  scattering  and  extinction  measurements  and 
does  not  suffer  from  the  complication  ul  uxidi/ci 
entrainment  at  the  rim  as  with  jet  ico-flou  1  flames 
Measurements  are  taken  along  the  forward  stagna¬ 
tion  streamline,  which  is  a  line  of  svniiuetrx 

The  burner  apparatus  is  similar  to  that  used  l>\ 
Yandsburger  et  al . 10  and  has  been  described  m  de¬ 
tail  in  an  earlier  work. 11  Briefly  ,  it  consists  of  a  po¬ 
rous  cylinder  mounted  horizontally  in  a  vertically 
oriented  wind  tunnel  that  has  a  uniform  velocity 
profile  across  its  test  section.  The  cylinder  is 
mounted  to  a  vertical  translator  which  allows  it  to 
be  moved  in  the  streamwise  direction.  The  oxidizer 
gas  is  supplied  by  the  wind  tunnel  while  the  fuel 
gas  is  ejected  from  the  porous  cylinder  The  cyl¬ 
inder  is  composed  of  a  sintered  stainless  steel  tube 
with  an  outside  diameter  of  2  5  cm.  a  length  of  5 
cm,  and  is  cooled  uniformly  along  its  axis 

The  fuels  used  are  commercial  grade  ethylene  and 
instrument  grade  propane.  The  inerts  used  are 
krypton,  argon,  neon  and  helium 

Temperature  measurements  are  obtained  with 
silica-coated  Pt/Ft- l2%"Rh  thermocouples.  The  wire 
diameters  used  are  0.076  mm  for  the  diliited-pm- 
panc  flames  and  0.127  mm  for  the  dilutcd-cthxh-ue 
flames  as  well  as  flames  where  the  composition  of 
the  oxidizer  inert  is  varied.  The  smaller  thermo¬ 
couple  is  preferred  because  it  offers  less  distui- 
bailee  to  the  flame  and  less  measurement  crroi  due 
to  radiation  loss.  However,  the  temperatures  in  these 
latter  flames  are  high  enough  that  the  thinner  tin  r- 
mocouple  would  melt,  thus  necessitating  the  use  ol 
the  thicker  thermocouple. 

Soot  measurements  are  obtained  by  using  stan- 
dard  light  scattering  and  extinction  techniques  with 
an  apparatus  and  procedure  identical  to  those  of 
Flower  and  Bowman12  and  described  in  greater  dr- 
tail  in  Ref.  11.  Soot  volume  fraction  is  obtained  from 
the  Lambert-Beer  law  for  light  extinction  with  the 
assumption  that  the  soot  particles  behave  as  spher¬ 
ical  Rayleigh  scatterers  and  have  a  complex  refrac¬ 
tive  index  of  ri  =  1.57-i0.56.  The  soot  volume  frac¬ 
tion  is  assumed  uniform  over  the  length  of  the 
burner,  with  no  correction  for  end  effects  Since  this 
work  is  concerned  with  relative  changes  instead  of 
absolute  values,  the  Rayleigh  assumption  is  accept¬ 
able 

Gas  sampling  ol  the  stable  species  is  performed 
with  an  uncooled,  quartz  microprobe  designed  with 
the  criterion  of  Fristrom  and  Wcstcnhorg  '  '  To  al¬ 
low  tor  sampling  in  the  sooting  region  of  the  flame 
a  probe  was  designed  with  an  internal  plunger11  that 
allows  soot  build-up  at  the  probe  orifice  to  be  in 
moved  periodicals  The  orifice  diameter  is  70  jjliii 
and.  to  ensure  that  reactions  are  quern  lied,  sam 
plmg  is  conducted  at  pressures  such  that  the  flow 
through  the  orifice  is  choked  I  he  gas  samples  ,ue 
obtained  by  batch  sampling  into  ex.uuatcd  vessels 
and  the  gas  analysis  is  performed  In  quantitative 


D-2 


SOOTING  COUNTERFLOW  DIFFUSION  FLAMES 


mass  spectrometry  Since  batch  sampling  is  used, 
data  are  not  available  for  water  and  results  are  re¬ 
ported  on  a  dry  basis. 

A  standard  single-component,  laser- Doppler  ve- 
loeinietrv  (LDV)  system  is  used  to  measure  the  ve¬ 
locity  component  along  the  forward  stagnation 
streamline.  Measurements  are  made  in  the  forward 
scattering  mode,  without  frequency  shifting.  Alu¬ 
minum-oxide  seed  particles  of  0.05  pm  nominal  di¬ 
ameter  are  supplied  on  the  oxidizer  side. 

The  soot  loading  of  the  propane  diffusion  flame 
in  air  is  too  low  to  measure  extinction  accurately. 
To  increase  the  soot  loading  of  propane  flames  the 
oxygen  volume  fraction  in  the  oxidizer  is  increased 
from  0.21  to  0.25,  with  the  balance  being  nitrogen. 
Ethylene  flames  do  not  require  enrichment  and  their 
oxygen  concentrations  are  maintained  at  that  of  air. 

For  the  experimental  conditions  in  this  study  the 
flame  iocation  is  a  few  millimeters  on  the  oxidizer 
side  of  the  stagnation  point  (Fig.  1).  The  soot  for¬ 
mation  and  growth  region  lies  between  the  flame 
and  the  region  of  the  stagnation  point.  A  conse¬ 
quence  of  this  is  that  it  is  the  flow  field  on  the 
oxidizer  side  that  is  of  primary  importance  tor  res¬ 
idence  time  considerations.  Therefore,  the  up¬ 
stream  velocity  of  the  oxidizer  is  maintained  at  58 
em/s  for  all  studies. 

An  additional  consideration  is  the  gas  flow  ejected 
from  the  burner  Two  approaches  are  used  in  this 
study.  One  is  to  maintain  the  volumetric  flow  rate 
of  the  fuel  constant.  In  this  case  the  addition  of  gas 
to  the  fuel  increases  the  ejection  velocity  from  the 
burner  The  other  approach  is  to  maintain  the  fuel 


Fit;.  1  Cross  section  of  the  counterflow  diffusion 
flame  established  in  the  forward  stagnation  region 
ol  a  porous  cvfinder. 


ejection  parameter  constant  15  The  fuel  ejection  pa¬ 
rameter  defined  lor  variable  gas  density  is  f„  — 
-fPui'.,/p*KKeJ'1/2.  where  p„  and  i„  are  the  gas 
density  and  velocity  at  the  burner  surface .  p,  and 
v,  the  freestreain  values  of  density  and  kinematic 
viscosity  for  the  oxidizer,  and  k  the  characteristic 
velocity  gradient  defined  as  k  =  2L3./R.  with  l  , 
being  the  freestreain  velocity  of  the  oxidizer  and  R 
the  burner  radius.  The  fuel  ejection  parameter  de¬ 
fined  in  this  way  allows  flow  similarity  to  be  main¬ 
tained  even  when  variable  density  inerts  are  added 
Only  results  for  constant  fuel  ejection  parameter.  fu 
=  —3.0,  will  be  presented  herein  because  main¬ 
taining  flow  similarity  allows  for  an  interpretation 
based  on  distance  from  the  burner  It  is  of  interest 
to  note  that  the  results  were  found  to  be  indepen¬ 
dent  of  whether  fuel  flow  rate  or  fuel  ejection  pa¬ 
rameter  is  held  constant 

Results  and  Discussions 

As  discussed  in  the  Introduction,  both  concen¬ 
tration  modification  and  temperature  broadening  will 
yield  similar  trends  for  inert  addition  to  the  coflow 
flame.  Therefore  it  is  difficult  to  determine  with  the 
coflow  flame  which,  if  either,  mechanism  is  the 
dominant  cause  lor  the  observed  trends.  This  is  not 
the  case  with  the  countcrflovv  diffusion  flame  Rv 
referencing  Fig.  1  it  can  be  seen  that  the  fuel 
reaches  the  flame  by  diffusing  across  the  stagnation 
plane  into  the  oxidiz.er  side.  The  streamlines  on  both 
the  fuel  and  oxidizer  sides  of  the  stagnation  plane 
are  such  that  lighter  species,  which  can  more  easily 
diffuse  across  the  streamlines,  will  preferentially  ac¬ 
cumulate  in  the  region  of  the  stagnation  streamline, 
while  heavier  species  will  be  preferentially  swept 
away.  Consequent)),  when  helium  is  added  to  the 
fuel,  concentration  modification  should  lead  to  a  re¬ 
duction  in  fuel  concentration  and.  thereby .  a  cor¬ 
responding  reduction  in  the  soot  loading  as  com¬ 
pared  to  the  situation  of  adding  argon  to  the  fuel 

U’c  further  note  that  soot  burn-up  does  not  oc¬ 
cur  in  the  forward  stagnation  region  ol  the  coun¬ 
terflow  diffusion  flame  Therefore,  il  helium  addi¬ 
tion  were  to  broaden  the  temperature  distribution 
the  effect  would  be  to  increase  the  soot  loading  .is 
based  on  the  reasoning  of  Ref  5  discussed  in  the 
Introduction.  Thus,  in  the  countcrflovv  diffusion 
flame,  the  effect  of  concentration  modification  is 
opposite  to  that  of  broadening  of  the  temperature 
distribution  and  we  can  compare  the  relative  im¬ 
portance  of  the  two  mechanisms  for  this  flame 

In  addition  to  helium  and  argon  we  have  also  used 
neon  and  krvpton  for  diluents  in  this  study.  Ml  di¬ 
luents  are  therefore  completely  inert  and  Mine  tliev 
are  monatomic  and  have  the  same  heat  capacitv. 
they  should  have  a  similar  influence  on  the  maxi¬ 
mum  flame  temperature 

The  soot  volume  fraction  drops  oil  around  the 
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stagnation  point  because  ol  the  inability  of  the  soot 
particles  to  propagate  past  the  stagnation  point  into 
the  fuel  side  (see  big.  1).  Data  in  this  region  have 
been  suppressed  as  they  are  not  relevant  to  this 
work. 

The  soot  volume  fraction  distributions  along  (he 
stagnation  streamline  for  the  propane/inert  mix¬ 
tures  are  plotted  in  Fig.  2.  The  curves  show  a  strong 
relationship  between  the  mobility  of  the  inert  and 
the  maximum  soot  loading  for  the  flame.  Krypton, 
being  the  least  mobile  inert,  yields  the  greatest  soot 
loading  (0.52  x  10~fi)  while  helium,  being  the  most 
mobile,  yields  the  lowest  soot  loading  (0.2fi  x  10  h). 

The  results  for  inert  addition  to  ethylene,  shown 
in  Fig.  3,  display  trends  similar  to  those  ol  the  pro¬ 
pane  flames  in  both  their  relationship  to  the  inert 
mobility  and  in  the  extent  ol  the  change  in  soot 
volume  fraction,  the  soot  loading  increases  with  de¬ 
creasing  inert  mobility,  and  the  maximum  soot  vol¬ 
ume  fraction  of  the  krypton-diluted  flame  is  about 
twice  as  high  as  that  of  the  helium-diluted  flame. 

The  above  results  indicate  that  the  soot  loading 
varies  with  the  mobility  of  the  inert  added.  To  de¬ 
termine  whether  the  primary  cause  for  the  change 
m  the  soot  loading  is  indeed  modification  of  the  fuel 
concentration,  we  have  determined  the  tempera¬ 
ture.  velocity,  and  species  concentration  profiles  for 
these  flames. 

Figure  -4  shows  the  temperature  distributions, 
uncorrected  for  radiation  losses,  for  the  propane/ 
inert  mixtures.  The  oxygen  concentration  for  the 
temperature  measurements  is  that  of  air,  or  21  per¬ 
cent  by  volume,  as  compared  to  25  percent  for  the 
soot  volume  fraction  measurements.  While  the  trends 


Fit;  3.  Soot  volume  fraction  as  a  function  of  dis¬ 
tance  from  burner  surface  for  ethv  Icne/air  cotui- 
terflosv  diffusion  flames  with  inert  addition  on  the 
fuel  side 

iii  maximum  temperature  and  temperature  distri¬ 
bution  will  not  lie  alleeted  bv  the  lower  oxygen 
concentration,  the  soot  loading  is  much  reduced 
Coiiso(|ucntly,  the  accuracy  and  icproducibilitv  ol 
the  measurements  are  greatly  unproved  due  to  the 
fact  that  soot  Iwild-up  on  the  thermocouple  iv  less 
The  maximum  flame  temperatures  measured  show 
slight  increases  with  decreases  in  the  mobility  (in¬ 
creases  in  the  molecular  weight)  of  the  inert,  kry  p¬ 
ton  yields  the  highest  temperature  il641cCl,  he- 


Fli;.  2  Soot  volume  fraction  as  a  function  of  dis¬ 
tance  from  burner  surface  for  propane/enriched  air 
.25°T  Op  counterflow  dilliision  flames  with  inert  ad¬ 
dition  nil  the  fuel  side 


l'T( .  1  Temperature  distribution  lor  propane /an 

counterflow  diffusion  flames  with  inert  addition  on 
the  fuel  side. 


D-4 


SOOTING  COl'NTKRKLOW  Dll-  Kl  SION  FLAMES 


Hum  the  lowest  (1607°  (7).  with  a  difference  of  34°  C. 
The  temperature  profiles  are  similar  with  no  sig¬ 
nificant  broadening  for  the  lighter  inerts. 

One  may  ask.  at  this  point  whether  the  34°  C 
difierence  in  the  maximum  measured  temperature 
seen  in  Fig  4  could  account  for  the  two-fold  dif¬ 
ference  in  soot  loading  The  influence  of  this  tem¬ 
perature  difference  has  been  investigated  by  vary¬ 
ing  the  flame  temperature  in  such  a  way  that  the 
reactive  species  concentrations  are  not  affected.11 
The  results  show  that  even  with  a  43°  C  increase 
in  the  maximum  measured  temperature  the  soot 
loading  increases  by  only  13%. 

For  the  diluted  ethylene  flames  the  maximum 
flame  temperatures  were  measured  and  were  found 
to  he  within  20°  C  for  flames  diluted  with  each  of 
the  four  inerts.  Considering  this  along  with  Fig.  4 
it  is  reasonable  to  expect  that  temperature  effects 
have  only  a  minor  role  in  explaining  how  the  inerts 
are  influencing  the  soot  volume  fractions. 

We  next  note  that  there  is  the  possibility  that 
the  velocity  fields  vary  for  the  different  flames  and 
consequently  the  residence  times  are  not  similar. 
■Since  the  soot  formation  and  growth  region  lies 
nominally  between  the  flame  and  the  stagnation 
point,  it  is  the  velocity  field  in  this  region  that  in¬ 
fluences  the  residence  times  of  the  soot  particles. 
Figure  5  is  a  plot  of  velocity  along  the  stagnation 
streamline  versus  distance  from  the  burner  surface 
tor  the  propane  flame  when  diluted  with  the  inerts. 
The  velocity  fields  are  almost  identical,  which  is 
reasonable  because  the  temperature  fields  are  al¬ 
most  identical  and  we  have  maintained  fluid-dy¬ 
namic  similarity  by  holding  2U./H  and  the  fuel 
ejection  parameter  /„.  constant. 


A  close  inspection  of  Figs  2-3  reveals  that  the 
locations  ol  the  flame  and  stagnation  point  are  shifted 
away  from  the  burner  surface  as  the  mobilitv  of  the 
inert  is  increased.  The  shift  in  stagnation  point  is 
most  accurately  indicated  by  the  locations  ol  max¬ 
imum  soot  in  Figs.  2  and  3  while  the  shift  in  flame 
location  is  indicated  by  both  peak  temperatures  and 
velocities  in  Figs.  4  and  5.  The  shifts  are  small 
(<l/2  mm)  and  will  not  influence  soot  particle  res¬ 
idence  times  because  the  flame  and  stagnation  point 
shift  together. 

This  shift  can  be  understood  In  considering  (In¬ 
tact  that  in  these  experiments  wc  have  maintained 
/„  and  k  fixed  so  that,  in  the  nonreactive  limit.  (In¬ 
flow  fields  are  similar  in  terms  of  -q  —  us/p.p.,  1 
JKpdi/  Iherefore.  while  the  stagnation  points  occur 
at  the  same  T|  for  all  flows,  they  are  shifted  in  tin- 
physical  coordinate  y  For  example,  since  the  den¬ 
sity  is  less  in  the  helium-diluted  flame,  the  phvsical 
distance  to  the  stagnation  point  is  greater  Esti¬ 
mates  using  mean  densities  show  agreement  be¬ 
tween  the  anticipated  and  observed  shifts. 

With  the  velocity  and  temperature  fields  being 
largely  independent  of  the  additive  mobilitv  the 
implication  is  that  the  local  concentrations  of  the 
reactive  species  have  been  modified.  To  verify  this, 
a  measurement  of  the  species  concentrations  within 
the  sooting  flames  is  obtained.  Probe  sampling  of 
the  stable  species  is  performed  at  conditions  similai 
to  those  ol  the  temperature  measurements,  in  that 
the  oxidizer  is  21  percent  oxygen  as  compared  to 
2-5  percent  used  for  the  soot  data.  The  lower  oxi¬ 
dizer  concentration  reduces  the  extent  of  sooting 
while  preserving  the  important  concentration  fea¬ 
tures. 


|-|C  5  Yclocitv  as  a  function  of  distance  from  the 
humci  surface  for  propane/ciu iehed  air  (25%  0.1 
counterflow  diffusion  flames  with  inert  addition  on 
the  fuel  side 
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Fig.  7.  Mole  fraction  of  parent  fuel  as  a  function 
of  distance  from  the  flame  for  propane/air  coun¬ 
terflow  diffusion  flames  with  inert  addition  on  the 
fuel  side. 

Figure  6  compares  the  mole  fraction  profiles  of 
helium  and  argon  for  the  helium-  and  argon-diluted 
flames,  respectively.  The  mole  fraction  of  the  inert 
diluent  in  the  flame  region  is  higher  in  the  helium- 
diluted  flame  than  in  the  argon-diluted  flame,  a 
consequence  of  the  disparity  in  the  inert  mobilities. 

Since  the  local  concentrations  of  the  additive  dif¬ 
fer,  the  local  fuel  concentrations  could  be  affected. 
The  relevant  coordinate  to  compare  propane  con¬ 
centrations  is  distance  from  the  flame,  rather  than 


distance  from  the  burner  because  concentration 
gradients  in  the  flame  region  are  steep  and  the  flame 
locations  are  not  precisely  the  same.  In  Fig.  7  the 
mole  fraction  of  the  parent  fuel,  propane,  is  plotted 
as  a  function  of  distance  from  the  flame.  Near  the 
flame,  the  concentration  gradients  are  found  to  be 
less  for  the  helium-diluted  flame  which,  in  itself, 
delineates  that  the  mole  fraction  of  the  fuel  will  be 
locally  less,  and  Fig.  7  illustrates  this.  However,  a 
caveat  is  necessary  with  regard  to  interpreting  Fig 
7.  The  concentration  gradients  are  so  steep  near  the 
flame  that  spatial  measurement  errors  of  just  0. 1 
mm  could  influence  the  conclusions  drawn.  Al¬ 
though  care  was  exercised  to  maximize  the  spatial 
accuracy,  there  is  an  alternative  means  of  demon¬ 
strating  the  effects  of  preferential  diffusion  that  is 
not  as  sensitive  to  spatial  accuracy. 

Concentration  modification  of  the  fuel  will  nec¬ 
essarily  lead  to  a  modification  of  the  concentrations 
of  the  fuel-based  species  in  the  flame  and  these 
species  will  attain  a  maximum  somewhere  within 
the  flame.  The  maximum  is  insensitive  to  location 
and,  where  applicable,  is  a  good  indicator  ot  the 
relative  availability  of  soot  precursors.  Table  I  lists 
the  maximum  mole  fractions  for  all  measured  spe¬ 
cies  except  for  ethane  and  propadiene  which  did 
not  display  obvious  maximums.  Without  exception, 
the  mole  fractions  are  higher  in  the  argon-diluted 
flame,  and,  in  the  case  of  potentially  important  soot- 
particle  precursors,  notably  butadivne  (diacetylene), 
butadiene  and  benzene,  the  mole  fractions  are  15- 
25%  higher.  The  results  demonstrate  that  the  local 
mixture  composition  is  indeed  modified  in  the  pres¬ 
ence  of  preferential  diffusion  and  this  modification 
can  influence  soot  formation. 


TABLE  I 

Maximum  mole  fraction  of  stable  species  in  propane/air 
counterflow  diffusion  flames  with  inert  addition  on  the  fuel  side. 


Fuel  Mixture 

Maximum  Mole  Fraction 

Hydrogen 

’  (H2) 

Carbon 

Monoxide 

(CO) 

Carbon 

Dioxide 

(CO,) 

Methane 

(CH4) 

Acetylene 

(C..IL) 

Ethylene 

<C,H4! 

20%  He/80%  C,Hh 

.038 

.052 

.086 

.017 

.018 

030 

20%  Ar/80%  CjH* 

.042 

.057 

.091 

.018 

020 

032 

%  Difference 

10 

9 

6 

6 

10 

6 

1-3 

1-3 

Cvclopen- 

Propene 

Butadiyne 

Butadiene 

2-Butene 

tadiene 

Benzene 

<C,H6) 

(C4H2) 

(C,hb) 

(C,H,) 

(C5H„) 

(C6Hfi) 

20%  He/80%  QH, 

.0066 

.00042 

.00066 

.00063 

.00021 

00048 

20%  Ar/80%  C5H» 

.0070 

.00049 

.00089 

.00068 

.00024 

.00060 

%  Difference 

6 

14 

26 

7 

12 

20 
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Fig.  8.  Temperature  distribution  for  propane/air 
counterflow  diffusion  flames  with  20%  of  the  nitro¬ 
gen  in  the  oxidizer  replaced  with  argon  or  helium. 


however,  is  coupled  to  the  inert  concentration  so 
that  stratification  of  the  inert  can  have  a  significant 
effect  on  the  concentrations  of  fuel-related  species. 

The  effects  of  inert  mobility  on  the  oxidizer  side 
of  a  propane  flame  are  considered  in  Figs.  8  and 
9.  Twenty  percent  of  the  nitrogen  in  the  oxidizer 
is  replaced  with  either  argon  or  helium  so  that  the 
concentration  of  these  test  inerts  in  the  oxidizer  is 
15  percent.  Figure  8  reveals  that  the  maximum 
temperatures  of  the  helium-  and  argon-diluted  flames 
are  still  quite  similar  even  though  the  helium-di¬ 
luted  flame  shows  some  broadening  of  the  temper¬ 
ature  distribution.  Furthermore,  the  maximum  soot 
volume  fractions  seen  in  Fig.  9  are  almost  identical 
Referring  to  Fig.  1  it  is  seen  that  the  transport  of 
fuel  from  the  stagnation  point  to  the  flame  is  totally 
dependent  on  diffusion  while  the  transport  of  the 
oxidizer  to  the  flame  is  assisted  by  convection. 
Therefore  it  is  reasonable  that  the  mobility  of  the 
inert  added  to  the  fuel  would  have  a  strong  influ¬ 
ence  on  fuel  concentrations  while  adding  the  same 
inert  to  the  oxidizer  side  would  have  little  effect  on 
fuel  concentrations. 


It  may  appear  somewhat  surprising  that  the 
stratification  of  the  inert  concentrations  would  affect 
the  composition  without  affecting  the  temperature 
distributions.  Although  the  thermal  diffusivities  of 
the  inerts  used  vary  widely,  these  additives  are  only 
a  small  portion  of  the  total  gas  composition  in  the 
region  of  the  flame.  Consequently,  the  thermal 
transport  properties  are  still  dominated  by  the  ma¬ 
jor  species,  e  g.  nitrogen.  The  fuel  concentration, 


Distance  from  Burner  (mm) 

Fig  9  Soot  volume  fraction  as  a  function  of  dis¬ 
tance  from  burner  surface  for  propane/enriched-air 
(25%  CL)  counterflow  diffusion  flames  with  20%  of 
the  nitrogen  in  the  oxidizer  replaced  with  argon  or 
helium 


Concluding  Remarks 

In  the  present  investigation  we  have  provided 
experimental  evidence  that  gives  strong  support  to 
the  concept  that  soot  formation  in  diffusion  flames 
can  be  significantly  affected  by  concentration  mod¬ 
ification  due  to  preferential  diffusion  of  the  reac- 
tive/inert  species.  We  have  further  shown  that, 
white  this  result  is  reasonable  and  possibly  not  un¬ 
expected,  a  convincing  demonstration  of  this  fact 
requires  careful  consideration  and  ultimately  elim¬ 
ination  of  other  possible  causes,  such  as  those  due 
to  the  potential  modification  of  the  temperature  and 
velocity  fields. 

Recognizing  the  existence  and  importance  of  dif- 
fusional  stratification  of  the  concentrations  of  the 
reactive  mixtures,  it  would  be  of  interest  to  extend 
the  present  investigation  to  diffusion  flames  of  other 
configurations,  to  premixed  flames,  and  to  the  for¬ 
mation  of  other  pollutant  species  which  are  sensi¬ 
tive  to  concentration  modifications. 
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COMMENTS 


R  J.  Santoro,  Pennsylvania  State  Univ..  USA.  In 
the  presentation,  you  indicated  that  the  tempera¬ 
ture  was  held  constant  to  within  about  40K.  Could 
the  variation  in  temperature  account  for  the  small 
changes  (approximately  20%)  you  observed  in  the 
species  concentration. 

Author's  Reply.  The  concentration  modification 
resulting  from  preferential  diffusion  influences  flame 
temperature,  and  any  difference  in  temperature 
could,  in  turn,  affect  the  concentrations  of  the  py¬ 
rolysis  products.  For  the  helium-  and  argon-diluted 
flames  the  measured  difference  in  maximum  tem¬ 
perature  is  30°C  and  a  similar  temperature  differ¬ 
ence  is  iound  throughout  the  pyrolysis  region.  (Note 
that  if  there  was  a  noticeable  broadening  of  the 
temperature  distribution  due  to  helium  addition, 
the  temperature  difference  of  the  two  flames  would 
have  been  even  less  in  the  pyrolysis  region.) 

This  difference  in  temperature  certainly  plays  a 
role  in  the  concentration  differences,  however  we 
do  not  believe  that  temperature  alone  can  account 
for  the  measured  differences  in  either  soot  or  spe¬ 
cies  concentrations.  While  the  authors  are  not  aware 
of  any  data  describing  how  flame  temperature  in- 
lluences  species  concentrations  for  this  flame,  some 
insight  can  be  gained  by  considering  results  where 
the  influence  of  temperature  on  soot  concentration 
of  a  propane  flame  was  studied.1 '  When  flame  tem¬ 
perature  was  increased  43°C  without  changing  either 


the  fuel  or  oxidizer  concentration,  the  maximum  soot 
concentration  increased  only  13%.  far  less  than  the 
almost  60%  increase  observed  between  the  helitim- 
and  argon-diluted  flames  of  the  present  study.  The 
variation  in  species  concentrations  resulting  from  a 
30°C  increase  in  temperature  would  not  appear  to 
account  for  the  observed  variations  in  soot  concen¬ 
tration. 

The  species  measurements  shown  in  Table  1  and 
Fig.  7  also  imply  that  factors  other  than  tempera¬ 
ture  are  responsible  for  the  measured  difference). 
Not  only  are  the  products  and  by-products  of  the 
fuel  greater  in  the  argon-diluted  flame,  the  parent 
fuel  concentration  is  also  locally  greater — an  un¬ 
likely  result  il  temperature  is  responsible  for  the 
concentration  differences. 
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EXPERIMENTS  ON  THE  SOOTING  LIMITS  OF 
AERODYNAMIC  ALLY-STRAINED  DIFFUSION  FLAMES 
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University  of  California 
Davis,  California  95616 


An  experimental  study  has  been  performed  with  axisymmetric  eounterflow  diffusion  (lames 
to  investigate  the  influence  of  aerodynamic  straining  on  the  relevant  sooting  limits  of  the 
lower  alkanes.  The  limits  are  defined  by  the  critical  strain  rate  at  which  either  soot  lumi¬ 
nosity,  soot  particle  light-scattering,  or  fluorescence  is  negligible  compared  to  the  appropriate 
background  signal.  The  critical  strain  rates  are  found  to  be  greatest  for  the  sooting  limit  based 
on  the  fluorescence  signal,  with  those  based  on  luminosity  and  light-scattering  being  similar 
The  fluorescence  signal,  if  attributed  to  polycyclic  aromatic  precursors,  yields  a  limit  that 
can  be  interpreted  as  the  extinction  of  soot  precursors  and  is  suggested  to  he  a  possible  limit 
for  identifying  a  completely  nonsooting  flame  condition.  The  separate  effects  of  flame  tem¬ 
perature  and  fuel  concentration  on  the  critical  strain  rates  for  soot  extinction  have  also  been 
studied.  The  results  are  indicative  of  how  temperature  and  concentration  influence  the  soot 
particle  inception  process  and  they  show  that  both  are  potentially  important  parameters.  The 
critical  strain  rates  display  an  Arrhenius  temperature  dependence  and  this  dependence  is 
similar  for  all  alkanes  considered. 


Introduction 

Soot  formation  is  a  rate  process,  and  as  such  is 
subjected  to  residence  time  considerations.1,2  Since 
most  flames  are  under  the  influence  of  aerodynamic 
straining,  the  relevant  time  scale  for  such  flames  is 
the  rate  of  strain  experienced  by  the  flame,  ex¬ 
pressed  in  the  unit  of  inverse  time.3  Indeed,  Tsuji 
and  Yamaoka’1  have  shown  that  disappearance  of  the 
luminous  yellow  zone  in  a  counterflow  diffusion  flame 
can  be  achieved  by  increasing  the  strain  rate  of  the 
approach  flow. 

The  present  study  aims  to  provide  a  detailed  ex¬ 
perimental  investigation  of  the  influence  of  aero¬ 
dynamic  straining  on  particle  formation/extinction 
along  the  following  directions.  First,  we  note  that 
the  strain  rate  of  Ref.  4  is  a  global  one,  valid  in 
the  limit  of  an  inviscid,  incompressible,  nonre ac¬ 
tive,  unbounded  flow.  It  is  reasonable  to  expect  that 
the  actual,  local  strain  rate  as  experienced  by  the 
flame  can  be  substantially  different  from  the  global 
definition.  We  shall  therefore  use  the  local  strain 
rate,  determined  by  laser  Doppler  velocimetry,  as 
the  parameter  to  quantify  the  various  soot  extinc¬ 
tion  phenomena  to  be  studied. 

Our  second  objective  is  to  understand  and  quan¬ 
tify  the  definition  of  the  sooting  limit,  i.e.,  soot  in¬ 
ception/extinction  limit,  in  diffusion  flames  and  the 
various  sub-processes  associated  with  it,  using  the 


local  extinction  strain  rate  as  the  quantifying  pa¬ 
rameter.  Specifically,  it  is  of  interest  to  explore 
whether  there  is  a  close  relation  between  the  visual 
disappearance  of  the  yellow  luminous  region  and 
the  actual  disappearance  of  the  soot  particles  This 
is  a  relevant  concern  because  for  diffusion  flames 
the  temperature  at  the  location  of  the  peak  soot 
volume  fraction  is  typically  many  hundreds  of  de¬ 
grees  lower  than  the  flame  temperature.5  6  There¬ 
fore  soot  luminosity,  with  it’s  T3  dependence,  may 
not  be  an  accurate  indication  of  soot  extinction  in 
these  flames.  Furthermore,  because  of  this  tem¬ 
perature  dependence,  luminosity  may  not  be  ap¬ 
propriate  for  studying  the  relative  effects  of  tem¬ 
perature  on  the  sooting  limits. 

An  additional  consideration  with  the  definition  of 
soot  extinction  is  that  since  the  formation  of  soot 
particles  is  preceded  by  the  formation  of  soot  pre¬ 
cursors,  absence  of  the  particles  does  not  necessar¬ 
ily  imply  the  corresponding  absence  of  the  precur¬ 
sors  and  thereby  the  inability  to  soot.  We  shall 
therefore  also  provide  experimental  data  and  evi¬ 
dence  which  attempt  to  distinguish  between  the 
sooting  limits  of  the  particles  and  precursors. 

Having  defined  the  strain  rate  and  the  sooting 
limits,  our  third  objective  is  to  use  these  definitions 
to  quantify  the  influences  of  temperature  and  inert 
addition  to  the  fuel  on  the  sooting  limits.  Since  in¬ 
ert  addition  changes  both  the  flame  temperature  as 
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well  as  the  fuel  concentration,  it  is  necessary  to  iso¬ 
late  these  two  factors  and  separately  identify  their 
individual  effects,  as  has  been  done  in  Ref.  6.  In 
this  previous  work  it  was  suggested  that  fuel  dilu¬ 
tion  can  have  a  significant  influence  on  soot  volume 
fraction  by  changing  the  total  surface  area  in  the 
particle  inception  zone.  Since  the  inception  limit  is 
concerned  with  whether  particles  are  produced,  not 
how  many  are  produced,  it  is  not  clear  a  priori 
whether  dilution  affects  the  inception  limit  as  well. 
Furthermore,  while  the  importance  of  temperature 
on  soot  formation  is  well  documented  as  it  relates 
to  soot  yield,  ‘  there  are  no  known  studies  on  how 
relevant  flame  temperature  is  to  the  sooting  limits. 
In  our  study,  temperature  adjustment  will  be  ac¬ 
complished  by  substituting  a  portion  of  the  nitro¬ 
gen  in  the  oxidizer  with  an  equal  quantity  of  argon. 
With  this  technique  the  flame  temperature  can  be 
adjusted  without  affecting  the  concentration  of  the 
reactants  on  either  the  fuel  or  oxidizer  sides  of  the 
flame. 

The  experimental  methodology  is  outlined  in  the 
next  section,  which  is  followed  by  presentation  and 
discussion  of  the  experimental  results 


Flo.  I  Schematic  ot  the  burner  coiitiguration  ami 
the  optical  system  used  to  measure  light  scattering 


Experimental  Methodology 

The  opposed-jet  counterflow  diffusion  flame  is 
employed  in  this  study.  This  flame  is  stable  and  has 
a  simple,  well-defined,  axisym metric  flow  field. s  The 
flame  is  amenable  to  light-scattering  measurements 
and  does  not  suffer  from  the  complication  of  oxi¬ 
dizer  entrainment  at  the  rim  that  occurs  in  coflow 
flames.  Measurements  are  taken  along  the  center- 
line,  which  is  the  stagnation  point  streamline. 

The  burner  configuration  is  shown  in  Fig.  1.  The 
tlairie  is  stabilized  in  the  region  of  the  stagnation 
plane  established  between  two  converging  nozzles 
(following  nitrogen  is  employed  to  eliminate  oxi¬ 
dizer  entrainment  and  extinguish  the  outer  flame. 

Figure  1  also  shows  the  optical  system  employed 
to  measure  light  scattering  The  beam  of  an  argon- 
ion  laser  operating  at  488  tint  is  modulated  by  a 
mechanical  chopper  and  then  focused  with  a  300 
mm  focal  length  lens  at  the  midspan  of  the  burner. 
After  passing  through  a  standard  set  of  collection 
optics,  the  detection  of  scattered  light  is  accom¬ 
plished  with  a  photomultiplier  tube  whose  output 
is  processed  by  a  lock-in  amplifier  interfaced  to  a 
microcomputer 

Close  to  the  sooting  limit  there  are  three  signif¬ 
icant  sources  of  light  scattering  in  these  diffusion 
flames:  gaseous  Rayleigh  scattering,  particulate 
Rayleigh  scattering,  and  broadband  fluorescence. 
Since  the  species  responsible  for  the  fluorescence, 
assumed9- 11  to  be  polycyclic  aromatic  hydrocarbons 
(PAH),  are  usually  associated  with  the  soot  precur¬ 


sor  molecules, 10  ,2  11  then  the  disappearance  of  the 
fluorescence  signal  may  be  an  indication  of  the  limit 
at  which  precursor  extinction  has  been  achieved. 
Therefore,  determining  a  viable  method  of  discrim¬ 
inating  between  the  three  sources  of  scattering  is 
an  important  aspect  of  this  work:  the  details  are 
discussed  in  the  Results  section. 

Temperatures  are  determined  with  BeO-VM)* 
coated  Pt-20%Rh/Pt-40%Rh  thermocouples  with  a 
0.05  mm  wire  diameter.  The  lead  wires  are  ori¬ 
ented  parallel  to  the  flame  to  minimize  conduction 
errors.  Radiation  corrections  are  performed  by  as¬ 
suming  a  convection/radiation  balance  on  a  spher¬ 
ical  bead  with  a  value  of  0.38  for  the  cmissivity.  1 
in  the  calculation,  the  Nusselt  number  is  assumed 
to  be  constant  and  equal  to  2.0. 

A  standard  single-component,  laser- Doppler  \e- 
locimetry  (LDV)  system  is  used  to  measure  the  ve¬ 
locity  along  the  forward  stagnation  streamline  The 
system  is  operated  in  the  forward  scattering  mode, 
without  frequency  shifting  Aluminum-oxide  seed 
particles  with  a  1.0  pm  nominal  diameter  are  sup¬ 
plied  on  the  oxidizer  side 

For  the  present  axisymmetric  counterflovv  flame 
we  define  the  strain  rate  K  as  the  negative  ol  the 
measured  velocity  gradient,  -di/thj.  ahead  of  the 
flame.  In  Fig.  2  a  typical  curve  showing  the  mea¬ 
sured  velocity  field  along  the  stagnation  streamline 
is  shown.  The  data  is  taken  on  the  oxidizer  side  ot 
the  stagnation  point  where  the  flame  and  soot  for¬ 
mation  regions  occur  The  velocity  gradient  is  found 
by  fitting  the  data  upstream  of  the  flame. 
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aerodynamically-strained  diffusion  flames 


FlC.  2.  A  typical  velocity  profile  for  the  axisym- 
inetric  counterflow  diffusion  flame,  illustrating  the 
method  of  determining  the  local  strain  rate  from  the 
streamwise  velocity  gradient  ahead  of  the  flame. 

Results  and  Discussions 

Particle-Inception  and  Fluorescence  Limits: 

In  obtaining  a  particle-inception  limit,  our  inten¬ 
tion  is  to  find  the  critical  strain  rate  where  particle 
inception  is  so  small  that  light  scattering  from  soot 
particles  is  negligible  compared  to  background  scat¬ 
tering.  Near  this  limit  there  are  three  significant 
sources  of  light  scattering:  Rayleigh  scattering  due 
to  gases  and  particles,  and  fluorescence,  presum¬ 
ably  due  to  PAH’s.  The  soot  particles  are  spherical 
in  this  early  stage  of  their  life  and  therefore,  as  with 
Rayleigh  scattering  from  gases,  the  light  remains 
polarized  Furthermore,  since  the  scattering  is  elas¬ 
tic,  the  wavelength  is  unchanged  On  the  other  hand 
the  inelastic  fluorescence  signal  is  depolarized  and 
broadband  with  both  stokes  and  anti-stokes  fluo¬ 
rescence.15"1,  In  fact,  for  a  488  nm  excitation  wave¬ 
length,  the  spectral  peak  of  the  fluorescence  signal 
is  located  close  to  the  excitation  wavelength  A„,  with 
a  large  contribution  of  the  signal  being  at  A  =  488 
nm.  Consequently,  near  the  sooting  limit  the  flu¬ 
orescence  signal  can  be  of  the  order  of  the  Rayleigh 
signal  even  though  a  laser  line  filter  and  polariza¬ 
tion  analyzer  are  employed  ahead  of  the  detector 
to  diminish  the  fluorescence  signal  (Fig.  1). 

To  obtain  a  true  particle  limit  it  is  necessary  to 
discriminate  between  the  Rayleigh  signal  and  the 
residual  fluorescence  signal.  To  do  this  we  exploit 
the  fact  that,  unlike  Rayleigh  scattering,  fluores¬ 
cence  is  depolarizing.  Total  depolarization  of  the 
fluorescence  signal  has  been  checked  by  measuring 
the  vertical  (S, ,)  and  horizontal  (S /,v)  scattering 


components  with  an  edge  filter  used  to  remove  the 
elastic  scattering  signal.  It  was  found  that  S, ,  =  S,-„ 
for  A  >  530  nm.  Then  it  is  reasonable  to  assume 
that  S,,  =  S/„  at  A  =  488  nm,  and  a  measurement 
of  S;„(488  nm)  will  reveal  how  much  of  S,,(488  nm) 
is  due  to  fluorescence. 

Figure  3  illustrates  the  method  of  obtaining  the 
particie  inception  limit.  First  a  spatial  profile  of 
S,,,,(488  nm)  is  obtained  (Fig.  3a).  This  signal  con¬ 
tains  contributions  due  to  fluorescence  and  Ray¬ 
leigh  scattering  from  gases  and  particles.  A  profile 
of  S(,,.(488  nm)  is  then  obtained  (Fig  3a)  and  this 
profile  is  subtracted  from  the  S',, (488  nm)  profile  to 
yield  a  profile  that  is  only  due  to  elastic  scattering 
(Fig.  3b).  The  peak  in  Fig.  3b  is  due  to  scattering 
from  particles  while  the  balance  is  due  to  scattering 
from  molecules.  The  signal  rises  alter  the  peak  be¬ 
cause  we  are  entering  the  cold  fuel  zone.  The  crit¬ 
ical  strain  rate  for  the  suppression  of  soot  is  ob- 


Distance  Irom  Stagnation  Point  (mm) 
(a) 


Distance  from  Stagnation  Point  (mm) 

(t>) 

Fig.  3.  Profiles  of  the  measured  scattering  inten¬ 
sity,  in  arbitrary  units,  illustrating  the  method  of 
determining  the  particle-inception  limit.  In  Fig.  3a 
the  horizontally  and  vertically  polarized  compo¬ 
nents  of  detected  light  are  plotted  as  functions  of 
distance  from  the  stagnation  point.  The  detected 
light  has  passed  through  a  1  nm  laser  line  filter 
centered  on  the  incident  wavelength  488  nm  ami 
some  fluorescence  (S*,,)  is  observed  to  pass  through 
the  filter.  In  Fig.  3b  the  fluorescence  signal  has  been 
subtracted  from  SVf  so  that  the  remaining  signal  is 
elastic  scattering  from  molecules  and  soot  particles. 
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tamed  by  increasing  the  strain  rate  until  the  peak 
is  eliminated.  This  condition  is  identified  as  the 
particle-inception  limit  Kp  which  is  quite  sensitive 
to  changes  in  strain  rate  when  K  is  less  than  Kp, 
and  the  limit  can  easily  be  determined  within  around 
±3  s  For  K  greater  than  Kp  the  scattering  pro¬ 
file  is  only  weakly  affected  because  the  remaining 
signal  is  due  to  scattering  by  molecules. 

It  is  worthwhile  to  note  that,  as  with  follow 
flames,'"  the  peak  of  the  fluorescence  signal  just 
precedes  the  particle  inception  region;  this  region 
can  be  identified  by  the  initial  rise  of  the  peak  lo¬ 
cated  around  1.5  mm  from  the  luminous  zone.  This 
sequential  relationship  has  been  used  to  support  the 
view  that  the  fluorescing  molecules  are,  if  nothing 
else,  by-products  of  the  soot  inception  chemistry. 

It  is  desirable  to  base  the  particle  and  fluores- 
ence  limits  on  the  same  background  signal,  i.e.,  the 
molecular  scattering  signal.  To  accomplish  this  we 
continue  to  measure  SM  but  the  laser  line  filter  is 
removed  so  that  a  large  fluorescence  signal  is  ob¬ 
tained  At  this  strain  rate,  scattering  from  particles 
is  negligible  because  we  are  already  at  the  incep¬ 
tion  limit.  By  again  increasing  the  strain  rate,  the 
base-line  light  scattering  signal  from  molecules  is 
retrieved  at  another  critical  strain  rate,  Kp,  which 
we  shall  refer  to  as  the  fluorescence  limit  at  which 
production  of  the  fluorescing  PAH  s  is  assumed  to 
be  suppressed.  This  limit  might  be  interpreted  as 
a  representative  precursor  limit  although  it  is  rec¬ 
ognized  that  there  are  important  precursors  at  dif¬ 
ferent  stages  of  the  soot  formation  process  that  do 
not  fluoresce,  and  that  the  fluorescing  species  for 
=  -IBS  nm  may  not  actually  be  the  precursors 
themselves. 

The  sooting  limit  based  on  soot  luminosity  is  also 
measured  Here  the  limit,  indicated  by  Kp.  is  de¬ 
fined  as  the  condition  where  the  yellow  luminous 
zone  is  no  longer  visibly  perceptible  when  viewed 
through  a  5X  stereoscope. 

The  flame  extinction  limit,  Kk,  has  also  been  de¬ 
termined  by  increasing  the  strain  rate  until  the  flame 
is  extinguished.  We  then  measure  the  velocity  gra¬ 
dient  in  a  flame  with  a  slightly  lower  strain  rate 
where  a  stable  flame  can  exist  and  define  this  point 
as  the  flame  extinction  limit. 

Table  I  lists  the  limits  of  luminosity  (K(J,  particle 
inception  (Kp),  fluorescence  (Kp),  and  flame  extinc¬ 
tion  (Kp)  for  flames  of  methane,  ethane,  propane 
and  butane  with  air.  For  a  given  fuel/air  system, 
the  relation  Kp  =  Kp  <  Kp  <  Kp  is  evident.  Lu¬ 
minosity  appears  to  be  an  adequate  ir  thod  of  de¬ 
termining  the  particle  limit.  However,  it  should  be 
noted  that  the  luminosity  measurement  is  some¬ 
what  subjective  and  potential  errors  are  much  greater 
than  lor  the  measurement  of  Kp.  In  fact,  an  ac¬ 
ceptable  measurement  of  Kp  for  methane  could  not 
be  made  The  methane  flame  lias  a  noticeable 
amount  ol  emission  in  the  red  due  to  the  overtone 
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bands  ol  water16  and  this  emission  also  occurs  in 
the  region  of  the  luminous  yellow  zone.  Therefore, 
it  is  extremely  difficult  to  determine  the  limit  at 
which  the  luminous  yellow  has  disappeared  when 
the  background  is  red.  The  red  emission  was  ev  i¬ 
dent  to  a  lesser  extent  in  the  other  flames  and  did 
not  cause  as  great  an  ambiguity. 

Assuming  that  the  fluorescing  species  are  indic¬ 
ative  of  soot  precursors,  the  result  of  Kp  <  Kp  sug¬ 
gests  that  a  substantial  amount  of  precursors  may 
still  exist  even  though  the  production  ol  soot  par¬ 
ticles  has  been  suppressed.  It  is  reasonable  to  an¬ 
ticipate  that,  in  more  complex  flame  configurations 
these  precursors  could  survive  to  other  parts  ol  the 
flame  where  the  residence  time  and  temperature 
are  more  conducive  to  particle  generation  Thus  at¬ 
tainment  of  Kp  is  not  necessarily  sufficient  to  guar¬ 
antee  the  suppression  of  soot  generation. 

The  result  Kp  <  Kp,  is  interesting  in  that  it  in¬ 
dicates  the  possibility  of  completely  eliminating  soot 
generation  by  a  sufficiently  high  strain  rate  without 
extinguishing  the  flame.  We  hasten  to  add,  how¬ 
ever,  that  this  result  is  obtained  only  tor  the  pres¬ 
ent  relatively  non-sootv  fuels  of  the  lower  alkanes 
It  is  ol  interest  to  explore  whether  Kp  — *  Kp  tor 
more  sooty  fuels  such  that  total  soot  suppression 
through  straining  cannot  be  achieved  without  ex¬ 
tinguishing  the  flame.  Obviously  il  the  fuel  is  a 
highly -fluorescing  PAH,  then  Kp  -  Kp  would  lie 
expected  to  hold. 

Effects  of  Concentration  and  Temperature- 

While  soot  growth  in  diffusion  flames  has  been 
studied  in  detail,16'12  soot  inception  has  not.  F,.\- 
perimcntal  difficulties  are  largely  rcs|x>nsihic  for  this 
lack  of  information  about  inception.  However,  as 

*  .  .  n  i  .  .  i  i 

noted  by  Kennedy, the  critical  strain  rate  char¬ 
acterizes  the  inception  chemistry,  therefore  a  study 
of  how  flame  temperature  and  fuel  concentration 
influence  the  sooting  limits  of  diffusion  flames  will 
yield  information  about  how  these  parameters  in¬ 
fluence  inception  chemistry 

Figure  I  shows  the  various  sooting  limits  for  pro- 
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Fic.  -4.  Sooting  limits  based  on  soot  luminosity  (KJ, 
soot  particle  light-scattering  (KP)  and  PAH  fluores¬ 
cence  (Kf)  for  propane  diluted  with  nitrogen. 


pane  diluted  with  nitrogen.  These  results  again 
demonstrate  the  differences  in  the  limits  based  on 
the  various  definitions.  Furthermore,  the  results 
show  that  the  critical  strain  rates  K|_,  Kj.  and  Kp 
all  decrease  with  inert  addition  and  the  relative 
trends  ol  the  decrease  are  similar  for  all  three  def¬ 
initions.  In  fact,  the  curve  for  KL  is  similar  to  that 
for  Kp  even  though  inert  addition  decreases  the 
temperature  (e  g.,  for  X/  =  0.4,  AT  =  113  K)  and 
this  decrease  will  reduce  the  luminous  intensity. 

Since  the  chemistry  of  soot  particle  inception  is 
generally  a  function  of  both  temperature  and  re¬ 
active  species  concentrations,  the  results  of  Fig.  4 
do  not  directly  delineate  whether  the  reduction  in 
the  critical  strain  rates  with  fuel  concentration  is  a 
consequence  of  temperature  reduction,  concentra¬ 
tion  reduction,  or  both  To  isolate  these  two  ef¬ 
fects,  we  adopt  the  technique  of  Ref.  6.  That  is, 
after  dilution  of  the  fuel  by  a  certain  amount  of  in¬ 
ert,  we  bring  the  flame  temperature  back  to  the 
nondiluted  value  by  replacing  an  appropriate  molar 
amount  of  nitrogen  in  the  air  with  an  equal  amount 
of  argon. 

Figures  5  and  6  show  the  critical  strain  rates  for 
the  diluted  flames  with  and  without  adjustment  of 
the  flame  temperature.  The  difference  between  the 
nondiluted  value  and  the  diluted  value  with  tem¬ 
perature  adjustment  is  viewed  as  the  effect  due  to 
"pure  dilution,"  while  the  difference  between  the 
diluted  values  with  and  without  temperature  ad¬ 
justment  is  viewed  as  the  effect  due  to  flame  tem¬ 
perature  variation  The  curve  for  Kp  of  Fig  4  has 
been  duplicated  in  Fig.  5  to  illustrate  this  tech¬ 
nique. 

Referring  to  Fig.  5,  it  is  seen  that  the  effect  due 
to  pure  dilution  is  not  large  for  small  dilutions  but, 
for  fuel  mole  fractions  less  than  0.6,  the  effect  can 


be  as  large  as  that  due  to  llame  temperature.  In 
Ref.  6  we  have  demonstrated  a  similar  influence  of 
inert  addition  on  soot  volume  fraction  for  flames  with 
low  strain  rates  and  high  soot  loadings  It  was 
suggested1’  that  the  limited  influence  of  flame  tem¬ 
perature  when  an  inert  is  added  to  a  fuel  is  a  con¬ 
sequence  of  stoichiometry  and  the  presence  of  the 
large  amount  of  nitrogen  in  air  These  facts  imply 
that  the  flame  temperature  is  only  minimally  re¬ 
duced  even  with  a  moderate  amount  of  inert  ad¬ 
dition  to  the  fuel.  Reference  6  supports  this  expla¬ 
nation  with  calculated  adiabatic  flame  temperatures 
for  different  amounts  of  inert  dilution.  In  the  pres¬ 
ent  study  we  further  substantiated  these  calculated 
results  by  experimentally  determining  the  depen¬ 
dence  of  flame  temperature  on  inert  dilution. 

The  results  of  Fig.  6b  show  that  the  critical  strain 
rates  for  methane  are  noticeably  more  dependent 
upon  inert  addition  than  those  for  either  propane 
or  butane.  In  Ref.  6  we  showed  that  inert  addition 
influences  the  flame  temperature  of  methane  more 
than  it  does  that  of  propane  or  butane,  so  inert  ad¬ 
dition  (filled  symbols)  would  be  expected  to  have  a 
greater  influence  on  the  critical  strain  rates  of 
methane.  Figure  6b  suggests  that  even  when  the 
temperature  is  not  allowed  to  vary  with  inert  ad¬ 
dition.  i.e.  pure  dilution  (open  symbols/,  the  criti¬ 
cal  strain  rates  are  more  dependent  on  fuel  con¬ 
centration  for  methane  than  for  the  other  fuels.  By 
comparing  Kp  in  Figs.  5  and  6a,  it  is  also  seen  that 


FlC.  5  The  effect  of  nitrogen  addition  on  the 
particle-inception  limit  (KP)  for  propane.  The  open 
symbols  represent  data  tor  flames  where  the  tem¬ 
perature  ol  the  diluted  llaine  has  been  adjusted  to 
be  equal  to  that  of  the  undiluted  flame,  so  that  these 
data  represent  pure  dilution  effects. 
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Flo.  6  Tlu-  effect  of  nitrogen  addition  on  a)  the 
soot  luminosity  and  particle-inception  limit  for  bu¬ 
tane  and  b)  the  PAH  fluorescence  limit  for  methane 
and  butane.  The  open  symbols  represent  data  for 
flames  where  the  temperature  of  the  diluted  flame 
has  been  adjusted  to  be  equal  to  that  of  the  un¬ 
diluted  flame,  so  that  these  data  represent  pure  di¬ 
lution  effects. 


the  inception  limits  for  propane  and  butane  display 
similar  dependencies  on  fuel  concentration. 

Figures  5  and  6  show  the  influence  of  pure  di¬ 
lution  on  the  critical  strain  rates.  In  Fig.  T  we  pre¬ 
sent  the  critical  strain  rates  for  flames  with  a  fixed 
concentration  of  fuel  and  oxidizer  but  with  different 
flame  temperatures.  The  fuel  mole  fractions  for 
ethane,  propane  and  butane  are  1.0,  0.5  and  0.4, 
respectively.  The  flame  temperature  is  increased  by 
substituting  a  portion  of  the  nitrogen  in  air  by  an 
equal  amount  of  argon.  When  the  temperature  of 
the  flame  is  increased,  it  can  exceed  the  operating 


limits  of  the  thermocouple  Fortunately,  experi¬ 
ments  as  well  as  adiabatic  flame  temperature  cal¬ 
culations  show  that  the  ratio  of  AT  to  A  i %  An  is 
almost  a  constant,  and  this  constant  is  almost  in¬ 
dependent  ol  fuel  concentration.  Therefore  this  ra¬ 
tio  is  determined  at  lower  fuel  concentrations  where 
temperatures  are  measurable,  and  is  then  used  to 
determine  AT  when  temperatures  are  too  high  to 
measure. 

Figure  7  shows  that  when  concentration  effects 
are  eliminated,  the  flame  temperature  can  have  a 
strong  influence  on  soot  extinction  by  straining.  A 
200  to  300  K  increase  in  the  flame  temperature  can 
cause  a  corresponding  increase  in  the  critical  strain 
rates  by  factor  of  two  to  three.  Figure  8  shows  the 
Arrhenius  plots  of  the  critical  strain  rates.  A  linear 
behavior  is  observed,  implying  that  the  various  ex¬ 
tinction  states  can  be  quantified  by  a  corresponding 
extinction  Damkohler  number  Da  which  varies  with 
K"1  exp  t  — E„/R"T/),  where  F,„  is  an  overall  acti¬ 
vation  energy  and  TV  the  flame  temperature  The 
fact  that  an  Arrhenius  correlation  exists  lor  T f  is  in 
agreement  with  the  expectation  that  the  flame  tem¬ 
perature  is  useful  in  characterizing  global  temper¬ 
ature  effects. 

The  overall  activation  energies  for  particle  and 
fluorescence  extinction  are  shown  in  Fig.  S  for  the 
various  fuels.  It  is  of  interest  to  note  that  the  tem¬ 
perature  dependence  of  KP.  depicted  by  F„.  is  sim¬ 
ilar  for  all  of  the  alkanes  studied.  In  addition,  the 
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FtC.  7.  The  effect  ot  temperature  on  the  sooting 
limits  tor  the  various  fuels.  Temperature  adjust¬ 
ment  is  attained  without  changing  the  concentra¬ 
tions  of  either  the  fuel  or  oxidizer  The  tempera¬ 
tures  have  been  corrected  for  radiation  loss 
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Fit.'.  8.  Arrhenius-type  plot  for  the  sooting  limits 
for  various  fuels,  including  the  global  activation 
energies  obtained  from  least-squares  fit  to  the  data. 


values  of  E„  for  the  particle  and  fluorescence  limits 
of  butane  are  quite  close.  This  is  not  surprising  if 
the  chemistry  of  the  fluorescing  species  is  related 
to  that  of  particle  formation. 

Concluding  Remarks 

In  the  present  investigation  we  have  distin¬ 
guished  the  sooting  limits  based  on  luminosity,  light 
scattering  by  particles,  and  fluorescence.  Each  of 
these  extinction  phenomena  describes  a  different 
feature  of  the  soot  particle  formation  process  as 
characterized  by  the  respective  critical  strain  rates. 
The  luminous  limit  is  convenient  and  a  relatively 
good  indicator  of  the  absence  of  soot  particles  al¬ 
though  the  results  are  strongly  dependent  on  back¬ 
ground  conditions  which  can  vary  greatly  with  fuel 
type  and  flame  geometry.  The  inception  limit  based 
on  light  scattering  is  an  accurate  indicator  of  the 
absence  of  soot  particles,  yet  it  does  not  necessarily 
imply  suppression  of  the  soot  formation  process  be¬ 
cause  gaseous  soot  precursor  species  could  still  be 
present.  It  is  therefore  believed  that  a  sufficiently 
rigorous  indication  that  the  soot  formation  process 
is  totally  suppressed  requires  the  extinction  of  the 
soot  precursors.  The  fluorescence  limit  may  be  a 
suitable  method  of  depicting  this  condition,  al¬ 
though  further  work  must  be  performed  to  deter¬ 
mine  if  the  species  that  fluoresce  when  excited  at 
488  nm  are  indeed  critical  to,  or  at  least  are  by¬ 


products  of,  the  chemical  processes  leading  to  soot 
particle  inception. 

Results  further  show  that  soot  particle  formation 
is  sensitive  to  both  fuel  dilution  and  flame  temper¬ 
ature.  It  is  also  demonstrated  that  flame  tempera¬ 
ture  influences  both  the  fluorescence  and  particle- 
inception  limit  in  an  Arrhenius  manner. 

Finally,  the  finding  that  the  parameters  impor¬ 
tant  for  soot  inception  are  similar  to  those  for  soot 
formation  and  growth  is  in  agreement  with  the  well- 
accepted  observation  that  the  final  soot  loading  in 
sooty  flames  is  heavily  dependent  on  the  particle 
inception  process. 
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